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This  work  presents  a  full-Coulomb,  multi-electron  formalism  of  line  broadening 
due  to  perturbation  by  plasma  electrons.  This  work  is  based  on  the  relaxation 
theory  of  Smith  and  Hooper,  which  was  extended  to  treat  highly  ionized  radia- 
tors by  Tighe  and  Hooper.  A  further  extension  by  Woltz  and  Hooper  introduced  a 
full-Coulomb  expression  of  the  perturber-radiator  interaction  in  the  electron  broad- 
ening operator,  rather  than  applying  the  dipole  approximation.  The  term  of  the 
electron  broadening  operator,  first  order  in  the  radiator-perturbing  electron  inter- 
action is  reintroduced,  causing  a  red  line  shift  and  asymmetries.  This  formalism 
is  also  generalized  to  treat  multi-electron  radiating  ions,  using  the  formalism  to 
calculate  general  transition  arrays  presented  by  Cowan.  The  effects  of  the  plasma 
ions  are  treated  dynamically  using  an  ion  microfield  probability  distribution  with 
the  model  introduced  by  Boercker,  Iglesias,  and  Dufty. 


The  line  shapes  calculated  using  this  model  are  dramatically  shifted  and  dis- 
torted with  respect  to  lineshapes  calculated  in  the  dipole  approximation  of  the 
radiator-perturbing  electron  interaction.  They  also  are  narrower  and  the  differ- 
ences between  the  two  models  increase  as  the  principal  quantum  number  of  the 
upper  state  of  the  transition  increases.  This  can  be  attributed  to  the  fact  that  the 
orbital  size  is  greater  relative  to  the  average  plasma  electron  spacing.  The  asym- 
metries are  attributed  to  the  dependence  of  the  level  shifts  on  angular  momentum 
state.  There  is  good  agreement  between  the  shift  calculations  compared  with  those 
of  Griem  et  al.  and  Nguyen  et  al.  The  line  profiles  are  also  compared  to  those 
calculated  using  a  semi-classical  model  that  is  carried  to  all  orders  in  the  radiator- 
perturbing  electron  interaction.  Analysis  of  K-shell  Ar  spectra  using  the  shifted 
line  spectra  is  shown.  Finally,  the  impact  of  the  line  shift  on  merging  of  adjacent 
lines  in  a  Rydberg  series  at  ultra-high  densities,  <  3  x  10^'*cm~^,  is  considered. 


vi 


CHAPTER  1 
INTRODUCTION 

1.1    Plasma  Line  Broadening:  A  Brief  Survey 

The  study  of  plasma  line  broadening  has  been  employed  for  years  to  characterize 
hot,  dense  plasmas  [1,2].  As  a  radiator  interacts  with  a  plasma,  the  energy  levels 
of  its  bound  states  are  modified  by  interactions  with  the  surrounding  ions  and 
electrons,  resulting  in  broadened  line  spectra  that  is  dependent  on  the  plasma 
density  and  temperature.  Since  the  1970's,  high-powered  lasers  have  been  used 
to  implode  gas-filled  micro-balloons  [3] .  Usually  the  gas  fill  consists  of  deuterium 
and/or  tritium  doped  with  trace  amounts  of  elements  with  greater  nuclear  charge 
Z.  As  technology  improved,  progressively  larger  and  more  powerful  lasers  have 
become  available,  and  this  has  led  to  the  creation  of  plasmas  with  ever  increasing 
densities  and  temperatures  [4,  5,  6,  7]. 

The  understanding  of  these  plasmas  has  necessitated  continuing  advances  in 
theoretical  plasma  line  broadening,  beginning  with  Stark  broadening  of  neutral  ra- 
diators due  to  the  ion  microfield  as  well  as  broadening  due  to  plasma  electrons  [8,  9]. 
The  line  broadening  formalism  was  extended  for  ionic  radiators  when  temperatures 
and  densities  become  too  high  to  study  spectra  from  neutral  radiators  [10,  11,  12], 
as  well  as  full-Coulomb  calculations  when  the  dipole  approximation  was  no  longer 
valid,  again  at  high  densities  [13,  14].  Then,  with  the  availability  of  the  Nova-  and 
24-beam  Omega  lasers,  it  was  possible  not  only  to  see  K-shell  Ar  lines  but  the 


1 


2 

satellites  of  the  resonance  lines  as  well.  At  the  same  time  it  was  possible  to  ob- 
serve broadened  L-shell  lines  from  krypton.  This  was  propitious  since  we  had  the 
theoretical  capability  to  calculate  the  spectra  from  multi-electron  radiators  with 
the  computer  code  MERL  [15,  16].  Multi-electron  calculations  could  treat  L-shell 
lines  [17]  and  also  lines  from  several  ionization  stages  [5,  6,  18]. 

As  plasma  electron  densities  on  the  order  of  1  x  lO^^cm"^  were  reached,  the 
controversial  Plasma  Polarization  Shift  was  unambiguously  observed  and  diag- 
nosed [19,  20,  21,  22].  The  radiators  were  so  highly  ionized  and  the  plasma  electron 
densities  were  so  great  that  the  electrons  penetrating  radiator  orbitals  could  shield 
the  nucleus  from  the  radiator  electron,  thus  lowering  its  energy  state.  Although 
spectral  line  broadening  due  to  plasma  ions  and  electrons  has  been  observed  for 
many  years  and  has  been  shown  to  be  consistent  with  existing  theories  at  lower 
plasma  densities  [23,  24],  the  existence  of  plasma-induced  line  shifts  has  only  re- 
cently been  widely  accepted. 

As  early  as  twenty  years  ago  [1,  25],  interest  in  spectra  of  ionized  radiators 
in  dense  plasmas  motivated  ad  hoc  calculations  of  red  shifts  induced  by  perturb- 
ing electron  shielding  of  the  nucleus.  Early  attempts  to  calculate  the  shift  em- 
ployed linearized  Debye-Huckel  potentials,  while  others  used  self-consistent  static 
charge  distributions  to  describe  the  shielding  of  the  nucleus  [24,  26,  27,  28,  29]. 
These  calculations  of  what  was  called  the  Plasma  Polarization  Shift  (PPS)  fre- 
quently overestimated  the  shifts  observed  experimentally,  particularly  those  using 
the  Debye-Huckel  model.  This  shift  was  often  written  as  a  first  order  average  inter- 
action between  the  radiator  and  the  plasma  electrons,  and  Cooper,  Kelleher  and 
Lee  [30]  cautioned  against  introducing  an  ad  hoc  shift  term  in  conjunction  with 
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the  electron  broadening  operator,  pointing  out  that  a  shift  term  should  naturally 
arise  as  the  real  part  of  the  broadening  operator,  while  the  imaginary  part  serves 
to  broaden  the  line-shape.  Boercker  and  Iglesias  [31]  had  stated  that,  for  neutral 
radiators,  the  first  order  average  interaction  was  negligible  and  showed  that  there 
was  significant  cancelation  in  the  terms  second  order  in  the  interaction  potential. 
This  result  was  sometimes  misinterpreted  to  mean  that  the  entire  shift  experienced 
substantial  cancelation.  This  misunderstanding  and  the  poor  agreement  of  early 
theories  with  the  little  available  data  led  some  to  meet  any  line  shift  calculations 
with  skepticism. 

Ultimately,  line  shift  calculations  for  highly  ionized,  moderate  Z  radiators  that 
employed  scattering  theory  first  order  in  the  radiator-perturbing  electron  interac- 
tion provided  results  that  agreed  well  with  each  other  and  eventually  with  experi- 
ment [23,  24].  In  order  to  obtain  a  consistent  result  for  both  line  broadening  and 
a  line  shift,  Nguyen  tt  al.  used  quantum  mechanical  impact  theory  to  calculate 
shift  and  width  terms  for  Lyman  series  members  of  high-Z  radiators  ~  10  [24]. 
The  calculation  of  Nguyen  et  al.  included  all  terms  of  the  multipole  expansion  of 
the  radiator-perturbing  electron  interaction,  up  to  the  octopole.  These  researchers 
discovered  that  the  monopole  term  of  the  interaction  in  which  the  plasma  electron 
penetrates  the  radiator  orbital  contributed  the  most  to  the  shift.  This  is  significant 
because  the  electron  broadening  operator  was  usually  calculated  in  the  dipole  ap- 
proximation, which,  although  valid  for  higher  order  terms,  allowed  the  first  order 
shift  term  to  reduce  to  zero  erroneously.  Griem  et  al.  [23]  applied  a  distorted  wave 
scattering  approach  to  the  problem,  employing  perturber  wavefunctions  distorted 
by  only  the  monopole  term  of  the  radiator-perturbing  electron  interaction.  This 


first-order  impact  calculation  summed  the  phase  shifts  of  the  partial  waves  to  ob- 
tain the  line  shifts  for  He-like  Ar.  The  results  were  in  good  agreement  (within 
~  10%)  with  Nguyen  et  al.  [24]  and  the  calculation  demonstrated  that  the  bulk  of 
the  shift  could  be  obtained  by  including  only  the  monopole  term  of  the  radiator- 
perturbing  electron  interaction  potential. 

However,  when  researchers  attempted  to  present  shift  measurements,  many 
objections  were  raised.  Only  recently  has  a  consensus  been  reached  as  line  shift 
observations  confirmed  the  calculations  by  Nguyen  et  al.  [24]  and  Griem  et  al  [23]. 
One  obstacle  to  an  early  experimental  verification  of  shifts  of  spectra  emitted  from 
highly  charged  radiators  was  the  attention  paid  to  K-shell  a-lines,  because  we  now 
know  a-lines  experience  very  small  shifts  relative  to  higher  lying  series  members 
and  which  were  also  most  likely  to  be  obscured  by  opacity  effects.  In  the  following 
paragraphs  we  will  outline  two  important  experiments  where  shifts  were  observed 
and  the  concerns  that  were  raised. 

Hammel  et  ai,  [5]  measured  a  red  shift  in  the  He-/3  line  of  Ar,  approximately 
lleV  at  an  electron  density  of  approximately  1.2  x  10^^  cm~^,  which  was  roughly 
consistent  with  the  theories  of  Nguyen  et  al.  and  Griem  et  al.  and  was  greater  than 
the  instrumental  width  of  their  high-resolution  spectrometer.  However,  the  pres- 
ence of  dielectronic  satellites  on  the  red  wing  complicated  the  shift  measurement 
and  obscured  asymmetries  arising  from  the  differing  shifts  of  the  various  angular 
momentum  states. 

Leng  et  al.,  [32]  observed  small  red  shifts  and  a  marked  asymmetry  of  the  Ly-7 
of  CVI.  The  bulk  shifts  of  the  lines  seemed  to  be  comparable  with  those  calcu- 
lated by  Nguyen,  However,  the  theoretical  line-shapes  were  calculated  using  an 


average  shift,  as  opposed  to  incorporating  shifts  dependent  on  the  angular  mo- 
mentum states  of  the  radiator  into  the  initial  line-shape  calculation.  Therefore, 
Leng  et  al.  employed  theoretical  line-shapes  that  neglected  the  asymmetries  due 
to  these  differing  shifts.  In  addition,  existing  plasma  density  gradients  could  have 
obscured  details  associated  with  line  shifts,  particularly  shift  asymmetries.  Further 
objections  arose  because  the  lines  for  different  densities  were  convolved  together 
to  compensate  for  plasma  gradients,  which  would  have  interfered  further  with  a 
definitive  shift  measurement.  Finally,  relying  heavily  on  the  correctness  of  com- 
puter modeling  results  for  the  density  inferences  from  the  widths  of  the  lines  cast 
doubt  on  the  application  of  any  particular  shift  theory. 

In  the  past  few  years,  more  evidence  of  plasma-induced  line  shifts  has  begun 
to  appear.  First,  Renner  et  al.  observed  the  Lyman  series  lines  of  aluminum  shift 
to  the  red  during  high  resolution  flat  target  experiments  [19].  In  this  case,  the 
density  inferences  from  the  widths  of  the  lines  relied  heavily  on  the  correctness 
of  computer  modeling  results  and  cast  doubt  on  the  application  of  any  particular 
shift  theory.  Second,  at  Lawrence  Livermore  National  Laboratory  and  the  Labora- 
tory for  Laser  Energetics,  in  implosion  experiments  conducted  on  deuterium  filled 
microballoons,  doped  with  trace  amounts  of  argon,  an  apparent  red  shift  of  the 
He-^  and  He-7/Lyman-/?  complex  of  K-shell  Ar  was  observed,  even  after  satellites 
were  included.  [18]  Since  some  of  these  shifts  were  far  in  excess  of  the  experimental 
error,  the  University  of  Florida  group  incorporated  a  first-order  shift  calculation 
for  highly  ionized  radiators  into  our  line-shape  code,  MERL  [20].  Line  shifts  due  to 
plasma  electron  interactions  with  a  radiator  had  not  been  incorporated  into  multi- 
electron  line-shape  codes  such  as  MERL  since  the  existence  of  line-shifts  had  not 
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been  conclusively  observed  in  experiments.  More  recently,  Woolsey  et  al.  and 
Saemann  et  al.,  reported  red  shifts  in  the  He-/3  of  K-shell  argon  and  aluminum, 
respectively,  both  roughly  consistent  with  Nguyen's  results.  [21,  22] 

1.2    The  Focus  of  this  Work 

We  will  extend  the  theoretical  foundation  of  plasma  line  broadening  by  pre- 
senting a  generalized  full-Coulomb  formalism  of  line  broadening  for  multi-electron 
ions,  where  the  line  shifts  arise  consistently  from  the  relaxation  theory. 

There  are  several  reasons  to  study  the  dense  plasma  line  shift,  despite  the 
controversy.  First,  since  accurate  plasma  diagnostics  has  been  the  primary  goal  of 
plasma  line  broadening,  it  is  significant  if  we  can  provide  more  accurate  descriptions 
of  analyzed  spectra.  Secondly,  we  will  have  a  more  complete  understanding  of 
plasma-induced  modifications  of  energy  level  structure  of  atoms  and  ions.  This  is 
not  a  purely  academic  pursuit:  indeed,  obtaining  a  complete  model  of  the  plasma- 
induced  shift  provides  a  basis  from  which  to  understand  the  effect  of  the  shift  on 
plasma  phenomenon  including  continuum  lowering  and  ionization  balance  [33,  34, 
35,  36,  37],  which  in  turn  impact  opacities  and  line-ratio  diagnostics. 

Despite  the  success  of  Griem  et  al.  [23]  and  Nguyen  et  al.  [24]  improvements  on 
their  work  were  possible.  Neither  presented  a  dynamic,  in  other  words  frequency- 
dependent,  shift  and  width  operator.  Also,  both  calculations  were  effectively  lim- 
ited to  first  order  expressions  of  the  the  radiator-perturber  electron  interaction. 
No  attempt  was  made  by  Nguyen  et  al.  to  calculate  line-shapes  from  the  results. 
Finally,  their  work  was  limited  to  K-shell  radiators.  Multi-electron  models  are 
necessary  for  spectral-line  analysis,  since  composite  lines,  that  is,  resonance  lines 


together  with  their  accompanying  satellites,  are  an  effective  diagnostic  of  temper- 
ature as  well  as  electron  density. 

This  formalism  extends  their  work  by  enabling  one  to  calculate  the  shift  and 
width  terms  simultaneously  and  self-consistently  using  the  relaxation  model  [9, 
38].  The  full-Coulomb  expression  of  the  radiator-perturbing  electron  interaction 
is  necessary  not  only  because  the  dipole  approximation  of  the  radiator-perturbing 
electron  interaction  breaks  down  at  the  high  densities  of  interest,  but  also  because 
the  shift  is  dominated  by  the  monopole  term  of  the  multipole  expansion  of  this 
interaction,  which  is  lost  if  only  the  dipole  approximation  is  employed.  Finally, 
this  work  generalizes  the  shift  and  width  operator  to  treat  arbitrary  multi-electron 
radiators  in  order  to  address  composite  line  spectra  as  well  as  L-shell  line  spectra. 

1.3  Outline 

In  Chapter  II,  a  theoretical  basis  for  spectral  line  broadening  will  be  presented. 
A  set  of  assumptions  will  be  put  forth  that  are  appropriate  for  the  plasma  condi- 
tions under  consideration.  This  will  allow  an  application  of  kinetic  theory  tech- 
niques to  obtain  a  form  of  the  electron  shift  and  width  operator,  as  described  by 
Hussey,  in  order  to  describe  the  effects  of  the  plasma  electrons  on  the  lineshape. 

Chapter  III  will  focus  on  the  electron  shift  and  width  operator.  We  will  intro- 
duce a  useful  form  of  the  radiator-perturbing  electron  interaction  for  multi-electron 
ions.  We  will  then  justify  the  application  of  the  second  order  approximation,  and 
a  generalization  for  multi-electron  radiators  will  be  presented.  Finally,  we  will 
obtain  calculable  forms  of  the  electron  shift  and  width  operator  by  addressing  the 
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issue  of  calculating  matrix  elements  of  the  interaction  for  complex  configurations 
which  experience  intermediate  angular-momentum  coupling. 

We  will  present  our  results  in  Chapter  IV.  The  behavior  of  the  shift  as  a  func- 
tion of  various  parameters  will  be  explored,  as  will  the  impact  of  the  level  shift 
on  line  shape  and  position.  The  shifts  obtained  from  this  formalism  will  be  com- 
pared to  the  theories  of  Nguyen  and  Griem,  and  lineshapes  calculated  using  the 
current  formalism  will  be  compared  to  lineshapes  generated  from  a  semi-classical 
formalism  that  includes  all  orders  in  the  radiator-perturbing  electron  interaction. 
Multi-electron  line  profiles  will  be  shown,  comparing  the  full-Coulomb  calculation 
with  some  simpler  models.  Experimental  data  will  be  fit  with  theoretical  spectra 
with  and  without  the  dense  plasma  line  shift.  The  effect  of  the  level  shift  on  the 
phenomenon  of  line  merging  will  also  be  demonstrated. 

In  the  Conclusion,  we  will  summarize  our  findings  and  consider  the  effects  of 
including  the  dense  plasma  line  shift  on  other  aspects  of  the  field  of  dense  plasma 
physics. 


CHAPTER  2 
PLASMA  LINE  BROADENING 

In  this  chapter,  we  will  outline  the  theory  of  plasma  line  broadening.  This 
work  is  based  on  a  dynamical  model  of  ion  broadening  and  a  relaxation  model  of 
electron  shifts  and  broadening.  We  will  first  discuss  the  sources  of  line  broadening 
with  shifts  in  dense  plasmas.  Then,  we  will  consider  the  time  scales  over  which 
interactions  between  the  radiator  and  the  plasma  take  place  in  order  to  determine 
the  appropriate  manner  in  which  the  plasma  effects  should  be  addressed.  Useful 
scaling  parameters  will  also  be  introduced.  Finally,  we  present  the  line  broad- 
ening formalism  including  Doppler  broadening  due  to  radiator  motion,  a  model 
of  electron  shifts  and  broadening  using  kinetic  theory  techniques  [38,  39],  and  a 
dynamical  model  of  ion  broadening  [1,  18,  40]. 

2.1    Model  of  the  Plasma 

In  order  to  gain  a  complete  understanding  of  all  the  interactions  that  a  highly 
ionized  radiator  experiences  in  a  plasma  and  how  these  interactions  might  affect 
line  spectra  emitted  from  such  radiators,  we  need  to  develop  a  model  of  the  plasma. 
As  we  indicated  in  the  introduction,  high-powered  lasers  have  been  used  to  im- 
plode microballoons  filled  primarily  with  deuterium  and/or  tritium  doped  with 
trace  amounts  of  elements  with  greater  nuclear  charge  Z.  The  experiments  under 
consideration  involve  trace  amounts  of  argon  or  krypton  in  deuterium  or  plasmas 
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of  pure  Ar  or  Kr,  or  mixtures  of  the  two.  At  electron  densities  from  1  x  10^^ /cc 
to  better  than  1  x  10^'' /cc  and  temperatures  of  approximately  1  keV  or  greater, 
line  spectra  from  K-shell  Ar  and  L-shell  Kr  can  be  expected.  At  even  higher 
temperatures,  K-shell  Kr  has  been  observed  [41]. 

At  such  high  temperatures,  the  deuterium  would  be  completely  ionized.  As  a 
result,  we  anticipate  that  a  radiator  interacted  with  fast  moving  plasma  electrons 
that  had  once  been  bound  to  the  deuterium  ions,  as  well  as  with  the  deuterium 
ions  themselves.  The  deuterium  ions,  by  virtue  of  their  greater  mass,  would  not 
be  expected  to  have  as  great  an  average  velocity  as  the  electrons,  and  the  Ar  or 
Kr  ions  would  move  still  more  slowly.  Since  the  gases  inside  the  microballoons 
undergo  spherically  symmetric  implosions,  we  can  assume  the  bulk  velocity  of  the 
plasma  is  zero. 

As  the  electron  density  increases,  we  can  anticipate  close  collisions  between  a 
radiator  and  plasma  electrons.  If  the  density  is  sufficiently  high,  we  can  expect 
an  increasing  probability  of  perturbing  electrons  penetrating  the  orbitals  of  bound 
radiator  electrons.  The  nature  of  this  interaction  will  necessitate  a  detailed  un- 
derstanding of  individual  radiator-plasma  electron  collisions.  On  the  other  hand, 
the  plasma  ions  would  be  repelled  as  strongly  as  the  electrons  are  attracted  to  the 
highly  charged  radiator. 

Also,  since  the  ions  move  more  slowly  than  the  electrons,  we  can  anticipate  that 
their  interactions  with  the  radiator  would  be  screened  by  the  plasma  electrons.  As 
a  result,  we  will  be  able  to  treat  the  effects  of  the  plasma  ions  collectively  in 
terms  of  a  fluctuating  electric  field  due  to  all  the  plasma  ions,  taking  into  account 
screening  due  to  plasma  electrons. 
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Next,  we  will  consider  the  mechanisms  by  which  spectral  lines  might  be  altered 
in  such  a  plasma. 

2.2    Sources  of  Line  Shifts  and  Broadening 

There  are  several  means  by  which  spectral  lines  are  broadened  and  shifted. 
First,  there  exist  "natural"  widths  of  spectral  lines  which  can  be  obtained  from 
the  rates  of  spontaneous  emission  between  initial  and  final  states  [1,  42].  However, 
for  the  plasma  conditions  under  consideration,  natural  line  broadening  provides  a 
neglible  contribution  to  the  total  line  width,  a  fraction  of  an  electron-volt  for  the 
Ly-P  of  Ar+^'',  for  example. 

The  Doppler  effect  is  another  source  of  broadening  and  shifts.  As  radiators 
move  within  the  plasma,  the  emitted  photons  are  Doppler  shifted.  If  there  is  only 
thermal  motion,  the  collection  of  photons  will  be  shifted  both  to  higher  and  lower 
energies,  thereby  broadening  the  line.  There  can  also  be  a  line  shift  corresponding 
to  bulk  motion  of  the  plasma.  Doppler  broadening  is  incorporated  into  the  current 
model;  however,  the  plasmas  of  interest  do  not  experience  much  bulk  motion  and 
do  not  necessitate  including  a  bulk-Doppler  shift. 

Another  source  of  broadening  arises  from  absorption  and  re-emission  of  pho- 
tons as  they  travel  through  the  plasma,  i.e.  opacity  broadening.  Photons  are 
preferentially  absorbed  at  frequencies  where  the  intensity  is  the  greatest,  but  can 
be  re-emitted  over  a  wide  range  of  angles  and  frequencies.  This  has  the  effect  of 
lowering  the  peak  of  a  spectral  line  and  distributing  the  intensity  to  the  wings. 
The  effects  of  opacity  broadening  will  be  included  in  the  analysis  of  experimental 
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data;  however  this  phenomenon  is  described  in  detail  in  other  works  and  will  not 
be  pursued  here  [1,  18,  43]. 

We  will  be  examining  plasma  regimes  that  are  sufficiently  dense  that  the  domi- 
nant source  of  line  broadening  is  due  to  the  Stark  effect.  The  Stark  effect  refers  to 
the  modification  of  the  energy  level  structure  of  the  bound  radiator  electrons  due 
to  an  external  electric  field.  Although  both  plasma  electrons  and  ions  contribute 
to  this  field,  they  cannot  be  treated  in  the  same  manner.  In  order  to  treat  each 
type  of  interaction  correctly,  one  needs  to  determine  the  time  scale  over  which  they 
act. 

2.3    Time  and  Length  Scales 

Whether  one  treats  the  interactions  of  the  radiator  with  the  plasma  with  a 
quasi-static  or  dynamical  model  is  dependent  on  the  time  scales  over  which  the 
interactions  occur,  relative  to  the  time  over  which  the  excited  ion  or  atom  radiates. 
A  radiative  lifetime  corresponds  to  the  frequency  half  width  of  a  lineshape  by 
using  a  relation  utilizing  the  Uncertainty  Principle,  AuAt  ~  1.  Portions  of  a 
lineshape  at  a  large  frequency  separation  from  line  center  correspond  to  transitions 
occurring  over  a  short  time.  If  a  perturber  acts  slowly  over  a  time  long  relative 
to  this  radiative  lifetime,  the  effect  of  the  perturber  on  the  radiator  does  not 
change  significantly  over  the  lifetime  and  may  be  treated  using  a  quasi-static  model. 
On  the  other  hand,  very  small  frequency  separations  from  line  center  correspond 
to  transitions  occurring  over  a  relatively  long  time.  If  the  time  over  which  a 
perturber  acts  is  short  relative  to  this  radiative  lifetime,  the  interaction  between 
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the  perturber  and  the  radiator  during  the  lifetime  is  time-dependent  and  the  effects 
of  the  perturber  must  be  treated  dynamically. 

We  will  now  determine  the  change  in  the  effect  of  the  perturbers  on  the  radiator 
for  an  average  radiative  lifetime  for  the  plasmas  under  consideration.  The  average 
radiative  lifetime  corresponds  to  the  experimentally  observed  full  width  at  half 
maximum  for  a  given  lineshape.  The  fwhm  of  a  K-shell  /?-line  is  approximately 
40  eV,  which  corresponds  to  At  =  10~^^  seconds.  We  will  determine  the  relative 
change  in  the  electric  field  of  a  perturber  changes  over  this  time,  6S/£, where  the 
electric  field  is  defined  as  follows: 


^^='^2  (2-1) 


where  r  is  the  distance  between  the  perturber  and  the  radiator.We  now  take  the 
differential  to  get  an  expression  for  S£, 


5r  5r 


Se  =  eS{r-^)  =  -2e—  =  -2—^.  (2.2) 


In  order  to  set  an  upper  bound  on  58 /S  we  set  6r  equal  to  v^ygAt.  An  aver- 
age velocity  of  a  perturbing  particle  of  a  given  species  can  be  obtained  from  the 
temperature  as  follows, 

/2A;rTA^/^ 

(2.3) 


where  ks  is  the  Boltzmann  factor,  Tj  is  the  temperature  of  that  species,  and  mj  is 
the  mass  of  the  species.  We  will  consider  the  case  of  highly  ionized  Ar  immersed 
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in  a  plasma  of  electrons  and  either  deuterium  or  argon  ions,  for  a  temperature, 
kT  =  lOOOeK  and  an  electron  density  of  1  x  lO^Vcc. 

Below,  we  show  the  values  of  S£/£  for  various  perturbers.  We  will  next  choose 
a  value  of  r  that  corresponds  to  the  distance  of  closest  approach. 

se        3A     ,  , 

<  —      for  electrons 
8  rc 

.05k  ^  , 

<    for  deuterium  ions 

rc 

.012A  ^ 

<    for  argon  ions  (2.4) 


In  order  to  determine  the  distance  of  closest  approach  for  the  plasma  electrons,  we 
begin  with  the  inter-electron  spacing,  roe, 

Toe  =  y—rie  j       ,  (2.5) 

where  Uc  is  the  electron  density.  This  gives  us  a  value  of  .62  A  for  r^e  at  an  electron 
density,  =  lO^V^c  •  The  relative  change  in  the  electric  field  becomes  5S/E  <  4.8. 
Since  the  radiators  under  consideration  are  highly  charged,  the  distance  of  closest 
approach  is  probably  even  smaller.  Therefore  we  can  assume  5£/£  will  be  even 
larger.  It  is  clear  that  we  must  treat  the  electrons  dynamically. 

If  we  consider  a  plasma  of  pure  argon,  we  must  employ  the  average  particle 
distance,  r^j  =  Z  roe,  where  Z  refers  to  the  average  charge  of  the  perturbing  ions. 
The  value  of  Z  will  be  on  the  order  of  16  for  the  plasmas  under  consideration.  This 
gives  us  a  value,  S£/S  <  .008,  which  suggests  that  it  is  appropriate  to  treat  argon 
ion  perturbers  quasi-statically. 
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If  the  radiator  is  immersed  in  a  deuterium  plasma,  the  deuterium  ions  will  be 
completely  ionized  and  their  density  will  roughly  correspond  to  the  electron  density. 
As  a  result  we  get  the  value  .08  for  bSfS,  which  suggests  employing  a  quasi-static 
model.  However,  if  we  consider  a  smaller  frequency  separation,  for  example  10 
eV  or  less,  the  relative  change  in  electric  field  begins  to  approach  one.  This  is 
consistent  with  calculations  of  argon  spectra  employing  a  dynamical  model  of  ion 
effects.  The  calculations  of  Ar  line  spectra  using  both  quasi-static  and  dynamical 
models  diflfered  significantly  only  near  line  center  [18].  Haynes  et  al.  considered 
plasma  mixtures  of  deuterium  and  argon  ions  and  found  that  as  the  concentration 
of  Ar  ions  was  increased,  the  lineshapes  calculated  in  the  dynamic  model  reduced 
to  the  static  limit,  with  little  difference  in  the  two  models  at  approximately  20% 
Ar  in  deuterium  [18]. 

Another  appropriate  time  scale  for  the  interactions  of  the  radiator  with  elec- 
trons would  be  the  time  an  electron  can  travel  close  to  an  ion  and  not  be  shielded 
by  the  intervening  plasma.  A  length  scale  appropriate  to  the  interaction  would  be 
the  Debye  length,  Xp,  because  it  represents  the  distance  past  which  a  charge  is 
effectively  screened  from  the  radiating  ion  by  the  plasma  environment. 


where  n  is  the  density  of  the  intervening  particles  and  e  is  the  charge  of  an  electron. 
We  will  attribute  most  of  the  screening  to  the  plasma  electrons. 

If  we  determine  the  time  it  takes  for  an  electron  to  travel  across  a  Debye 
sphere  when  it  has  the  most  probable  velocity  for  a  given  temperature,  we  obtain 


(2.6) 
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a  quantity  that  can  be  related  to  the  electron  plasma  frequency. 


1/2 


(2.7) 


1 


(2.8) 


Therefore,  the  electrons  must  be  treated  dynamically  when  the  frequency  sepa- 
ration is  less  than  the  electron  plasma  frequency,  which  coincidentally  is  approxi- 
mately 35  eV  for  an  electron  density,  Ue  =  10^'* /cc  .  For  all  the  spectral  lines  under 
consideration,  significant  portions  of  the  line  profile  lie  within  this  range. 

In  addition  to  considering  the  interactions  of  the  plasma  with  the  internal  states 
of  the  radiator,  we  must  consider  the  effect  of  the  plasma  on  the  center  of  mass 
motion  of  the  radiator.  We  will  be  making  the  approximation  that  the  velocity  of 
the  radiating  ion  will  be  constant  during  the  lifetime  of  the  transition:  therefore, 
we  must  determine  whether  this  approximation  is  valid.  We  can  safely  assume 
that  little  momentum  is  transfered  to  the  radiator  center  of  mass  as  it  collides 
with  plasma  electrons;  however,  we  will  estimate  the  momentum  transfered  to  the 
radiator  as  it  interacts  with  the  ion  microfield.  We  will  simply  use  the  definition 
of  impulse,  assuming  that  the  electric  field  can  be  approximated  by  the  field  due 
to  a  plasma  ion  at  the  ion  sphere  radius,  described  as  follows, 


respectively,  Zp  is  the  charge  of  the  perturbing  ion  and     is  the  distance  of  closest 


TUrAVr  =  FAt  < 


At, 


(2.9) 


where  m^,  Vr,  and  Zr  are  the  mass,  velocity,  and  charge  of  the  radiating  ion, 


approach. 
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We  will  compare  the  change  in  momentum  to  the  average  momentum,  deter- 
mined using  the  ion  temperature  in  the  following  expression: 

TTlrAVr  ZrZpC^  At 

rrirVr     ^       rl     {2mrkBTy^^  ^^'^^^ 

— — -   <   .004^,  2.11 
rrirVr  rj 

for  a  radiating  Ar+^^  ion  at  kT  =  1000  eV.  We  determine  this  quantity  for 
perturbing-deuterium  and  argon  ions  to  be  approximately  .003  and  .0075,  respec- 
tively, for  an  electron  density  of  rig  =  W^/cc,  which  suggests  that  we  can  treat  the 
center-of-mass  motion  independently  of  the  other  degrees  of  freedom  of  the  plasma 
system. 

There  is  one  more  quantity  we  will  consider  before  we  proceed  with  outlining 
the  general  theory.  The  ratio  of  the  radiator  orbital  size,  f,  to  the  average  inter- 
electron  spacing,  r^e,  taken  to  the  third  power,  is  a  measure  of  the  average  number 
of  plasma  electrons  'penetrating'  an  orbital  at  a  given  density.  This  number,  a, 
draws  attention  to  the  need  to  eliminate  the  dipole  approximation  of  the  radiator- 
perturbing  electron  interaction  at  high  electron  densities.  Also,  it  gives  a  strong 
indication  of  the  likelihood  that  one  will  observe  a  plasma  induced  shift.  The 
electron  density  per/cc  is  represented  by  n^,  refers  to  the  principle  quantum 
number,  /  represents  the  orbital  angular  momentum,  and  Z;  is  the  effective  charge 
of  the  radiating  ion.  We  define  a  as. 


a  — 


Jo,e)  ~     3     V  iz: 


(2.12) 


r 
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If  this  parameter  is  also  applied  to  ions,  one  might  draw  the  conclusion  that 
a  full-Coulomb  representation  of  the  radiator-perturbing  ion  interaction  is  also 
necessary,  especially  when  the  perturbers  are  deuterium  ions.  However,  we  must 
remember  that  the  radiator  has  a  positive  charge  and  would  repel  the  ions  as 
strongly  as  the  electrons  are  attracted.  We  do  not  need  to  be  concerned  with  close 
interactions  between  the  radiator  and  individual  perturbing  ions. 

2.4    The  Line  Shape  Expression 

In  this  section,  we  will  introduce  a  general  power  distribution  for  dipole  radi- 
ation. We  develop  the  expression  in  terms  of  the  autocorrelation  of  the  radiator 
dipole.  Since  the  power  distribution  will  be  written  in  terms  of  an  ensemble  av- 
erage over  plasma  states,  we  will  discuss  how  the  average  will  be  evaluated  over 
the  various  degrees  of  freedom  within  the  radiator-plasma  system.  We  will  then 
separate  out  the  effects  of  the  motion  of  the  radiator  center  of  mass,  which  leads 
to  Doppler  broadening. 

The  power  distribution  of  dipole  radiation,  as  a  function  of  frequency,  uj  can 
be  written  in  terms  of  an  ensemble  average.  The  derivation  of  this  expression  has 
been  described  in  many  other  sources  and  will  not  be  presented  here.[l,  2] 

=  ^  E  l(^|de-''^-|a)|Va^(a;  -  a;,,),  (2.13) 

ab 

where  the  radiator  dipole  is  represented  by  d,  the  indices  a  and  b  represent  the 
initial  and  final  states,  respectively,  of  the  radiator  together  with  the  plasma  en- 
vironment. The  position  of  the  center  of  mass  of  the  radiator  is  denoted  by  r. 
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and  the  vector,  k,  represents  the  wavevector  of  the  emitted  radiation,  uj/c.  pa 
represents  the  eigenvalue  of  the  density  matrix  for  the  initial  states  of  the  radiator 
and  the  plasma  environment.  The  energy  difference  of  the  initial  and  final  states, 
{Ea  —  Ei,)h  is  written  as  Uab- 

The  line  profile  can  be  rewritten  in  terms  of  its  transform  using  the  integral 
representation  of  a  Dirac  delta  function, 

^(^)   =    ^Pi^)  (2.14) 
=   $:(a|e^'^'-d|6)(6|de-"^'-pa|a)(5(a;-a;„6),  (2.15) 

ab 


1  r°° 

=   —  E  /     (ite^^'^-")'(a|e'''"-de'^'-'/'^|6)(6|de-''^-"e-'^"*/V|a)  (2.16) 

ab 
1  f°° 

-   —  /     dte'^''-''^'Tv{e'^  ''de'"'^^  ■de-'^  'e~'"'^''p},  (2.17) 

^TT  J — oo 


where  the  trace  over  the  plasma  radiator  system  is  denoted  by  Tr{-  •  •}. 

The  power  distribution  is  now  expressed  in  terms  of  in  time  development  oj)- 
erators,  which  are  derived  from  the  usual  equation  of  motion; 

ih-m  =  HT{t)  (2.18) 

which  has  the  solution 


T{t)  =  e-'^«/^      where     T(0)  =  1, 


(2.19) 
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and  H  is  the  Hamiltonian  of  the  combined  radiator  and  plasma  system.  The 
damping  factor  e  represents  the  natural  line  width;  although  it  is  a  small  factor,  it 
serves  to  provide  convergence  for  the  integral  in  eq.  2.17  as  t  approaches  infinity. 
If  we  define  the  operator  dk{t)  as  follows, 

dk^e-'^'d     and     e'^^^^dy^e-'"'^^  =  di,{t),  (2.20) 

we  can  express  eq.  2.17  as  the  transform  of  the  autocorrelation  function  of  the 
radiator  dipole,  as  shown  below: 

1  f°° 

n^)  =  ^  e'('^-''^)*Tr{dk^  •  dk(i)p}.  (2.21) 

Although  this  expression  is  interesting  from  a  physical  point  of  view,  we  need 
to  obtain  a  form  of  the  power  distribution  that  will  allow  us  to  obtain  ensemble 
averages  over  the  various  degrees  of  freedom  of  the  radiator-plasma  system.  Re- 
membering that  a  trace  is  invariant  under  cyclic  permutations,  we  arrive  at  the 
following  expression: 

1 

=   —       die^('^-^^)*Tr{dk-e-'^*/Vdk^e^^'/''}.  (2.22) 

Referring  to  eq.  2.22,  we  will  define  a  Liouville  operator  as  a  commutator  of  an 
arbitrary  operator  with  the  Hamiltonian,  such  that 


Lf  =  [H,  f]  and  e-^"/  =  e"'^*/'^  /e'^*/^  (2.23) 
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in  order  to  arrive  at  the  following  expression: 

1  r°° 

=  ^  /  ^  ^^Tr{dk  •  e-('^-^-")Vdkt}.  (2.24) 

Let  us  step  back  a  moment  and  consider  our  ultimate  goal.  We  currently  have 
an  expression  that  is  an  ensemble  average  over  the  entire  plasma  system.  We  wish 
to  obtain  ensemble  averages  over  the  various  degrees  of  freedom  of  the  radiator- 
plasma  system.  We  can  imagine  taking  the  transform  of  the  factor  e-'^'^"^-^^)'  to 
get  an  expression  of  the  form, 

=  (dk^i?(a;)dk),  (2.25) 

where 

R{u)-^  =  u-L-ie.  (2.26) 

We  would  prefer  to  take  the  average  over  individual  plasma  systems.  Therefore, 
we  will  now  write  the  full  Hamiltonian  of  the  radiator  and  the  plasma  electrons 
and  ions, 

H  =        +        +     +  H,  +  y,,  +  K,r  +  V,,r.  (2.27) 

Here,  the  first  term  refers  to  the  internal  states  of  the  radiator  while  the  second 
term  refers  to  the  center-of-mass  motion  of  the  radiator.  The  next  two  terms 
on  the  right  hand  side  of  the  equation  refer  to  the  kinetic  and  potential  energies 
of  the  respective  subsystems,  Hj  =  Tj  +  Vj.  The  interaction  potentials  between 
subsystems  j  and  k  are  given  by  Vj^k-  Each  term  of  the  hamiltonian  corresponds 
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to  a  Liouville  operator,  and  we  can  write  the  resolvent  R{u)  as  follows, 

R{uj)-'  =uj-  L;"'  +  L^""  +  L.  +  Le  +  Li,  +  Ler  +  ^>ir  "  i^.  (2.28) 

Since  the  radiator  dipole  is  not  dependent  on  the  coordinates  of  plasma  ions 
and  electrons,  we  can  take  an  ensemble  average  over  the  plasma  coordinates  of  the 
resolvent  R{u)). 

I{uj)  =  (dk^(i?(u;))p,a,dk)rad,  (2.29) 

where  the  subscripts  rad  and  plas  indicate  ensemble  averages  over  the  radiator 
and  plasma  degrees  of  freedom  respectively.  As  a  result  of  this  average,  we  will 
obtain  form  of  the  resolvent  of  the  form  below, 

{R{u))   =   l-uG{u;)-'G{uj)  (2.30) 
G{u)   =   ({a;-L-*  +  L7  +  L(r,5)  +  ^zv  +  (Le.)e  +  ^1e.-^c}-^)^(2.31) 

where  £  represents  an  electric  field  at  the  location  of  the  radiator  due  to  the 
plasma  ions  shielded  by  the  plasma  electrons  and  Qs  represents  an  average  over 
a  distribution  of  fields.  The  quantity  L{r,e)  represents  the  interaction  of  the 
radiator  with  this  ion  microfield.  The  factor  u,  along  with  the  form  of  R{uj)  in 
terms  of  G{uj),  represents  the  effects  due  to  the  fluctuations  of  the  ion  microfield. 
The  {Ler)e  represents  mean-field  eflFects  due  to  the  plasma  electrons,  the  term  Mer 
accounts  for  higher-order,  time-varying  efl!"ects  or  fluctuations  due  to  the  radiator 
interaction  with  the  plasma  electrons. 
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In  addition  to  averaging  over  the  plasma  degrees  of  freedom,  we  may  average 
over  the  center  of  mass  coordinates  of  the  radiator  in  order  to  obtain  an  expression 
for  Doppler  broadening.  We  will  address  this  issue  in  the  next  few  pages.  In 
the  following  section,  we  will  present  forms  of  the  operators  listed  in  eq.  2.31 
describing  the  effects  of  the  plasma  electrons  and  ions  using  certain  weak  coupling 
approximations. 

We  will  now  address  the  problem  of  radiator  motion.  We  will  treat  radia- 
tor velocity  as  being  statistically  independent  of  the  plasma  environment  and  the 
internal  radiator  state.  Another  way  to  express  this  is  to  factor  the  density  matrix: 

P       PcmPr,i,e^  (2.32) 

The  index,  r,  refers  to  the  internal  radiator  states,  cm  refers  to  the  the  center- 
of-mass  states,  while  i  and  e  refer  to  the  plasma  ions  and  electrons,  respectively. 
Although  we  have  shown  in  section  2.3  that  the  center-of-mass  motion  is  effectively 
independent  of  the  other  degrees  of  freedom  for  the  plasma  conditions  under  consid- 
eration, the  reader  can  find  quantum  theories  of  line  broadening  where  Stark  broad- 
ening is  assumed  to  be  statistically  dependent  on  radiator  motion  [44,  45,  46,  47]. 
For  cases  where  the  perturber  mass  is  much  less  than  the  radiator  mass  it  was 
shown  that  the  these  models  reduce  to  the  statistically  independent  case  [45,  46]. 

If  one  assumes  that  the  radiators  behave  like  an  ideal  gas,  their  velocity  distri- 
bution would  be  Maxwellian,  so  that 

''■^^{^y'^-""'''''  =  m-  (2.33) 
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The  mass  of  the  radiator  is  denoted  by  rrij.  and  its  velocity  by  v.  We  then  re- 
express  our  time  development  operators  in  order  to  separate  out  the  radiator's 
translational  motion  as  follows: 

where  H'  represents  the  combined  hamiltonian  for  the  radiator  internal  states  and 
the  plasma  system,  and  Hem  represents  the  hamiltonian  for  the  radiator  transla- 
tional states.  With  our  previous  assumption,  the  operators  on  the  radiator  center 
of  mass,  k  •  r,  now  commute  with  hamiltonian  for  the  internal  radiator  states  and 
the  plasma  system. 

Since  the  trace  in  eq.  2.22  is  invariant  under  cyclic  permutations  we  get  the 
following  expression: 


1  f 
=  ^^^^''^"""^7  (2-35) 

1    r°°  c 


1  /-ao 

-00 

1 

-00 

Tr{de-'''  '-(*)  •  e-'^''/Vr,i,ee'*'  '"de'^'*/'^}.  (2.36) 


If  the  velocity  of  the  center  of  mass  is  treated  as  constant  over  the  lifetime  of  the 
excited  state,  we  get  the  following  expression: 

1  r 
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We  can  take  advantage  of  the  fact  that  the  expression  inside  the  trace  is  the 
autocorrelation  function  of  the  radiator  dipole,  C{t).  Since  d  is  hermitian,  C{t) 
has  the  property, 

C{-t)  =  C{ty,  (2.38) 

we  employ  the  relation  above  to  write  I{uj)  as  follows, 

1      r°°      ,     ^  r 
I{u;)   =   -Re/    rfte'^"^-")'/ di;/(t;)e-''^-^'TV{d-e-'^''/V,i,ede'^''/'^K2.39) 

TT       •/  0  J 

By  evaluating  the  integral  over  v  and  applying  the  Fourier  convolution  theorem 
we  arrive  at  the  following  expression  [48]: 

where 

The  lineshape  for  the  radiator  in  the  center  of  mass  frame  is  given  by 

1  r°° 

Is{uj)  =  -Re      rfte'('^-«)*Tr{d  •  e-'^''/Vr,,ede'^'*/'^}.  (2.42) 

In  the  next  section,  we  will  the  address  the  problem  of  Stark  broadening  due  to 
ions  and  electrons.  For  the  time  being,  we  will  refer  to  Is{uj)  in  eq.  2.42  as  simply 
licv). 
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2.5    Stark  Broadening  Theory 


Our  theoretical  development  of /s(a;)  follows  standard  methods  regarding  Stark 
broadening  by  perturbing  electrons  and  ions.  The  current  calculation  incorporates 
a  model  of  the  plasma  that  accounts  for  the  motion  of  the  electrons  and  ions;  the 
latter  can  have  a  significant  broadening  effect  at  the  center  of  Ar  lines  when  the 
Ar  radiator  is  immersed  in  a  plasma  of  relatively  lighter  ions,  such  as  deuterium, 
as  is  the  case  in  ICF  implosions.  For  cases  where  the  ion  perturbers  are  more 
massive,  such  as  targets  filled  entirely  with  higher-Z  ions  such  as  Ne,  Ar,  and 
Kr,  ion  dynamic  effects  are  reduced  greatly  and  line-shapes  closely  resemble  those 
generated  with  the  static-ion  approximation  [18]. 

We  must  now  write  the  full  Hamiltonian  of  the  radiator  and  the  plasma  elec- 
trons and  ions,  excluding  the  translational  motion  of  the  radiator, 


Here,  the  first  three  terms  on  the  r.h.s.  refer  to  the  kinetic  and  potential  energies 
of  the  respective  subsystems,  Hj  =  Tj  +  Vj.  The  interaction  potentials  between 
subsystems  j  and  k  are  given  by  Vj^k-  Referring  to  eq.  2.42,  we  will  define  a  Liouville 
operator  as  a  commutator  of  an  arbitrary  operator  with  the  Hamiltonian,  such  that 


H'  =  Hr  +  H,  +  H,  +  Vi,e  +  Ve,r  +  V^,,. 


(2.43) 


L'f  =  [H',  f] 


and 


(2.44) 
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in  order  to  arrive  at  the  following  expression: 

I{uj)  =  -Re      (iie'('^-^^)'Tr{d  •  e-'^'*/Vr,i,ed}.  (2.45) 

Since  the  dipole  operator  acts  only  on  the  radiator  internal  states,  we  can  bring 
the  trace  over  plasma  states  within  the  brackets. 

1  /■°° 

-  -ReTr,       dt  e^(-«)M  •  Tri,e{e-^  '/Vr,i,ed}  (2.46) 

Traditionally,  the  effects  of  the  plasma  ions  are  addressed  at  this  point.  Since  we 
will  include  the  effects  of  ion  dynamics,  we  find  it  simpler  to  examine  the  radiator- 
plasma  electron  interactions  first.  Operators  acting  on  plasma  ion  coordinates  will 
be  indicated  with  the  subscript,  £i.  The  interaction  of  the  plasma  electrons  with 
the  plasma  ions  Vie  will  be  included  in  the  plasma  ion  hamiltonian  to  account  for 
screening  of  the  plasma  ions, 

H'  =  Hr  +  V^^i^  +Hi  +  He  +  Ve,r,  (2.47) 

where 

KA=Vi,r  +  Vi^e.  (2.48) 

In  order  to  address  the  eflfects  of  the  plasma  electrons,  we  will  be  presenting  an 
electron  shift  and  width  operator  that  was  obtained  employing  standard  kinetic 
theory  techniques  as  described  by  Hussey  [38]  and  simplified  using  certain  weak 
coupling  limits  by  Dufty  et  al.  [49].  Due  to  the  sheer  length  and  complexity  of 
the  formalism,  we  will  present  the  results  of  Dufty  et  al.  and  proceed  from  that 
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point.  First,  we  must  describe  procedures  to  separate  out  the  effects  of  the  plasma 
electron  interactions  in  the  hamiltonian  and  the  corresponding  Liouville  operator. 
The  Liouville  operator  corresponding  to  a  single  plasma  electron,  can  be 

written  as  follows, 

L(l)/=[r(l)  +  Vo(l),/],  (2.49) 

where  Vo{l)  represents  the  monopole  interaction  of  the  electron  with  the  radiating 
ion,  which  is  dependent  on  only  the  perturber  coordinates  if  we  treat  the  center 
of  mass  of  the  radiator  as  an  inertial  frame.  The  interaction  potential  between  an 
individual  plasma  electron  and  the  internal  states  of  the  radiator,  Vi{r,  1),  can  also 
be  expressed  in  terms  of  a  Liouville  operator.  In  order  to  account  for  the  screening 
due  to  plasma  electrons  we  will  also  introduce  the  operator  C{r,  1), 

[Li{r,l)f  =  [V,ir,l)J]     and     [C{r,l)f  =  [V,{r,l),f],  (2.50) 

where  Vs{r,  1)  represents  the  screened  interaction  potential  between  the  radiator 
and  one  perturbing  electron. 

Using  the  preceding  definitions,  the  line  profile  can  be  expressed  as  follows: 

=  -^ImTr.d  •  Tr,,ei?(a;)/(r,  5)d},  (2.51) 

where 

R{u)-'  =u-  L{r)  -  Li{r,i)  -B-  M{uj),  (2.52) 

and  fir,£)  is  a  distribution  function  for  the  internal  radiator  states  and  plasma  ion 
states.  The  operators  B  and  M{uj)  make  up  our  electron  shift  and  width  operator, 
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which  correspond  to  a  mean-field,  (Ler)e,  and  a  fluctuation  term  Mer  respectively. 
Dufty  et  al.  obtained  forms  for  B  and  M{uj)  similar  to  the  ones  below  [49], 


M{cj)  =  neTri£(r,  l)f{rl){u  -  L{r)  -  L,{r,i)  -  L(l)  -  C{r,  l)}-'C{r,  1).  (2.54) 

In  Chapter  III,  we  will  discuss  the  operators  B  and  M{u;)  in  detail.  We  will 
now  address  the  issue  of  the  plasma  ion  interactions.  The  assumption  that  the 
ions  are  stationary  during  the  time  of  interest  will  not  be  made.  We  will  follow 
the  formalism  introduced  by  Boercker  et  al,  in  ref.  [40]  and  outlined  by  Griem 
in  ref.  [1].  The  procedure  to  separate  out  the  effects  of  the  plasma  ions  is  similar 
to  the  one  employed  to  obtain  B  and  M(a;);  however,  there  are  some  important 
distinctions.  The  plasma  ions  are  treated  collectively  and  their  interactions  with 
the  radiator  are  represented  by  an  electric  field  acting  with  the  radiator  dipole 


We  begin  by  defining  the  dipole  moment  operator,  ¥{r,S;u),  where  the  eff"ects 
the  radiator-plasma  electron  interactions  have  already  been  accounted  for, 


5  =  neTriLi(r,l)/(rl)/(r) 


-1 


(2.53) 


and 


moment. 


F(r,5;t) 


(2.55) 


where 


L'  =  L{r)^L^{r,e)+B^M. 


(2.56) 
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We  employ  a  dummy  variable  E  to  represent  the  microfield  due  to  one  of  the  plasma 
ion  configurations  so  that  the  trace  over  ion  coordinates  can  be  expressed  as  an 
integral  over  a  microfield  distribution  function,  Q{S),  which  is  shown  below, 

Tr^F(r,  t^t)   =   J  dSTv^  [F(r,  £;  t)S{£  -  £)] 

=   I  d£Q{£)F{r,£;t),  (2.57) 

where 

F(r,  £■  t)  =  Q{£)-'Tx^  [F(r,  5;  t)S{£  -  £)]  .  (2.58) 

Using  the  time  dependent  Shrodinger  equation,  we  obtain  the  following  dynam- 
ical equation: 


^F(r,  £;t)  +  i  [L{r)  +  L^{r,  £)  +  B  +  M]  F(r,  £;  t)  =  ^F(r,  £■  t) 


,  (2.59) 


where  Li{r,£)  is  the  Liouville  operator  associated  with  the  radiator  interaction 
with  a  given  microfield  £.  The  term  on  the  right  hand  side  is  analogous  to  the 
collision  term  in  Boltzmann  equation,  and  hence  can  be  rewritten  as 


=  I  d£'  [W{£'  ^  £)F{r,  £';  t)  -  W{£  ^  £')F{r,  £■  t)] ,  (2.60) 


where  W{£'  ->  £)  represents  a  transition  rate  from  microfield  £'  to  £.  If  we  take 
advantage  of  the  symmetry  relations  between  W{£'  £)  and  W{£  -»  £')  and  use 
the  simplest  choice  for  Wwe  get  the  following  expression: 


W{£'  ^£)  =  uQ{£'), 


(2.61) 
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so  that  we  may  re-express  eq.  2.59  as  follows: 


dt 


+  iL{r)  +  ill  (r,  S)  +  iB  +  iM 


=  —u 


'{r,£;t)-  I  d£'Q{£')F{r,£'-t) 


(2.62) 


where  the  quantity  nu  is  related  to  the  momentum  autocorrelation  function.  We 
now  take  the  Laplace  transform  of  F{r,S;  t)  to  arrive  at  the  following  expression: 

F(r,<f;u;)  =  [d  +  ul{uj)][-iu  +    +  iL{r)  +  iLi{r,£)  +  iB  +  iM{uj)]-\  (2.63) 


where 


[(a;)  =  j  d£'Q{£')f{r,£'-uj). 


(2.64) 


By  integrating  Y{r,£\uj)  over  the  microfield  distribution  and  some  manipula- 
tion ,  we  can  obtain  an  expression  for  F(r,  £;  w), 

F(r,£;  Aa;)  =  i[Au;  +  iv  -  L,{r,£)  -B-  M{^u,£)]-\\  -  uG{^uj)]-^d,  (2.65) 


where 


and 


G(Aa;)  =  j  d£Q{£)n{Au), 


(2.66) 


IZiAu)  =  (Au  +  w-Li  (r,  £)-B-  M{Au),  £))-\ 


(2.67) 
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The  quantity  Ao;  is  defined  as  u  -  L{r).  The  expression  for  F(r,  u)  can  now 
be  substituted  into  eq.  2.51  to  arrive  at  the  following  expression  for  /(a;): 

/(a;)  =  -hmTvrd  ■  [1  -  i^G{Au;)]-'G{Auj)f{r,  £)d.  (2.68) 

At  this  point,  the  overall  lineshape  expression  is  complete,  although  few  as- 
sumptions have  been  made  about  the  form  or  details  of  the  radiator,  plasma  ion, 
or  plasma  electron  Hamiltonian.  In  the  next  chapter,  we  will  explore  in  greater 
detail  the  electron  shift  and  width  operators,  B  and  M{uj).  We  will  obtain  useful 
expressions  for  the  operators  employing  a  full  Coulomb  form  of  the  radiator-plasma 
electron  interaction.  We  will  also  address  the  problem  of  calculating  B  and  M{u) 
for  multi-electron  radiators. 


CHAPTER  3 

THE  ELECTRON  SHIFT  AND  WIDTH  OPERATOR 

In  this  chapter,  we  examine  the  electron  shift  and  width  operator  in  more  detail. 
A  form  the  radiator-perturbing  electron  interaction  potential  is  introduced  and 
applied  to  the  expressions  for  B  and  M{uj)  introduced  in  Chapter  II.  We  will  apply 
an  approximation,  second  order  in  the  radiator-perturbing  electron  interaction, 
to  treat  highly  ionized,  moderate-Z  radiators.  Its  validity  will  be  discussed  and 
justified.  We  then  outline  the  formalism  with  which  we  address  the  problem  of 
line  radiation  from  multi-electron  radiators.  We  present  the  calculational  forms  of 
the  first  and  second  order  terms  of  the  shift  and  width  operator  for  very  simple 
multi-electron  cases  to  show  the  expressions  reduce  to  a  form  very  similar  to  the 
hydrogenic  case  presented  by  Woltz  [13,  14].  Specifically,  we  examine  radiators 
with  no  more  than  two  electrons  in  an  upper  manifold. 

3.1    Radiator-Perturbing  Electron  Interaction 

In  order  to  comment  further  on  the  operators  B  and  M{Auj),  we  first  write 
the  Hamiltonian,  H',  in  which  the  effects  of  the  radiator-perturbing  ion  interaction 
have  already  been  incorporated  into  the  ion  microfield: 

H'  =  Hr+e  £Zr  +  K,  +  Kr  +  Ke.  (3.1) 
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The  first  term  on  the  right  hand  side  is  the  Hamiltonian  for  an  isolated  radi- 
ator. In  the  second  term,  the  effects  of  the  perturbing  ions  have  been  expressed 
in  terms  of  the  radiator  dipole  interacting  with  the  ion  microfield.  It  may  seem 
counter-intuitive  to  suggest  that  the  dipole  approximation  is  valid  for  radiator  in- 
teractions with  plasma  ions,  but  invalid  for  interactions  with  plasma  electrons  at 
high  densities.  However,  for  highly  charged  radiators,  the  perturbing  ions  experi- 
ence a  strong  repulsion  due  to  the  monopole  interaction  with  the  radiator,  which 
is  included  in  the  microfield  calculation.  Therefore,  the  ions  remain  suflficiently 
distant  from  the  radiator  for  the  dipole  approximation  to  be  valid.  Higher  order 
corrections  to  the  dipole  approximation  have  been  evaluated  and  appear  to  have  a 
relatively  small  effect  [50].  This  is  not  the  case  for  plasma  electrons  which  experi- 
ence an  equally  strong  attraction  and  may  have  the  opportunity  to  penetrate  the 
radiator  orbitals.  Continuing  with  the  rest  of  the  Hamiltonian,  Ke  is  the  kinetic 
energy  of  the  plasma  electrons,  Ver  is  the  radiator-perturbing  electron  interaction, 
and  Vee  represents  the  perturbing  electron-electron  interactions. 

In  order  to  calculate  B  and  M(Aa;),  the  interaction  potential  between  the 
radiator  and  a  given  plasma  electron,  is  divided  into  a  long-range  monopole  term, 
and  a  short-range  interaction,  F(^'(r,  j),  as  follows,  where  V^^\j)  is  the 
monopole  interaction  of  the  jth  perturbing  electron  with  the  radiator,  and  Xj  is 
the  coordinate  of  that  electron, 

y(o)(,)  =  Jl^^  (3.2) 

Xj 
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and  V^^^r,j)  is  the  sum  of  the  interactions  between  Nr  individual  radiator  electron 
and  the  perturbing  electron, 

v""(..)  =  - 1)  .  - 1) 

^  -Ar^r-]-  (3.3) 

\\Xj      Xrj^^l       Xj  J 

The  coordinates  of  the  iih  radiator  electron  are  denoted  by  .  Later,  we  will  take 
advantage  of  the  fact  that  V^^^{r,j)  has  been  expressed  as  a  sum  of  symmetric 
operators  on  the  radiator  electron  coordinates.  is  independent  of  radia- 

tor electron  coordinates  and  may  be  included  in  the  Hamiltonian  of  the  plasma 
electrons, 

mj)  =  ^-V^'\j).  (3.4) 

This  form  lends  itself  to  employing  Coulomb  wavefunctions  to  describe  the  motion 
of  the  perturbing  electrons,  as  was  done  in  the  past  [12,  13,  14]. 

We  can  now  define  Liouville  operators,  L(r,  1)  and  Li(r,  1),  as  follows, 

L{r,  1)/  =  [H'{r,  1),  /]  L,{r,  l)f  =  [V^'\r,  1),  /],  (3.5) 

where  /  is  an  arbitrary  operator.  H'{r,  1)  represents  the  Hamiltonian  for  the 
radiator  interacting  with  the  ion  microfield  and  a  single  plasma  electron. 


H'{r,  1)  =  Hr.  +  eSZr  +  H{1)  +  V^'\r,  1) 


(3.6) 
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As  was  shown  in  the  previous  chapter,  the  electron  interaction  operators  B  and 
M{Au})  can  be  expressed  in  terms  of  L{r,  1)  and  Li{r,  1)  [12,  38,  49]. 


B   =  nTriLi(r,l)/(rl)/(r)-i 
M{lo)   =  nTri[£(r,l)/(rl)[a;-L(r)-Li(r,£)-L(l)-£(r,l)] 


-1 


(3.7) 


^r,  l)/(r) 


-1 


(3.8) 


where  /(r)  and  /(rl)  are  distribution  functions  for  the  isolated  radiator  and  the 
radiator-perturbing  electron  system,  respectively.  These  expressions  were  derived 
from  an  exact  expression  obtained  by  Hussey  [38,  39]  and  simplified  by  Dufty  et 
al.  [49].  In  eq.  3.8  the  effect  of  plasma  electron  correlations  are  incorporated  by 
Dufty  and  Boercker  [49],  into  £(r,  1),  which  is  simply  statically  shielded  forms 
of  Li(r,  1).  The  quantity,  B,  represents  the  time-independent  portion  of  the  in- 
teraction of  the  plasma  electrons  with  the  radiator  corresponding  to  the  initial 
correlations.  It  is  real  and  provides  the  bulk  of  the  shift  for  the  cases  under  con- 
sideration. M{u})  is  the  transform  of  the  autocorrelation  function  of  the  electron 
interaction  potential  Vi{r,  1)  .  At  this  point,  B  and  M(uj)  are  exact  to  all  orders 
in  Vi(r,l). 


We  will  calculate  B  and  M{u)  to  second  order  in  the  radiator-perturbing  elec- 
tron interaction,  Li{r,  1),  using  a  full-Coulomb  expression  of  the  interaction  poten- 
tial [12,  51].  There  has  been  some  confusion  concerning  the  relationship  between 
the  relative  size  of  the  first  order  shift  and  the  shift  contribution  found  by  expand- 
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ing  the  electron  shift  and  width  term  in  powers  of  the  radiator-plasma  electron 
interaction.  The  size  of  the  first  order  term  of  the  shift  has  little  to  do  with  the 
appropriateness  of  terminating  electron  shift-  and  width  operator  expansion  with 
the  second  order  term  for  a  given  temperature  and  density.  As  will  be  shown  later, 
the  first  order  term  has  no  contribution  from  the  multipole  expansion  beyond  the 
penetrating  monopole  term,  therefore  it  is  inappropriate  to  compare  its  contribu- 
tion with  that  from  terms  containing  contributions  from  the  dipole  interaction  and 
beyond,  in  order  to  determine  the  appropriateness  of  the  second  order  treatment. 
In  fact,  one  could  anticipate  a  case  where  the  first  order  shift  is  on  the  order  of  the 
second  order  term  or  smaller  and  yet  the  second  order  approximation  is  valid. 

Let  us  now  specify  that  B^^^  represents  the  effects  of  the  shielding  due  to  plasma 
electrons  penetrating  the  radiator  orbital  on  the  energy  of  the  bound  electron  with 
no  higher  order  corrections,  then  the  expression  B  -  B^^^  +  M{Auj)  represents  the 
corrections  to  that  term.  We  have  rephrased  the  question  to  ask  whether  truncat- 
ing the  electron  shift  and  width  operator  gives  not  only  a  good  result  for  the  width 
but  an  adequate  correction  to  the  shift.  In  order  to  apply  the  second  order  approx- 
imation, two  conditions  must  be  satisfied.  First,  large  portions  of  the  broadened 
fines  must  satisfy  the  inequality,  Au  ^  2LUpe,  where  Aa;  is  the  frequency  separation 
from  line  center  and  ujpe  is  the  plasma  frequency.  This  is  relevant  because  second- 
order  calculations  have  been  shown  to  be  valid  in  this  range  [51,  52].  Second,  the 
plasma  electrons  must  undergo  predominantly  weak  collisions,  which  are  defined 
as  collisions  where  the  change  in  the  electron's  momentum  is  much  smaller  than 
the  total  momentum.  Another  way  to  express  this  idea  is  to  state  the  impulse 
experienced  by  the  perturber  is  small.   The  plasmas  under  consideration  must 
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be  characterized  by  sufficiently  high  temperatures  so  that  even  collisions  where 
the  perturbing  electron  penetrates  the  bound  electron  orbital  can  be  described  as 
'weak.'  In  other  words,  the  average  plasma  electron  has  a  large  kinetic  energy 
and  spends  a  relatively  short  time  near  the  radiator  itself  and  thus  receives  only  a 
small  impulse,  relative  to  its  momentum.  The  appropriate  temperature  regime  is 
defined  by, 

A;r>^t?^(Ryd)      [30].  (3.9) 

This  expression  was  obtained  by  comparing  the  difference  in  the  kinetic  energy  of 
an  electron  far  from  the  radiator  to  its  value  at  the  point  of  closest  approach. 

There  are  two  reasons  to  believe  that  the  second-order  approximation  is  an 
appropriate  choice.  First,  large  portions  of  the  a-,  /?-,  and  7-lines  of  K-shell  Ar, 
for  the  plasma  conditions  under  consideration,  satisfy  the  condition  that  the  width 
of  line  be  no  more  than  two  times  the  plasma  frequency  [52].  Second,  the  and 
7-lines  satisfy  the  inequality  in  eq.  3.9  for  the  temperature  conditions  observed  in 
recent  experiments  [20,  22,  21]. 

It  will  be  shown  later  in  this  section  that  B,  the  static  shift,  can  be  expanded 
into  a  term  first  order  in  the  interaction  potential  and  a  term  containing  all  higher 
order  eflfects  [31].  We  will  outline  the  process  here.  We  can  expand  the  distribution 
function,  /(r,  1),  and  thus  B,  in  terms  of  V^^\t,  1).  By  successively  closing  kinetic 
equations  in  the  equilibrium  Bogliubov-Born-Green-Kirkwood-Yvon  (BBGKY)  hi- 
erarchy, Boercker  and  Dufty  obtained  the  following  expression  for  /(r,  1): 


/(r,  1)  =  Zexp[-m{r)  +  H{1)  +  Vs{r,  1)]}      [49],  (3.10) 
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where  Vs{r,  1)  is  a  statically  screened  radiator-electron  interaction  potential.  Boer- 
cker  and  Iglesias  expanded  /(r,  1)  in  orders  of  Vs{r,  1)  and  obtained  expressions  for 
B^^^  and  B^'^'>  [31,  53].  We  begin  by  making  the  following  definitions: 

H'o  =  H{r)  +  H{l),  V,  =  V,{r,l).  (3.11) 

We  express  the  distribution  function  in  terms  of  H'q  and  a  "development  operator" 
W,0), 

e-""'  =  e-^^'«iY(/3, 0)  and  V,{t)  =  e'"'oV,e-'"'o,  (3.12) 

where  U{(3,0)  is  defined  as  follows: 

mO)  =  Ei-^r-^  I  dn---      rfr„T,[V^i(ri)---Fi(r„)].  (3.13) 
U  will  be  approximated  to  first  order  in  Vi : 

U{P,  ^)  =  '^-  [  dre-'^oVie-'^o  +  0{V,)\  (3.14) 

This  approximation  is  possible  only  if  the  integral  over  r  is  much  smaller  than  one. 
Since  r  is  no  larger  than  kT-\  the  average  interaction  potential  must  be  no  larger 
than  kT. 

If  we  define  /o(r)  and  /o(l)  using  the  following  equations, 


/o(r)  =  Z,e-^^W 


and 


/o(l)  =  Zee-^«(^ 


(3.15) 
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where  Zr  and  Zg  are  the  respective  partition  functions  for  /o(r)  and  /o(l),  we  can 
re-express  /(r,  1)  in  terms  of  /o(r)  and  /(I)  and  0). 

/(r,l)-/o(0/o(l)W,0).  (3.16) 

If  we  substitute  this  expression  into  eq.  3.7  using  the  expression  for  W(/3,0)  from 
eq.  3.14,  we  arrive  at  the  following  expressions  for  B'^^^  and  B^^^: 

=   nTri[Li(r,l)/o(l)],  (3.17) 
5(2)   =  -nTri{Li(r,l)/o(r)/o(l) 

X  f  dre^[^W+^(i)]K(r,  i)e-^[^('-)+^(i)] /o(r)-^}.  (3.18) 

Starting  with  eqs.  3.17  and  3.18,  we  can  express  the  trace  over  perturber  states 
in  terms  of  integrals  over  perturber  momenta.  We  introduce  complete  sets  of 
radiator  eigenstates  to  obtain  a  matrix  representation,  and  evaluate  the  integral 
over  T.  All  of  this  is  described  in  detail  in  Appendix  B.  We  rewrite  eqs.  B.8 
and  B.13  below, 

B(^)   =   nA^  I rfke-^  [(/i;  k|Vi(r,  -  B^^},6^^\  (3.19) 

=   -n/rfki|  c/k2E(//ki|yi(r,l)|//"k2)(/k2|K(r,l)|Ai)/o(/i") 


d(2) 


xfo{k2)fo{fi')--  ;  -1^2  +  (3.20) 


where  //  and  v  represent  arbitrary  radiator  states,  and  k  represents  the  wavenumber 
of  the  perturber  state.  bI]^,  has  the  same  form  as  the  first  term  but  the  matrix 
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element  is  between  states  i/  and  u'.  Also,  O.T  refers  to  a  second  term  in  B^"^^  where 
the  one  of  the  matrix  elements  is  between  states  u  and  v'.  For  the  sake  of  brevity, 
we  will  examine  only  this  first  term  for  now.  We  next  make  an  approximation  in 
order  to  express  the  matrix  elements  in  terms  of  initial  and  final  radiator  states. 

Since  B  is  a  tetradic  operator  acting  on  the  radiator  dipole,  d,  it  can  be  written 
in  the  following  form: 

[Bd]f,^  =      ^M-'^^v'd^v  (3.21) 

We  will  assume  that  initial  and  final  states  of  the  radiative  transition  are  sufficiently 
separated  in  energy  so  that  the  matrix  elements  of  V^i  (r,  1)  between  initial  and  final 
states  are  negligible  compared  to  matrix  elements  among  initial  states  or  among 
final  states.  This  is  the  "No  Quenching"  approximation;  we  assume  there  are 
no  non-radiative  transitions  between  initial  and  final  states.  Therefore,  since  the 
radiator  dipole  must  act  between  an  initial  and  final  state,  we  will  substitute  \i} 
for  I//)  and  |/)  for  The  matrix  elements  of  the  radiator-perturbing  electron 
interaction  Vi{r,  1)  between  initial  and  final  states  will  be  assumed  to  equal  zero. 
Below,  we  have  re-expressed  eqs.  B.8  and  B.13  using  a  partial  wave  expansion  for 
the  perturber  states, 

1=0  m=-l 

-(/;  k,  I,  m\V^'\r,         k,  I,  m)4'}  (3.22) 


and 


42 


^ifi'f'   =   -Sff'nXT  J  dki  J  dk2j2Di,i",i'{ki,k2) 

—  hui') 

(3.23) 


exp  I -P^  -  exp  +  hu,.  -  huj,) 


where 


oo  /l  CX)  <2 

Di,v'Ah,h)   =   klklY,    J2    E    11  {i;kulum,\V,{r,l)\i"-k2j2,m2) 

/i=0  mi  =— ii  /2=0  m2  =— (2 

{i";  k2,  h,  Tn2\Vs{r,        ki,h,mi).  (3.24) 


All  the  information  about  the  matrix  elements  of  the  interaction  potential  for  the 
second  order  term  is  contained  in  A,i",t'(^i,  ^2),  which  we  will  put  in  calculable 
form  in  Section3.5.  The  product  'kef  of  the  initial  radiator  state,  i,  and  the 
perturbing  electron  state  are  represented  by  \i;k,l,m)  The  plasma  electron  state 
is  represented  by  a  partial  wave  expansion,  where  /  and  m  describe  the  angular 
momentum  state,  and  k  is  the  wave  number  defined  as  follows: 

H{l)\k,  I,  m)  =  -—\k,  I,  m).  (3.25) 

Returning  to  eq.  3.8,  we  see  that  M{uj)  is  already  at  least  second  order  in 
£(r,  1),  the  screened  form  of  Li(r,  1).  For  a  strictly  second  order  form  of  M(a;), 
the  expression  L{r)  +  L(l)  +  £(r,  1)  in  the  resolvent  should  be  substituted  with 
L(r)+L(l).  However,  we  will  instead  employ  L(r)+L(l)+S(^),  since  the  first  order 
shift  provides  a  non-perturbative  modification  of  the  energy  levels,  with  respect  to 
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the  plasma  frequency;  therefore,  this  is  not  strictly  a  second  order  expression.  In 
the  form  presented  below,  the  term,  M{Au)^',  which  describes  the  interaction  of 
the  initial  states,  is  re-expressed  as  a  integral  over  time.  There  is  a  corresponding 
term  for  the  interaction  of  the  final  states,  as  well  as  interference  terms  that  are 
written  in  eq.  B.23, 

M{u)ii.   =   -inX^  j  dkie-'^'-2^  I  dk2 1  dtJ2Di,i",4ki,k2) 

.26) 


xexp{i{hAuji»f  -  B-"/  +        -        +  ie)t}Sff>.  (3. 


While  the  quantities,  B^^^  and  S^^)  are  real  and  contribute  only  to  the  shift, 
the  dynamic  term,  M{Auj),  is  complex.  It  contributes  to  both  the  shift  and  width. 
We  now  evaluate  the  integral  over  time,  remembering  that  we  are  taking  the  limit 

as  e  — )^  0, 


fOO  1 

-I  /    die^[/(-)+«l*   =    lim  —4  

f{uj)  -  ie 


fOO  1 

Jo  f{uj) 


(3.27) 


where  P  stands  for  the  Cauchy  principle  part.  It  is  relatively  straightforward  to 
evaluate  the  imaginary  part  of  M(Aa;).  Following  the  methods  of  O'Brien  and 
Hooper  [10],  we  will  use  the  following  relation,  where  F  is  a  function  of  k,  to 
evaluate  the  integral  of  the  real  part  of  M{I\oj): 

P/o         J:ki-k^  ->0-  (3-28) 
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Using  the  expression  above,  we  get  the  following  expressions  for  Aa;j" /  >  0,  with 
similar  expressions  for  Awi" /  <  0: 


Acjj"/  >  0, 


(3.29) 


ReM{u)u'  =  ^      [dki  f  dk^ihAco^.j  +  ^{k^  -  k^))-' 
IT''    ~z  J       J  2m 


e  2m 


9.m.\A.ii)-    I  \ 

Di,i",i'  {ki ,  k2)  -  Di^.^  \  \kl+       ^^"^^  j     ,  k2 


Aui>'f  <  0, 


(3.30) 


ImM(a;)i 


+hAw. 


[kj+'-^^y^" 

Aa;i"/>0,  (3.31) 


ImM(a;)ii/  =  - 


AtUj// /  <  0. 


(3.32) 


We  now  have  expressions  for  the  width  and  shift  operator  to  second  order  in 
the  radiator-perturber  interaction  potential  in  terms  of  the  matrix  elements  of  that 
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interaction.  The  next  step  involves  obtaining  calculable  forms  of  the  single  matrix 
element  in  B^^)  and  the  expression  ^2)  for  B^"^^  and  M(Aa;) 

Boercker  and  Iglesias  showed  for  neutral  radiators  that  there  is  significant  can- 
celation in  the  portion  of  the  shift  term  which  is  second  order  in  the  interaction 
potential,  between  B^^)  ^^^^  ReM(a;).  They  find  the  real  portion  of  M(0)  equals 
This  is  important  for  the  neutral  case,  where  the  first  order  term  is  be- 
lieved to  vanish  [31].  It  can  be  shown  that  the  second-order  behavior  is  the  same 
for  charged  radiators;  however,  the  resulting  cancelation  is  a  relatively  small  effect 
for  highly-charged  radiators,  since  the  first  order  shift  is  dominant. 

So  far,  our  expressions  have  not  specified  the  details  of  the  interaction.  In 
the  next  section,  we  will  address  the  challenges  of  calculating  the  matrix  element, 
(z,  A;i,/i,mi|V(i)(r,  l)|z';  A;2,/2,"^2),  for  a  multi-electron  radiator. 


3.3    Multi-Electron  Formalism 

In  equations  3.23  and  3.26,  the  ket  \i;k,l,m)  represents  a  direct  product  of 
the  energy  eigenstates  for  the  radiator  and  perturbing  electron.  However,  since 
the  interaction  potential  Vi(r,  1)  does  not  act  on  the  spin,  the  matrix  elements 
of  the  interactions  are  easiest  to  calculate  in  the  CS  representation  using  product 
wavefunctions  for  the  radiator  electrons.  Some  multi-electron  atom-  and  ion  config- 
urations can  experience  intermediate  coupling  between  angular  momentum  states. 
Hence,  in  order  to  calculate  the  matrix  elements  of  the  interactions,  it  is  necessary 
to  first  transform  from  the  energy  eigenstates  to  specific  electron  configurations  in 
the  CS  representation.  A  transformation  matrix  can  be  obtained  by  diagonalizing 
the  matrix  of  the  isolated  radiator  Hamiltonian  but  retaining  the  eigenvectors  for 
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the  particular  representation  most  useful  for  the  model.  Since  Cowan  does  not 
provide  a  specific  notation  for  this  transformation  matrix,  we  shall  refer  to  it  as 
Ta,i.  The  subscript,  a,  refers  to  all  the  relevant  quantum  numbers,  except  J  and 
M  which  represent  the  total  angular  momentum  and  its  azimuthal  component, 
of  a  given  configuration  in  the  appropriate  representation.  We  shall  employ  the 
geneological  CS  basis  to  express  the  states  of  the  angular  momenta.  This  repre- 
sentation is  described  in  Appendix  C  and  in  Section  4.10  of  Cowan's  book,  The 
Theory  of  Atomic  Structure  and  Spectra,  1981  [54].  We  note  that  for  an  isolated 
radiator  the  potential  is  spherically  symmetric,  and  the  total  angular  momentum 
is  conserved.  Therefore,  Ta,i  is  block  diagonal  in  J  and  invariant  in  A^.  We  will 
designate  individual  blocks  by  Ta,i{J)  and  the  relevant  basis  states  by  \<ifa)jM- 
Therefore,  the  ith  wavefunction,  \i)jM,  in  the  energy  basis  can  be  expressed  as 
follows: 

\i)jM  =  E  \'^a)jMTaAJ)-  (3.33) 
a 

This  results  in  the  following  matrix  elements  of  the  radiator-perturbing  electron 
interaction: 

jM{i;  ki,li,mi\Vs{r,        k2,  k,  rn2)j'M' 
=   E  Tl,a(^)^^(*a;  A;i,  /i,  mi|K(r,  l)|*a';  A;2,  m2)j'M'T^a',AJ'){^-^^) 

a, a' 

In  order  to  calculate  the  matrix  elements  of  the  interaction,  it  is  important 
to  remember  that  V^^\r,  1)  was  defined  in  eq.  3.3  as  a  sum  of  symmetric,  one- 
electron  operators.  This  allows  us  to  employ  tools  described  in  detail  in  Cowan  [54]. 
The  first  step  is  to  express  our  interaction  in  terms  of  reduced  matrix  elements. 
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F^^^(r,  1)  can  be  expressed  as  a  multipole  expansion  in  terms  of  reduced  tensor 
operators  for  each  radiator  electron  coordinate,  Xrj , 

00  t 

V^'Hr,,  1)  =  E  E  A\xr^,x,)C^'\x,^)Cl'>{x,),  (3.35) 

t-O  g=-t 

where  A''  and  C^'^  describe  the  radial  and  angular  portions  of  the  multipole  ex- 
pansion. 

where  xi  is  the  plasma  electron  coordinate.  x<  and  x^  are  the  lesser  or  greater, 
respectively,  of  the  perturbing  electron  coordinates  and  relevant  radiator  electron 
coordinates.  In  order  to  account  for  screening  due  to  electron  correlations,  we 
will  follow  the  usual  practice  of  introducing  a  cutoff  at  the  Debye  length  when  we 
integrate  over  the  perturbing  electron  coordinates.  One  may  also  use  alternative 
methods,  such  as  introducing  an  exponential  dependent  on  the  Debye  length  to 
the  potential,  but  the  choice  of  a  cutoff  allows  us  to  compare  our  results  more 
consistently  with  other  calculations  that  also  employed  a  cutoff.  Please  note  that 
the  interaction  in  B^^^  is  not  screened  and  does  not  need  a  cutoff.  In  our  next  step, 
we  employ  the  Wigner-Eckart  theorem  to  obtain  reduced  matrix  elements  [55, 
56].  We  will  express  our  potential  as  a  product  of  operators  dependent  on  either 
perturbing  electron  or  radiator  coordinates  so  that  we  may  evaluate  the  matrix 
elements  independently,  this  may  not  always  be  possible  for  other  potentials,  but 
it  is  necessary  to  be  able  to  evaluate  the  angular  portions  of  matrix  elements 
separately.  We  have  already  expressed  our  perturber  wavefunctions  in  terms  of  a 
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partial  wave  expansion  and  now  apply  the  Wigner-Eckart  theorem  [55,  56], 


jx(*a;  fci,  /i,  mi\Vs{r,  k2,  k,  rn2)j'M' 

h,l2  t,q  j=l 


^_-^^J-M+li-mi+q 


J       t  J 

\^-M    q   M'  j 


^  -mi    -q   7712  j 


where  we  have  the  used  the  relation. 


(3.37) 


(3.38) 


The  quantity  represents  the  portion  of  ^  dependent  only  on  the  coordinates  of 
the  jth  radiator  electron  and  A'\  is  the  portion  of  A*  ^  dependent  on  the  perturber 
coordinates.  The  manner  in  which  the  first  reduced  matrix  element  is  evaluated 
depends  on  whether  states  a  and  a'  are  in  the  same  electron  configuration.  The 
configuration  of  a  given  states  is  designated  only  by  the  number  of  electrons  in 
each  subshell  Wg,  in  which  the  subshells  are  in  turn  are  denoted  by  their  principle 
quantum  number  n^and  orbital  angular  momentum  for  example  As2pls  or  3dls^. 
If  two  states  are  in  the  same  configuration,  the  matrix  elements  between  them  are 
not  strictly  diagonal,  since  the  sums  of  the  angular  momenta  will  vary.  In  the  cases 
where  we  evaluate  a  matrix  element  between  states  with  the  same  configuration,  we 
will  obtain  contributions  to  our  calculational  result  from  every  occupied  subshell. 
We  refer  to  eq.  9.80  from  Cowan  [54],  since  our  interaction  potential  is  represented 


49 

as  a  sum  over  symmetric,  one-electron  operators, 

(*a,j|  E^il^a,j)  =  (3-39) 

j=l  s=l 

where  s  refers  to  the  sum  over  relevant  bound  subshells,  Wg  represents  the  occupa- 
tion number  of  the  subshell,  q  is  the  total  number  of  relevant  subshells,  acts  on 
the  last  electron  coordinate  in  subshell  s,  and  '^a,j  is  a  coupled  product  function 
with  basic  unpermuted  coordinate  ordering.  In  order  to  evaluate  matrix  elements 
over  the  individual  subshells  we  must  first  uncouple  the  total  spin  <S,  from  the 
orbital  angular  momentum  Cq.  In  order  to  understand  how  the  angular  momen- 
tum of  the  individual  subshell  Lg  is  uncoupled  from  the  total  we  must  know  how 
the  momenta  are  originally  summed.  We  employ  the  geneological  LS  representa- 
tion, where  the  orbital  angular  momenta  of  the  electrons  in  a  given  subshell  are 
summed  to  obtain  the  total  for  that  subshell  Lg.  Then  the  total  momentum  for 
that  subshell  is  added  to  the  intermediate  total  orbital  angular  momentum  includ- 
ing the  previous  subshell  £^-1  to  obtain  the  new  intermediate  total,  In  order 
to  evaluate  the  matrix  element  for  a  specific  subshell,  we  reverse  the  process  by 
successively  uncoupling  the  orbital  angular  momentum  of  each  subshell  from  the 
total  orbital  angular  momentum  and  the  intermediate  totals.  Finally,  we  uncouple 
the  last  electron  coordinate  Ig  from  the  total  for  that  subshell. 

For  the  case  in  which  a  and  a'  are  not  in  the  same  configuration,  we  may 
apply  the  methods  used  by  Cowan  to  calculate  general  transition  arrays.  This  is 
achieved  by  employing  a  process  similar  to  that  described  for  the  previous  case.  It 
is  a  tedious  but  tractable  process  of  accounting  for  coordinate  permutations  and 
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uncoupling  or  recoupling  angular  momenta  that  can  be  expressed  by  a  series  of 
factors,  which  Cowan  denotes  as  Di  through  D7.  One  may  refer  to  Appendix  C 
for  more  details.  Here,  we  show  only  an  abbreviated  form  similar  to  eq.  14.78  in 
Cowan  [54], 

{^a,j\\E^rCi'm^a',j')   =  D,  ■      ■  ■  ■  Dj  ■  {m,  mA"C^'^\\njJ,)  (3.40) 
j=i 

where 

<n'z'(^i)   =  {nuh\\A"\\n„l,).Af{x,) 

=  /rf:rra;XAW(^-^)/2„,,,(x,).  (3.41) 

The  indices,  i  and  j,  refer  to  the  subshells  relevant  to  the  one-electron  transitions. 
The  quantity,  .4|j'/„,,,(2;i),  in  eq.  3.41,  is  an  expression  for  the  matrix  element  of 
the  tth  term  in  the  multipole  expansion  of  the  interaction,  V^^\r,l),  integrated 
over  the  radiator  state  of  a  single  bound  electron.  Rn,i{xr)  is  the  relativistic  wave 
function  of  the  bound  electron  extracted  from  Cowan's  atomic  physics  code. 

At  this  point,  we  will  explore  the  details  of  evaluating  the  simplest  term  of  our 
electron  shift  and  width  operator,  B^^K 

3.4    Evaluation  of  B^^^ 

We  will  now  attempt  simplifications  of  the  expression  for  B^^^  made  possible  by 
the  trace  over  perturbing  electron  states.  Let  us  begin  with  eq.  3.37  and  substitute 
this  expression  into  B^^)  using  eq.  3.22,  disregarding  for  the  moment  the  term 
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acting  on  the  lower  states.  We  will  apply  the  Wigner-Eckart  Theorem  to  the 
matrix  elements  of  the  perturbing  electron  states, 


l,m  t,q 


^  —m  —qrrij 


X 


Nr 


Since  the  3-j  symbol  can  only  be  non-zero  when  mi  -  rui  +  q  =  0,  q  must  equal 
zero,  and  we  can  employ  Eq.  A.8  to  sum  over  m, 


E(-i)' 


~m+2l+t 


m   —m  0 


(3.43) 


This  expression  shows  that  the  index,  t,  of  the  multipole  expansion,  must  equal 
zero.  In  other  words,  only  the  monopole  interaction  between  the  radiator  electrons 
and  the  perturbing  electrons  contributes  to  the  first  order  shift.  We  can  now  re- 
express  B(^)  in  terms  of  the  reduced  matrix  element  of  the  monopole  interaction. 

'^Uj)i'^^,JM\t.A'\j)CTU)\'^a',J'M'){kA\<^^^  (3.44) 
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We  now  apply  the  Wigner-Eckart  theorem  to  the  matrix  elements  of  the  radiator 
states,  and  simplify: 


^   J     0    J'  ^ 


V 


-M  0  M' 


Nr 


TljJ){^^,j\\Y,A'\j)C^'\jm^,,j>){kJ\\A'',C^^^ 
We  use  the  identity  in  eq.  A. 5  to  simplify  the  expression  further, 


^    J     0  J'^ 
-M   0  M' 


{-iy-^6jj,SMM'[J]-'^',  (3.46) 


so  that 


41'/'   =   ^ff'^JJ'^MM'TiX'r  I  dhk^e-^'^Y^^to^' 


Nr 


(3.47) 


TU(J)(*a,j|lE^'°0)<^^°Hi)ll*a'.j)(A:,/||^';C'ri|A:,OT.-.(J). 


3=1 


We  must  now  consider  the  two  cases  for  which  we  calculate  the  reduced  matrix 
elements  over  the  radiator  states.  First,  when  the  matrix  elements  are  between 
states  in  the  same  configuration  we  will  define  an  operator  similar  to  Di  •  •  •  Dr, 


{^a,j\\^sA'\s)C^'\s)\\^^.,j,)  =  Dis){nMwsA'\s)C^'\s)\M.  (3.48) 
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Using  eqs.  3.47  and  3.48,  we  get  the  following  expression: 


t:'^Uj)D{s){nMw,A'\s)C^'\s)\\nJs){kJ\^^^^^  (3.49) 


5=1 


We  refer  to  Appendix  D  regarding  our  calculation  of  the  angular  term.  For 
the  case  when  a  and  a'  are  in  the  same  configuration,  we  arrive  at  very  simple 
solution,  Q  =  a',  and  the  factor  obtained  from  eq.  D.9  is  written  below: 

f  711/2 

^5(^)-^aa'^P^.  (3.50) 

If  we  now  factor  out  the  reduced  matrix  elements  of  C(°)(s)  and  C^^^l),  we 
can  use  the  identity  written  in  eq.  A. 12.  Also,  we  combine  A'°{s)  and  A'1  into 
'^nsU,n,i,i^^  1)>  as  it  was  defined  in  3.41.  So  we  can  write, 

Biflf   =   Sff>5aa'SjJ'^MM'StonXT  J  dkk'^e'^^^  x 

;,  ;  .      EWE"^.TtjJ)(A:,/|K,^„^Jl)||A;,/)T,,,(J')-  (3-51) 

Now  we  will  consider  the  second  case,  where  states  a  and  o.'  are  not  in  the  same 
configuration: 

B\%r   =   hj'SMM'StoSfrnX'r  fdkk'e-^''J2M!u^'' 

T^U'^)Di---Dr{nM\A''C^'^\\njlj){kJ\\A'^^^^^^  (3.52) 
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If  we  factor  out  the  reduced  matrix  element  (/i||C^°^||/j),  according  to  eq.  A. 12, 
we  find  that  Ij  =  Ij,  where  /j  is  the  orbital  angular  momentum  of  the  electron 
undergoing  the  transition  in  state  a  and  Ij  is  the  angular  momentum  of  the  same 
electron  in  state  a'.  Since  the  electron  must  be  in  the  same  angular  momentum 
state,  it  must  change  principle  quantum  number  in  order  to  undergo  a  transition. 
Since  all  the  other  radiator  electrons  do  not  undergo  transitions  a  and  a'  are 
members  of  a  Rydberg  series,  for  example,  3sls^  and  4sls^.  If  we  do  not  consider 
matrix  elements  where  electrons  undergo  transitions  between  manifolds  of  differ- 
ing principal  quantum  number,  only  the  diagonal  matrix  elements  remain.  The 
second  case  will  be  useful  for  examining  cross  terms  between  manifolds,  and  will 
be  examined  in  detail  in  Appendix  D.  For  the  hydrogenic-  or  near  hydrogenic  case, 
in  the  sense  that  there  may  be  no  more  than  one  electron  in  a  non-closed  subshell 
and  no  more  than  one  electron  in  the  upper  manifolds,  the  factors  Di  -  Dj  reduce 
to 

[7-11/2  9 

=  n  ^«»a;.  (3.53) 

Although  the  factor  [JY^'^/HY/'^  is  common  to  both  cases,  the  diagonal  terms  con- 
tain contributions  from  each  occupied  subshell,  while  only  the  subshells  containing 
the  electron  undergoing  the  transition  contribute  to  the  off-diagonal  terms.  We  can 
now  express  B^^)  in  a  more  familiar  form;  we  will  represent  the  perturber  states  in 
terms  of  Coulomb  wavefunctions,  Fi{r],  kxi),  where  /  is  quantum  number  represent- 
ing the  angular  momenta  in  the  partial  wave  expansion  and  k  is  the  wavenumber 
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of  the  perturber  momentum  [57].  Then, 


where 


Z  —  1 

a{l,  k)  =  argr(/  +  l  +  ir])  t]  =  — ;  .  (3.55) 

kao 


With  some  simplifications,  we  arrive  at  the  following  expressions  from  which  B^^^ 
was  calculated  in  this  chapter;  first, the  diagonal  term, 

BY/,f.   =   -nX^TT,T.^c.,a'St,oSff'    dke-^'\2l  +  l)x 
TtjJ)  I  rf:riFKr7,A;xi)^u;,4,^,„Jxi)F((77,A;xi)T,,,(J')  -By;^,,.  (3.56) 


Next,  the  cross  term, 

2 


41'/'    =   -"4EEM//'  n  /rffce-^*'(2/  +  l)x 

Q,a'  1=0  si^ii^i 

la(^)  / dx.Fiir),  kxMl,u,n,uMFi{r),  kx,)T^,,,{J)  -  Bfj,5i,  (3.57) 


where 


<^,„,/,(^l)     =     [dXrXlR*^^i^{Xr)(^--)Rn^,i.{Xr).  (3.58) 

\Xf       X\  / 


Xi  <  Xr 


We  have  arrived  at  useful  expressions  for  B^^^i  for  multi-electron  cases,  including 
a  cross  term  for  simple  cases.  There  are  two  results  one  should  remember  from 
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this  section.  First,  only  the  monopole  term  of  the  multipole  expansion  contributes 
to  the  first  order  shift,  and  only  if  the  plasma  electron  is  within  the  radiator 
orbital.  This  suggests  that  the  first  order  term  represents  an  average  shielding  of 
the  radiator  electron,  lowering  the  effective  ion  nuclear  charge,  Z.  Second,  there 
can  be  non-negligible  interactions  between  states  of  differing  principal  quantum 
number.  If  adjacent  levels  are  sufficiently  close,  the  shift  of  the  composite  lines 
may  not  be  adequately  described  by  the  linear  shifts  of  their  individual  energy 
levels. 

To  anticipate  the  relative  size  of  the  first  order  term,  one  can  calculate  the 
approximate  change  in  the  nuclear  charge  Z  experienced  by  the  radiator  electron 
as  it  is  shielded  by  plasma  electrons.  This  change  can  be  found  by  determining 
the  average  number  of  plasma  electrons  within  the  radiator  electron  orbital: 


(3.59) 


From  this  expression,  if  we  begin  with  the  expression  for  the  bound  energy  levels" 
and  obtain  a  difference  equation,  we  can  obtain  an  approximate  value  for  I\E^/E, 
due  to  the  first  order  shift: 


E 


Z^_Z^ 
nf  nj 

2Z/\Zi  _  2ZAZf 
nf  nj. 


(3.61) 


(3.60) 
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Using  eq.  3.59  and  the  expressions  above,  we  obtain  the  following  expression: 


This  rough  approximation  is  provided  as  an  aid  to  experimentalists.  It  suggests 
that  the  first  order  shift  is  strongly  dependent  on  principal  quantum  number,  as 
well  as  the  ion  charge  Z.  However,  the  dependence  on  Z  is  misleading  in  the  sense 
that  line  spectra  for  low  Z  ionic  radiators  are  usually  broadened  beyond  recognition 
before  an  appreciable  shift  can  be  observed. 


We  will  now  examine  A,t",i'(^i,  ^2),  which  contains  second  order  expression 
of  the  the  interaction  potential,  so  that  we  may  obtain  calculable  expressions  for 
5(2)  and  M{uj).  Employing  the  Wigner-Eckart  theorem,  we  use  eq.  3.37  and  the 
methods  introduced  in  Section  3.3  to  insert  the  expression  for  a  multi-electron 
matrix  element  from  eq.  3.34  into  the  expression  for  k-^  found  in  eq.  B.14: 

^1  '^^2  ^2) 


(3.62) 


3.5    Evaluation  of /)ii'/  i/(A;i,  A;2) 


E   E  E  E  E(-ir"^M^i,^i,mi|Al'Ci'^;(l)|A:2,/2,m2) 


a,a'  /l,mi  /2,"l2  <2,92 
a",a"' 


x(A;i,/i,mi|^t<^-^](l)l^2,/2,m2)- 


xTI,JJ)(*,,^^|X:A'(*'^(jOC(^)(j)|^,,,^,,^,,)T,",."(J'') 


(3.63) 
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Once  again  we  will  have  to  differentiate  between  cases  where  the  matrix  element 
is  between  states  of  the  same  configuration  and  those  that  are  not.  We  will  attempt 
to  simplify  the  angular  portions  as  we  did  for  B^^\  but  we  must  proceed  carefully. 
For  example,  the  ordering  of  the  matrix  indices  is  important,  and  the  angular 
momenta  for  state  a'  needs  to  maintain  the  same  relationship  with  a'"  as  state  a 
has  with  a",  so  indices  may  seem  to  be  reversed,  but  it  is  correct.  Since  we  will 
reverse  the  order  of  the  second  matrix  element,  we  must  remember  to  take  the 
complex  conjugate  of  the  factors.  This  will  be  seen  mostly  in  the  powers  of  -1, 
where  (-1)^  will  simply  be  replaced  by  (-1)"'^.  The  3-j,  6-j,  and  9-j  symbols  are 
real  valued,  so  we  only  need  to  be  concerned  with  having  the  correct  ordering. 

We  will  first  obtain  the  reduced  matrix  elements  of  the  perturbing  electron 
states, 


x{h,h,mM'iC^-i{m2,l2,m2r 
=    E  E  E  E(^i>^ill<C'^)||A:2,/2)(A:i,/i||AtC^^||A;2,/2)- 

h,mi  l2,m2  i2,92 


h  ti  I2 

-mi     -Qi  1712 


-mi    -q2  1712 


(.3.64) 


The  advantage  of  this  ordering  is  that  we  can  employ  the  sum  rule,  eq.  A.6,  by 
summing  over  mi  and  m2  without  reordering  the  expression. 


mim2 


h  ti  I2 

y  -mi    -qi   1712  J 


^    h       ^2      h  ^ 


-nil    -92  m2 


'5fit2«^gi92j^  (3.65) 
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We  now  know  that  ti  =  t2  and  qi  =  q^.  We  may  now  apply  the  Wigner-Eckart 
Theorem  to  the  matrix  elements  of  the  radiator  states, 


J=l  j=l 


=  (-1) 


J-M-J'+M' 


J       t  J" 

y-M    q   M"  j 


y  ~M'   q  M"  ^ 


Nr 


{^.M\T.A''  '\j)C^''\M^a'',J''){^a',J'\\  E  A'^'\f)C^Hfma"',J"r-  (3-66) 
J=l  j'  =  l 

Once  again  we  can  simplify  eq.  3.66  by  summing  over  q  and  M"  using  eq.  A. 7. 
This  is  possible  because  we  are  summing  over  states  i"  and  since  M"  is  a  good 
quantum  number  of  state  i".  Also,  both  the  reduced  matrix  elements  and  the 
transition  matrices  are  invariant  in  M".  We  will  refer  to  the  sum  over  i"  with  a 
prime  to  indicate  that  we  have  summed  over  M".  Referring  to  eq.  3.66,  we  arrive 
at  the  following: 


9  M" 


J       t  J" 
[-M    q    M"  J 


^  ^    J'      t     J"  ^ 


y  -M    q    M"  j 
which  leaves  us  with  the  following  expression: 


—  ^JJ'^MM'[J]  \ 

(3.67) 
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i" 

i"       a,c'    luh  t 
a", a'" 

x{k^,h\\A'\'C[''^\\k2,l2r 

><TU(J)(*a.^||E^'^'^^(i)<^^*'^0)ll*a'',^'')T.",a"(J'') 
j=l 

xTl„,,,„(J")(*a',j'||  E  A'^'\j')C^'Hj'mo''',j'nAJ')-  (3.68) 

j'=i 

We  cannot  go  any  further  with  this  expression,  unless  we  make  some  simplifying 
approximations.  We  will  attempt  to  reduce  the  transformation  matrices  between 
the  matrix  elements  of  the  interaction  potential  to  the  identity  by  summing  over 
the  energy  eigenstates  represented  by  i".  This  would  seem  simple  except  that  the 
expressions  for  5^2)  and  M{u)  in  eqs.  3.23,  3.29,  3.30,  3.31,  and  3.32,  are  dependent 
on  the  quantity  hui"  Let  us  make  the  assumption  that  the  states  are  degenerate 
within  a  given  manifold.  This  is  a  standard  assumption  as  long  as  level  splitting  due 
to  the  microfield  is  small  relative  to  the  energy  spacing  between  manifolds  as  well  as 
small  relative  to  the  plasma  frequency.  This  approximation  is  equivalent  to  setting 
Aw"  equal  to  Aa;'  +  a;„"„/,  where  uJn"n'  is  the  average  energy  difference  between 
the  n"th  and  n'th  manifolds.  We  will  also  make  the  assumption  that  Ta,i{J) 
is  block  diagonal  in  the  sense  that  the  cross  terms  of  Ta,i{J)  between  states  of 
differing  manifolds  should  be  negligible.  This  is  reasonable  since  the  basis  states 
are  determined  by  the  atomic  physics  code  in  such  a  way  that  the  states  of  differing 
manifolds  are  orthogonal  regardless  of  the  angular  momentum  representation  [54]. 
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This  allows  us  to  substitute  hui't  with  huii'  +  hun"n',  which  commutes  with  the 
block-diagonal  transformation  matrices.  We  sum  Ta",i"{J")Tl„^^,„{J")  over  i"  to 
obtain  6a"a"'-  We  now  rewrite  A,i",t'{A;i,  with  just  two  transformation  matrices 
bracketing  the  reduced  matrix  elements.  We  may  now  take  advantage  of  the  sum 
over  a". 

i" 

=  hj-^MM-  E  E  E  E[^'  JV{k,,  h\\A!^i^c?\\k,,  h) 

a,Q'  h,l2  t 

x(*a',^||  E  A'^'Hj')C^'\j'ma",j'-r  T,,AJ')  (3.69) 
j'=l 

As  we  now  evaluate  the  angular  portions  of  the  matrix  elements,  we  will  pro- 
ceed carefully,  considering  several  cases.  This  problem  is  described  in  detail  in 
Appendix  E.  For  a  multi-electron  atom  with  no  more  than  two  electrons  having 
the  same  principal  quantum  number  and  no  interactions  between  states  of  differing 
principal  quantum  number,  we  arrive  at  the  following  solution  when  states  a,  a', 
and  a"  are  in  the  same  configuration.  The  sum  over  the  states  a"  that  are  in  the 
same  configuration  as  the  two  other  states  is  indicated  by  a  prime. 

a"   s=l  r=l 

+2  E  E  ^-^f^^D'is,  r)w.w,A'^l^,^  ^t.n.,  (l.  \\C^'^  I  lU  (M \C^^^  \\lr) 
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This  expression  is  diagonal  for  the  multi-electron  cases  considered  and  is  similar  to 
the  diagonal  term  of  5^^^  in  the  sense  that  we  sum  over  the  contributions  from  each 
subshell.  Although  this  expression  appears  rather  formidable,  it  is  analogous  to 
the  Law  of  Cosines;  refer  to  eqs.  E.18  and  E.20.  If  one  considers  the  monopole  term 
of  the  multipole  expansion,  t  =  0,  D'{s,r)  reduces  to  (-!)'»+''•  and  the  equation 
above  can  be  expressed  as  a  square  of  the  sums.  We  arrive  at  the  following  solution 
for  the  second  case  involving  transitioning  matrix  elements  from  eq.  E.70. 


where  the  sum  is  executed  over  all  configurations,  except  the  configuration  of  state 
a,  as  opposed  to  individual  states  with  differing  intermediate  angular  momentum 
sums.  This  simplifies  the  problem  tremendously.  The  quantity  la  is  the  orbital 
angular  momentum  of  the  electron  undergoing  the  transition  while  it  is  in  state 
a,  and  n  is  the  occupation  number  of  the  subshell  containing  the  transitioning 
electron  in  state  a  and  k  is  the  occupation  number  of  the  subshell  containing  the 
transitioning  electron  in  a".  The  matrix  elements  of  all  electrons  not  undergoing 
transitions  reduce  to  delta  functions,  and  we  need  only  to  sum  over  the  angular 
momentum  of  the  electron  undergoing  the  transition  in  state  a".  We  now  express 
Di,i",i'{ki,k2)  in  terms  of  Coulomb  wavefunctions  [57],  and  express  the  reduced 
matrix  elements  (/||C(')||/')  in  terms  of  3-j  symbols.  The  following  expression  is 
the  sum  of  three  terms.  The  first  two  are  obtained  summing  over  the  states  a" 
which  are  in  the  same  configuration  as  states  a  and  a',  while  the  third  term  is 


a"  config' 


[J] 


(3.71) 
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obtained  by  summing  over  all  other  states  a". 

(^1 '  ^2)  =  A,i'  (^1,^:2) 


<^aa'  E  E  E  E 


a  t    ^  r=l 


[t] 


0     0  0 


V 


J  0 


0    0  0 

2 

Ti',a(J) 


+  ^^^'EEE;^E_E  K^.)^^"^^^|^'^^^''' 


a  li,l2    t    "  s=lr=s+l 


( 


D'{s,r) 


h  t  I2 
0   0  0 

Ad 


0   0   0  ^ 

TU(J)  l^''dx,Ft,{rjuk,x,)A'il^(x,)Ft,{ri2,k2X^) 


J 


\ 


0   0  0 


/  "  dX2Fi,  (772,  k2X2)Alll^  {X2)FI,  (771 ,  ^1X2) 


T.',a(J) 


+  '^-'E  E  EE^(n-^) 


"  config"h,/2  « 


[/l,  ^2,  ^q"] 
[t] 


2 

(     i2  ^ 

0  0  ^ 

0   0  , 

T.',a(J)  (3.72) 


We  have  now  obtained  the  factor  A;2)  which  contains  all  the  information 

regarding  the  second  order  expression  of  the  radiator-perturbing  electron  interac- 
tion, which  is  used  in  B^^)  ^nd  the  real  and  imaginary  parts  of  M{u)).  It  is  diagonal 
for  the  multi-electron  cases  under  consideration.  It  is  useful  for  He-like  and  Li-like 
satellites  of  K-shell  spectra  as  well  as  Li-like  and  Be-like  L-shell  spectra.  It  has 
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been  reduced  to  a  sum  over  configurations  with  special  consideration  for  the  term 
when  the  matrix  elements  act  between  the  same  configuration.  For  a  hydrogenic 
radiator,  n  and  k  and  Wr  reduce  to  one  and  there  is  no  sum  over  subshell,  thereby 
eliminating  the  second  term.  Therefore,  eq.  3.72  reduces  to  the  result  obtained  by 
Woltz  [13,  14]. 


CHAPTER  4 
RESULTS 


In  this  section,  we  will  discuss  the  results  of  our  calculations  using  the  full- 
Coulomb,  multi-electron  formalism.  We  will  show  that  significant  line  shifts  and 
asymmetries  are  revealed  by  employing  the  full-Coulomb  expression  of  the  radiator- 
perturbing  electron  interaction  potential.  Also,  we  will  compare  our  results  with 
other  models  such  as  a  multi-electron  formalism  employing  the  dipole  approxima- 
tion and  a  semi-classical  model  that  contains  all  orders  in  the  radiator-perturbing 
electron  interaction  potential.  Finally,  we  will  address  the  impact  the  line  shift  on 
line  merging,  especially  the  importance  of  including  interactions  between  states  of 
differing  principal  quantum  number,  the  so-called  "mixing"  terms. 

4.1    Computational  Aspects 

In  order  to  obtain  computed  lineshapes,  a  great  deal  of  information  concerning 
the  relativistic  atomic  physics  of  the  radiator  must  be  obtained.  The  existing 
multi-electron  line  shape  code  MERL  [15,  16]  uses  as  input,  energy  levels  and 
reduced  matrix  elements  of  dipole  and  quadrupole  interactions  generated  by  the 
suite  of  atomic  physics  codes  written  by  Robert  Cowan  at  Los  Alamos  National 
Laboratory.  In  addition,  the  author  of  this  manuscript  also  extracted  relativistic, 
multi-electron  radial  wavefunctions  for  the  bound  electron  states  as  well  as  mixing 
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coefficients  between  the  LS  representation  and  energy  eigenbasis.  These  coefficients 
form  the  basis  of  our  transformation  matrices,  Ti^a{J). 

The  reduced  dipole  matrix  elements  obtained  from  the  atomic  physics  had 
previously  been  used  to  calculate  the  electron  broadening  operator  and  are  still 
employed  to  obtain  the  radiator-ion  microfield  interaction.  These  dipole  values, 
should  not  be  confused  with  the  the  dipole  terms  in  the  multipole 
expansion  of  the  radiator-perturbing  electron  interaction,  which  account  for  the 
possibility  of  the  plasma  electron  penetrating  the  radiator  orbital.  The  radial 
wavefunctions  are  written  out  by  the  code  as  a  set  of  functions  for  each  occupied 
subshell,  with  the  intention  of  using  them  as  product  wavefunctions.  The  Coulomb 
wavefunctions  are  obtained  using  the  Numerov  method  to  solve  the  appropriate 
differential  equation. 

4.2    Shifts  and  Asymmetries 

In  this  section,  we  will  present  the  results  of  our  calculations,  focusing  on  the 
first  order  shift  B^^).  From  Eqn.  3.17,  we  know  that  B^^^  is  linear  in  density  and  is 
also  dependent  on  temperature,  principal  quantum  number  n  and  orbital  angular 
momentum  /.  In  Figure  4.1  we  see  the  first  order  shift,  B^^\  as  a  function  of 
temperature  for  At+16  immersed  in  a  plasma  (any  kind,  it  is  not  dependent  on 
concentration)  characterized  by  an  electron  density,  A^^  =  1  x  lO^'^cm"^.  In  the 
subsections  that  follow,  we  consider  various  factors  governing  the  nature  of  the 
calculated  shift. 
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Figure  4.1:  Comparison  of  level  shifts,  first  order  in  the  interaction  potential,  of 
Helium-like  Ar  as  a  function  of  temperature;  electron  density,  1  x  lO^^cm"^.  Note 
the  increase  in  the  magnitude  of  the  level  shifts  with  increasing  principal  quantum 
number  and  dependence  on  angular  momentum  state. 
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4.2.1  Dependence  on  principal  quantum  number,  n 

It  is  apparent  from  Figure  4.1  that  the  shift  is  strongly  dependent  on  the 
initial  states  of  the  radiator,  particularly  the  principal  quantum  number,  n.  In 
physical  terms,  this  can  be  understood  by  remembering  that  a  higher  shell  presents 
more  volume  to  be  penetrated  by  plasma  electrons.  Also,  for  a  given  principal 
quantum  number,  subshells  with  lower  angular  momentum  have  a  greater  volume, 
which  explains  why  the  s  orbital  experiences  the  greatest  shift.  These  results  are 
consistent  with  the  difficulty  experienced  measuring  the  shift  of  lines  originating 
from  more  tightly  bound  levels,  such  as  the  a-  or  /?-lines,  where  the  calculated 
shift  can  be  less  than  the  instrumental  resolution.  The  plasma  spectroscopy  group 
at  University  of  Florida  first  noticed  an  experimentally  measurable  shift  of  a  He- 
7  line  of  Ar+16;  it  was  this  observation  that  subsequently  inspired  this  study. 
Another  benefit  of  the  strong  shift  dependence  on  the  principal  quantum  number 
is  the  ability  to  measure  differential  shifts  between  members  of  a  Rydberg  series. 
In  the  measurements  that  we  will  discuss,  time-resolved  x-ray  spectrometers  were 
not  absolutely  calibrated,  and  the  author  has  had  to  rely  on  evaluating  the  shift 
of  one  line  relative  to  another. 

4.2.2  Dependence  on  electron  density,  rie 

As  we  stated  earlier,  the  shift  and  the  width  are  linear  in  electron  density.  As 
a  result,  the  profiles  both  shift  and  broaden  with  increasing  electron  density.  The 
/?-  and  7-lines  in  Figure  4.2  are  calculated  using  the  full-Coulomb  expression  of 
the  radiator-perturbing  electron  interaction.  The  dipole  approximation  is  used  for 
the  same  lines  in  Figure  4.3. 
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Figure  4.2:  Comparison  of  broadened  line-shapes  using  a  full-Coulomb  expression 
for  the  radiator-perturbing  electron  interaction,  a)  The  Ar-l-16  He-^  is  shown  as 
a  function  of  density,  with  a  temperature,  kT  =  1000  eV.  The  increasing  density 
causes  both  broadening  and  shifting  of  the  lines;  b)  The  Ar-|-16  He-7  line  is  shown 
for  the  same  plasma  conditions  as  a). 
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Figure  4.3:  Line-shapes  of  He-like  Ar  calculated  in  the  dipole  approximation  of 
the  radiator-perturbing  electron  interaction,  are  shown  as  a  function  of  electron 
density,  with  a  temperature  kT  =  1000  eV.  There  are  no  shifts  and  the  lineshapes 
appear  broader  than  those  shown  in  Figure  4.4  for  the  same  densities  .  a)  The 
Ar+16  He-/3  is  shown;  b)  The  Ar+16  He-7  is  shown. 
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The  effect  of  shifts  increasing  with  both  increasing  electron  density  and  prin- 
cipal quantum  number  can  cause  theoretical  difficulties  when  calculating  lines  in 
a  Rydberg  series.  One  can  anticipate,  at  sufficiently  high  electron  densities,  that 
higher  lying  lines  will  shift  past  their  lower  n  neighbors:  of  course  this  result  is 
based  on  the  shift  having  a  linear  density  dependence.  However,  since  the  mag- 
nitudes of  non-zero,  off-diagonal  matrix  elements  of  the  first-order  shift,  B^^\  are 
on  the  order  of  diagonal  elements,  we  must  consider  the  mixing  of  adjacent  level 
when  making  our  shift  calculations.  For  example,  the  matrix  element,  B^p^jp,  is 
about  half  that  of  B^iJ^^p.  Preliminary  calculations  have  shown  that  mixing  between 
states  within  a  Rydberg  series,  due  to  interactions  with  perturbing  electrons  and 
ions,  removes  the  possibility  of  higher  spectral  line  passing  lower  series  members. 
Rather  the  series  members  tend  to  'pile  up'  against  one  another  as  the  plasma 
density  increases.  This  phenomenon  will  be  addressed  in  subsequent  sections. 

4.2.3    Dependence  on  temperature,  kT 

Also,  from  Figure  4.1,  one  can  see  the  shift  is  weakly  dependent  on  the  plasma 
temperature.  This  has  led  some  to  use  an  "ion-sphere"  or  "uniform  electron"  model 
to  calculate  the  shift,  which  assumes  a  uniform  distribution  of  electrons  within  the 
ion  sphere  as  a  high  temperature  limit  [19,  36].  However,in  Figure  4.4,  we  see 
line  asymmetries  which  become  more  pronounced  as  the  temperature  decreases. 
Figure  4.4  shows  line  profiles  of  the  He-/?  and  He-7  of  Ar-f  16  as  a  function  of 
temperature  for  A^e  =  1  x  lO^^cm-^  The  shifted  profile  shapes  are  more  strongly 
dependent  on  temperature  than  line  profiles  calculated  in  the  dipole  approximation, 
as  seen  in  Figure  4.5. 
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Figure  4.4:  Comparison  of  broadened  line-shapes  using  a  full-Coulomb  expression 
for  the  radiator-perturbing  electron  interaction,  a)  The  Ar-M6  He-^  is  shown  as 
a  function  of  temperature,  with  an  electron  density,  A^^  =  1  x  lO^'^cm^^  The  shift 
in  the  full-Coulomb  calculation  causes  increasing  asymmetry  as  the  temperature 
decreases;  b)  The  Ar-|-16  He-7  line  is  shown  for  the  same  plasma  conditions  as  a). 
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Figure  4.5:  Line-shapes  of  He-like  Ar  calculated  in  the  dipole  approximation  of  the 
radiator-perturbing  electron  interaction,  are  shown  as  a  function  of  temperature, 
with  an  electron  density,  A^e  =  1  x  lO^^cm-^  There  are  no  shifts  or  asymmetries  in 
the  line-shapes,  and  the  temperature  dependence  is  not  as  pronounced  as  is  seen 
in  Figure  4.4.  a)  The  Ar-l-16  He-/?  is  shown;  b)  The  Ar-l-16  He-7  is  shown. 
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4.2.4  Dependence  on  orbital  angular  momentum,  / 

The  differing  shifts  of  the  angular  momentum  states  observed  in  Figure  4.1 
can  result  in  a  distinct  asymmetry  of  the  line  profile.  The  calculated  line-shape 
asymmetries  arise  from  the  splitting  of  angular  momentum  states,  due  primarily  to 
the  interaction  with  the  ion  microfield.  The  /?-  and  7-lines  in  Figure  4.4  represent 
the  ls3p  —  Is^  and  ls4p  —  Is^  transitions,  respectively;  however,  the  p  states  are 
mixed  with  other  angular  momentum  states  by  Stark  splitting  of  the  line.  Each 
Stark  component  shifts  differently  according  to  its  constituent  angular  momentum 
states.  However,  sum  rules  regarding  the  dipole  oscillator  strength  must  still  be 
obeyed.  As  a  result,  the  centers  of  mass  of  the  line  profiles  in  Figure  4.4  shift  at  the 
rate  of  the  Isnp  -  Is^  transition,  but  the  "peaks"  shift  less,  reflecting  the  mixing 
with  the  higher  angular  momentum  states. 

4.2.5  Dependence  on  ion  charge,  Z 

Finally,  the  shift  is  very  weakly  dependent  on  the  radiator  ion  charge  Z.  The 
shift  increases  with  decreasing  ion  charge;  in  other  words,  a  /?-line  from  Helium-like 
Ar  will  shift  slightly  more  than  the  the  Ly-p  of  Ar.  This  property  justifies  compar- 
ing shifts  between  members  of  the  same  Rydberg  series  as  opposed  to  attempting 
to  observe  a  shift  between  lines  from  the  same  transition  but  differing  ionization 
stages.  For  example,  observing  the  Ar  Ee-/3  line  and  its  Li-like  satellites  has  been  a 
useful  diagnostic  of  temperature  and  density  in  laser  implosion  experiments;  how- 
ever, the  complex  of  lines  all  shift  approximately  the  same  amount.  Therefore, 
since  absolute  calibration  of  spectral  data  was  not  available,  a  bulk  shift  of  the  line 
complex  might  be  interpreted  as  calibration  error. 
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4.3    Comparison  with  Other  Models 

In  this  section,  we  will  compare  our  results  with  those  based  on  other  models. 
We  find  that  the  first  order  shift  that  we  calculate  compares  well  with  two  sets  of 
calculations  based  on  impact  models  that  include  the  effects  of  the  monopole  term 
of  the  radiator-perturbing  electron  interaction.  Also,  we  compare  our  results  with 
a  semi-classical,  all-order  model  and  find  good  agreement  within  the  regions  of  the 
validity  of  the  second  order  approximation.  In  Figure  4.6,  the  current  model  is 
compared  with  calculations  by  Griem  et  al.  [23]  and  Nguyen al.  [24].  Griem  et 
al.  used  a  distorted  wave  method  to  determine  the  phase  shift  of  perturbing 
electron  due  to  the  penetrating  monopole  interaction  with  the  radiator.  Since 
they  employed  a  Born  approximation,  theirs  is  a  first  order  result  and  agrees  very 
well  with  the  first  order  result  of  the  current  model.  Nguyen  et  al.  used  a  quantum 
mechanical  impact  model  similar  to  the  current  model. 

In  the  electron-density  and  temperature  regime  under  consideration,  the  second 
order  theoretical  results  were  corroborated  by  independent  calculations  using  an 
all-order,  semi-classical  model  [58,  59].  Figure  4.7  shows  Lyman-/?  and  Lyman-7 
profiles  at  a  typical  temperature,  1000  eV  and  electron  density,  Ne=^lx  lO^^cm"^, 
which  corresponds  to  plasma  conditions  observed  in  ICF  experiments.  Ones  sees 
good  agreement  between  the  two  models,  although  the  second  order  calculation  has 
a  slightly  larger  shift.  Figure  4.8  shows  Lyman-/?  and  Lyman-7  lines  at  a  lower 
temperature,  600  eV,  A^e  =  1  x  lO^'^cm-^  Even  here  the  models  have  not  yet  begun 
to  diverge  significantly.  However,  we  can  notice  minor  differences  in  the  shifts  and 
the  width  is  slightly  greater  in  the  second  order  calculation.  Finally,  we  compare 
calculations  of  the  Ly-^  line  at  kT=300eV  in  Figure  4.9  in  Figure  4.9.  Here  we 
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Figure  4.6:  Comparison  of  shifts,  first  order  in  the  interaction  potential,  of  Helium- 
like Ar  with  a  shift  calculations  by  Griem  et  al.  (1990)  and  Nguyen  et  al.  (1986); 
electron  density,  1  x  lO^^cm-^  There  is  excellent  agreement  with  Griem  and  good 
agreement  with  Nguyen. 
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Figure  4.7:  Comparison  of  Ar+17  Lyman  series  lineshapes,  calculated  using  the 
second-order,  quantum  mechanical  formalism  and  an  all-order,  semi-classical  ex- 
pression, for  a  temperature  1000  eV  and  an  electron  density,  ^  1  x  lO^^cm-^ 
Notice  that  there  is  very  good  agreement,  a)  The  Ar+17  Ly-/?  is  shown;  b)  The 
Ar+17  Ly-7  is  shown. 
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begin  to  see  differences  between  the  two  models.  According  to  the  conditions  set 
forth  by  Cooper,  Kelleher,  and  Lee  [30],  the  temperature  limit  for  a  second  order 
calculation  would  be  and  460  eV  and  260  eV  for  the  Lyman  j3-  and  7-lines  of 
Ar+17,  respectively,  indicating  that  the  apparent  difference  is  consistent  with  this 
validity  condition.  As  was  stated  earlier,  this  may  confirm  that  strong  collisions 
have  become  more  important  at  lower  temperatures.  A  strong  collision  is  defined 
as  one  where  momentum  transfer  is  large  with  respect  to  the  initial  momentum 
of  the  colliding  electron.  It  may  seem  counter-intuitive  that  penetrating  electrons 
can  undergo  weak  collisions  but  one  must  realize  that  the  momentum  transfer 
is  dependent  on  the  duration  of  the  collision  as  well  as  on  the  proximity  of  the 
perturber.  Thus,  sufficiently  fast  electrons  can  still  experience  weak  collisions.  We 
expect  that  the  second  order  approximation  will  be  valid  under  the  conditions 
specified  by  CKL  [30].  It  is  unlikely  that  one  can  observe  or  7  lines  of  K-shell 
Argon  at  temperatures  below  600  eV,  but  they  might  be  observed  in  absorption. 

Although  the  second  order  approximation  is  sufficient  for  isolated  lines  of  highly 
ionized  radiators  under  most  conditions  where  such  lines  would  be  observed,  there 
are  imaginable  cases  where  the  second  order  approximation  is  not  valid.  It  has  been 
stated  that  the  second  order  approximation  breaks  down  for  frequency  separations 
from  line  center  that  are  much  greater  than  the  plasma  frequency.  Lines  emitted 
from  low-Z  radiators,  such  as  Hell,  can  have  widths  much  larger  than  the  plasma 
frequency.  This  condition  also  occurs  when  attempting  to  calculate  adjacent  lines 
in  a  Rydberg  series  simultaneously.  Such  a  calculation  is  necessary  when  studying 
line  merging.  We  will  discuss  such  calculations  in  a  later  section  under  conditions 
where  they  may  seem  inappropriate;  however  we  must  remember  that  the  plasma 
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Figure  4.8:  Comparison  of  the  Ar+17  Lyman  series  lineshapes,  calculated  using 
the  second-order,  quantum  mechanical  formalism  and  an  all-order,  semi-classical 
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expression,  for  a  temperature  600  eV  and  an  electron  density,  A^g  =  1  x  10 
As  the  temperature  decreases  there  is  still  good  agreement,  but  they  do  begin  to 
diverge,  a)  The  Ar+17  Ly-p  is  shown;  b)  The  Ar+17  Ly-7  is  shown. 
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Figure  4.9:  Comparison  of  the  Ar+17  Ly-/?  lineshape,  calculated  using  the  second- 
order,  quantum  mechanical  formalism  and  an  all-order,  semi-classical  expression, 
for  a  temperature  300  eV  and  an  electron  density,  A^e  =  1  x  lO^^cm"^.  We  now 
begin  to  see  significant  differences. 
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frequency  increases  with  electron  density  and  our  model  may  still  be  appropriate 
at  these  ultra-high  densities. 

4.4    Experimental  Analysis 

Using  theoretical  line-shapes  which  include  shifts,  significant  improvements 
have  been  achieved  in  interpreting  experimental  data.  Time-resolved  Ar  K-shell 
(3-4.5keV)  spectra  were  obtained  during  microballoon  implosions  performed  using 
the  OMEGA  laser  system  at  the  Laboratory  for  Laser  Energetics  at  the  University 
of  Rochester  [60].  20  micron  thick  CH  microballoons  of  440  micron  inner  radius 
were  filled  with  20  Atm  of  DD  doped  with  0.25%  Ar.  The  microballoons  were 
imploded  by  20-25  kJ  of  ultraviolet  energy  distributed  symmetrically  among  60 
beams.  Time-resolved  Ar  x-ray  line  spectra  provide  signatures  of  plasma  conditions 
during  the  implosion.  The  spectra  were  recorded  using  a  streaked  spectrograph 
with  a  RbAP  crystal  dispersing  the  spectrum  onto  a  250  Angstrom  Au  photo- 
cathode.  Spectral-  and  temporal  resolution  were  approximately  £'/A£'  ^  500  and 
At  ^  20  -  30ps,  respectively.  The  resulting  images  were  recorded  on  film,  then 
digitized.  Verifying  the  presence  of  shifts  first  requires  careful  determination  of 
isolated-ion  line  position,  the  digitized  images  were  corrected  for  film  sensitivity 
and  known  streak  camera  induced  distortions  including  streak  angle  and  curvature 
of  iso-temporal  lines  [61].  The  dispersion  relation  was  determined  by  the  line  po- 
sitions at  early  time,  when  the  low  densities  assured  negligible  shifts.  This  same 
dispersion  relation  was  then  used  for  all  times.  The  theoretical  line  positions  were 
those  obtained  using  Cowan's  [54]  atomic  physics  suite  as  input  into  our  spectral 
line  profile  calculation.  Figure  4.10  shows  time  resolved  spectra  for  a  single  implo- 
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sion  where  the  isolated  ion  line  positions  have  been  indicated  with  vertical  lines. 
The  line  shifts  of  the  7-lines  are  clearly  visible. 

Before  making  comparisons  with  experimental  data,  the  theoretical  spectrum 
was  modified  to  account  for  1)  shell  opacity,  2)  streak-camera  filtering,  3)  Au 
photocathode  response,  and  4)  the  instrumental  response  function  due  mainly  to 
source  size  broadening.  No  arbitrary,  wholesale  shift  of  the  theoretical  spectrum 
or  individual  lines  was  permitted.  In  order  to  analyze  the  data,  theoretical  spectra 
were  calculated  using  the  following  lines:  Ar  He-^  (ls3/  -)■  Is^),  He-7,  Ee-6,  Ly-/3, 
Ly-7>  Ly-<^  and  the  attendant  satellites  of  the  He-/?,  He-7,  and  Ly-p.  It  is  impor- 
tant to  include  satellites  since  their  appearance  in  the  spectra  can  be  incorrectly 
interpreted  as  broadening  and  shifting  of  the  resonance  lines.  Also  included  were 
the  effects  of  ion  dynamics,  ion-quadrupole  interactions,  Doppler-,  opacity-,  and 
instrumental  broadening. 

In  order  to  save  computer  time,  theoretical  lineshapes  were  calculated  in  the 
dipole  approximation  for  the  electron  width,  with  the  first  order  shift  included  in 
the  calculation.  This  approximation,  when  compared  with  our  more  exact  cal- 
culations, has  been  shown  to  be  adequate  for  data  analysis  under  the  conditions 
specified.  This  method  works  for  two  reasons.  First,  it  has  been  shown  that  using 
the  dipole  approximation  for  the  radiator-perturbing  electron  interaction  gives  a 
good  approximation  for  the  width  even  for  cases  where  the  inter-electron  spac- 
ing is  comparable  with  the  average  bound  electron  radius.  [13,  14]  Secondly,  we 
have  found  that  the  first  order  term  of  the  electron  shift  B^^)  provides  the  dom- 
inant contribution  to  the  shift.  In  Figures  4.11  and  4.12  we  show  a  comparison, 
respectively,  of  the  Ar+17  Ly-/5  and  Ly-7,  -5  and  -e  lines,  calculated  with  the 
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Figure  4.10:  Time-resolved  lineouts  from  LLE  shot  10778  (.125%  Ar  in  D2).  Ver- 
tical lines  represent  isolated  ion  line  positions  of,  going  from  left  to  right,  Ar  He-^, 
He-7,  Ly-^,  and  Ly-7.  Note  how  the  7  lines  shift  more  than  the  lines.  Time 
progresses  in  the  upward  direction. 
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full-Coulomb  model  of  the  electron  shift  and  width  operator  with  interactions  be- 
tween the  upper  manifolds,  and  a  profile  calculated  in  the  dipole  approximation 
for  the  electron  width  with  a  first  order  shift  included,  but  no  interactions  between 
the  upper  manifolds.  We  can  see  that  with  a  typical  resolution  and  for  a  high 
plasma  density,  the  full-Coulomb  calculation  of  the  Ly-f3  line  produces  a  slightly 
narrower  profile  than  does  the  dipole  approximation;  also  it  has  a  little  greater 
shift.  At  these  electron  densities,  the  /3-\me  can  be  satisfactorily  approximated 
by  the  dipole  approximation  with  the  first  order  shift.  However,  even  though  the 
7-line  experiences  a  great  deal  of  interaction  with  the  6  and  e  lines,  causing  it  to 
be  broader  than  it  would  appear  in  an  isolated  calculation,  the  dipole  calculation 
provides  a  noticeably  different  profile:  the  contrast  of  peak  height  to  wing  intensity 
is  substantially  greater. 

We  have  shown  lines  for  He-like  radiators,  but  not  for  more  complex  lineshapes, 
which  are  used  in  the  analysis  of  K-shell  Ar  spectra.  Figure  4.13  compares  the 
Ar+^^  Ee-p  and  its  Li-like,  n  =  2  satellite  calculated  in  the  full-Coulomb  formal- 
ism with  lineshapes  calculated  in  the  dipole  approximation.  We  can  see  that  the 
lineshapes  of  the  satellites  have  little  structure,  particularly  with  respect  to  the 
primary  resonance  line.  It  has  been  known  that  the  dipole  calculation  provides  a 
good  approximation  for  the  width.  While  here  it  may  be  possible  to  shift  the  entire 
lineshape  calculated  in  the  dipole  approximation  by  an  amount  equal  to  the  first 
order  shift,  for  lineshapes  with  more  structure  it  may  be  necessary  to  include  the 
detailed,  angular  momentum-dependent  shift  into  the  initial  calculation,  so  that 
the  asymmetries  apparent. 


Figure  4.11:  This  figure  compares  the  Ly-/?  of  1%  Ar+i^  in  D2,  where  the  effects 
due  to  plasma  electrons  are  calculated  in  the  dipole  approximation,  with  a  first 
order,  linear  shift  included,  and  with  a  full-Coulomb  expression  of  the  radiator- 
plasma  electron  interaction  potential,  with  interactions  between  the  n=3,4,5  and 
6  manifolds.  The  lineshapes  were  calculated  for  plasmas  characterized  by  a  tem- 
perature, kT=1000eV  and  electron  density,  =  1  x  lO^^cm-^  The  profiles  were 
convolved  with  a  Gaussian  with  fwhm  corresponding  to  an  instrumental  resolution 
E/l^E  =  500. 
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Figure  4.12:  This  figure  compares  the  Ly-7  through  -e  of  1%  Ar+^^  in  D2,  where 
the  effects  due  to  plasma  electrons  are  calculated  in  the  dipole  approximation, 
with  a  first  order,  linear  shift  included,  and  with  a  full-Coulomb  expression  of 
the  radiator-plasma  electron  interaction  potential,  with  interactions  between  the 
n=3,4,5  and  6  manifolds.  The  lineshapes  were  calculated  for  plasmas  characterized 
by  a  temperature,  kT=1000eV  and  electron  density,  rie  =  1  x  lO^'^cm^^  The  pro- 
files were  convolved  with  a  gaussian  with  fwhm  corresponding  to  an  instrumental 
resolution  E//\E  =  500. 
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Figure  4.13:  This  figure  compares  the  Ar-^^^  Li-like,  n  =^  2  satellite  of  the  He-/?. 
The  temperature,  kT=1000eV  and  electron  density,  iVg  =  1  x  lO^^cm-^  a)  These 
lines  were  calculated  in  the  dipole  approximation;  b)  The  full-Coulomb  expression 
of  the  radiator-perturbing  electron  interaction  was  used.  One  can  see  there  is  little 
difference  in  the  width  of  the  satellite  lines. 
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In  Figure  4.14,  spectroscopic  data,  lineout  170  from  Omega  shot  10778,  is  shown 
together  with  the  'best  fit'  unshifted-  and  shifted-theoretical  profiles,  where  only 
temperature  and  density  are  adjustable  parameters.  A  model  spectrum  character- 
ized by  the  temperature  kT  =  1250  eV  and  electron  density  Ne  =  1.5  x  lO^^cm"^ 
was  found  to  be  the  best  fit  to  the  data.  It  should  be  noted  that  because  the  a,  P, 
and  7  lines  shift  by  significantly  different  amounts,  one  cannot  shift  the  theoretical 
spectrum  by  an  arbitrary  amount  and  fit  the  entire  spectrum.  In  particular,  the 
region  containing  the  Ar  He-7  and  Ly-/?  lines  provides  an  excellent  test  of  any  the- 
ory of  shifts,  as  experimental  uncertainties  are  unlikely  to  vary  appreciably  over 
this  200  eV  spectral  range.  Figure  4.14  illustrates  the  size  of  the  shift  relative  to 
the  widths  and  demonstrate  a  good  agreement  of  theory  to  experiment  when  the 
shifts  are  included.  Also,  time  resolved  data  from  Nova  shot  26031209  is  shown 
in  Figure  4.15,  while  in  Figure  4.161ineout  e  from  Nova  shot  26031209  is  shown 
together  with  the  'best  fit'  unshifted-  and  shifted-theoretical  profiles  for  a  temper- 
ature of  kT  =  1000  eV  and  electron  density  iVg  =  6  x  lO^^cm^^  Notice  how  the 
shifted  Ly-/3  line  picks  up  the  asymmetry  of  the  spectral  data. 
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Figure  4.14:  Detail  from  lineout  170  from  shot  10778  (.125%  Ar  in  D2)  and  fits 
using  model  spectrum,  for  =  1.5  x  10'^'^  cm.-^ ,  kT  =  1250  eV.  Edge  effects  on 
the  data  due  to  the  photocathode  can  be  seen  just  to  the  right  of  the  Lyman- 
7  line,  where  the  spectrum  seems  to  drop  off  rapidly.  The  model  spectrum  has 
been  convolved  with  a  gaussian  with  a  fwhm  of  approximately  6-7  eV  to  account 
for  instrumental  broadening,  a)  Unshifted  lineshapes  were  used  for  the  model 
spectrum;  b)  Shifted  lineshapes  were  used  for  the  model  spectrum. 
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Figure  4.15:  Time-resolved  lineouts  from  Nova  shot  26031209  (.2%  Ar  in  D2). 
Vertical  lines  represent  isolated  ion  line  positions  of,  going  from  left  to  right,  Ar  He- 
P,  He-7,  Ly-p.  Note  how  the  7  lines  shift  more  than  the  p  lines.  Time  progresses 
in  the  upward  direction. 
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Figure  4.16:  Detail  from  lineout  e  from  shot  Nova  shot  26031209  (.2%  Ar  in  D2)  and 
fits  of  the  Ar  He-7  and  Ly-/3  using  model  spectrum,  for  rig  =  6  x  102^cm-^  kT  = 
1000  eV.  The  model  spectrum  has  been  convolved  with  a  Gaussian  with  a'fwhm 
of  approximately  6-7  eV  to  account  for  instrumental  broadening,  a)  Unshifted 
lineshapes  were  used  for  the  model  spectrum;  b)  Shifted  lineshapes  were  used  for 
the  model  spectrum. 
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4.5    Line  Merging  and  Continuum  Lowering 

At  very  high  electron  densities,  >  1  x  lO^^/cc,  higher  lying  series  members  of 
Ar+17  experience  plasma  induced  shifts  greater  than  those  of  more  tightly  bound 
lower  series  members,  allowing  for  the  possibility  for  line  to  shift  past  its  neighbor. 
However,  as  the  densities  increase,  the  splitting  of  adjacent  energy  levels  due  to 
the  ion  microfield,  together  with  the  line  shift,  allows  for  the  eventual  overlap  of 
states,  and  requires  these  previously  isolated  states  be  treated  as  composite  states. 
Mathematically,  the  magnitudes  of  non-zero  matrix  elements  of  the  first-order  shift, 
due  to  mixing  between  states  of  differing  principal  quantum  number,  are  on 
the  order  of  diagonal  elements.  For  example,  the  element,  BjJJgp,  is  about  half 
that  of  54p|4p.  We  will  show  that  the  inclusion  of  mixing  between  states  within  a 
Rydberg  series  produces  a  non-linear  shift  theory:  this  in  turn  tends  to  cause  the 
lines  to  'pile  up'  against  one  another  during  the  line  merging  process.  In  order  to 
demonstrate  the  effects  mixing  between  upper  manifolds,  we  must  calculate  the 
lineshapes  over  very  large  frequency  ranges  with  respect  to  the  plasma  frequency. 
It  has  already  been  stated  that  the  second  order  calculation's  validity  fails  for 
frequency  separations  greater  than  two  times  the  plasma  frequency.  Therefore, 
one  must  understand  that  these  results  offer  only  a  qualitative  understanding  of 
line  merging;  however,  at  some  point  the  levels  become  sufficiently  close  to  allow 
the  accurate  calculation  of  spectral  line  originating  from  these  composite  states. 
In  Figure  4.17,  we  show  the  Lyman  series  of  Ar+i^  from  the  p-line  to  the  e-line 
as  higher  lying  series  members  are  successively  added  to  the  calculation  of  the 
composite  line  profile;  the  electron  density  is  increased  from  1  x  lO^^cm"^  to  3  x 
lO^^cm-^  and  5  x  lO^^cm-^. 
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Figure  4.17:  These  figures  compare  the  Ly-^  through  -e  of  1%  Ar+i^  in  D2,  as  higher 
series  members  are  added  to  the  calculation  of  a  composite  line  shape  using  the  full- 
Coulomb  model  and  accounting  for  interactions  between  states  of  differing  principal 
quantum  number.  The  lineshapes  were  calculated  for  plasmas  characterized  by  a 
temperature,  kT=1000eV  and  electron  density,  a)  rie  =  1  x  lO^-^cm-^-  b)  n  = 
3  X  1024  cm-3;  c)  ne  =  5  X  10^4  cm-^ 
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In  Figure  4.18  we  show  the  full-Coulomb  calculation  of  the  Ar+^^Ly-/5  through 
-e  with  and  without  interactions  between  upper  manifolds  at  electron  densities, 
ne  =  lx  10^4  cm-^  and  rie  =  5  x  lO^^'cm-^  In  Figures  4.19  and  4.20  are  shown 
details  of  the  individual  lines  at  a  the  lower  density.  The  line  widths  and  shifts 
appear  to  be  reasonably  approximated  by  a  linear  model  at  the  lower  density; 
however,  the  peak  of  the  Ly-7  experiences  an  accelerating  shift  and  an  increase 
in  its  intensity  relative  to  the  Ly-p.  Currently,  there  is  not  data  available  at 
sufficiently  high  densities  to  confirm  this  behavior. 

We  also  show  a  comparision  in  Figure  4.21  with  the  dipole  calculation  with 
a  first  order,  linear  shift  included.  Once  again,  the  approximation  is  adequate 
for  the  lower  density,  but  fails  dramatically  at  the  higher  density.  This  clearly 
demonstrates  the  limitations  of  the  dipole  approximation. 
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Figure  4.18:  These  figures  compare  the  Ly-p  through  -e  of  1%  Ar+i^  in  D2,  cal- 
culated with  and  without  mixing  between  upper  manifolds,  the  former  of  which 
can  result  in  non-linear  shifts  at  sufficiently  high  densities.  The  profiles  were 
convolved  with  a  gaussian  with  fwhm  corresponding  to  an  instrumental  reso- 
lution E/AE  =  500.  The  lineshapes  were  calculated  for  plasmas  characterized 
by  a  temperature,  kT=1000eV  and  electron  density,  a)  =  1  x  lO^^cm-^; 
b)  We  =  5  X  10^4  cm-^;  differences  are  clearly  apparent  at  this  resolution  for  the 
higher  density.  The  lineshapes  calculated  with  the  linear  model  were  displaced 
vertically  by  an  arbitrary  amount  for  display  purposes. 
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Figure  4.19:  This  figure  compares  the  Ly-/3  of  1%  Ar+i^  in  D2,  where  the  effects  due 
to  plasma  electrons  are  calculated  in  the  full-Coulomb  expression  of  the  radiator- 
plasma  electron  interaction  potential,  with  and  without  interactions  between  the 
n=3,4,5  and  6  manifolds.  The  lineshapes  were  calculated  for  plasmas  characterized 
by  a  temperature,  kT=1000eV  and  electron  density,  Ue  =  1  x  lO^^cm-^  The  pro- 
files were  convolved  with  a  gaussian  with  fwhm  corresponding  to  an  instrumental 
resolution  E/AE  -  500. 
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Figure  4.20:  This  figure  compares  the  Ly-7  through  -e  of  1%  Ar+^^  in  D2,  where 
the  effects  due  to  plasma  electrons  are  calculated  in  the  full-Coulomb  expres- 
sion of  the  radiator-plasma  electron  interaction  potential,  with  and  without  in- 
teractions between  the  n=3,4,5  and  6  manifolds.  The  lineshapes  were  calculated 
for  plasmas  characterized  by  a  temperature,  kT=1000eV  and  electron  density, 
rie  =  1  X  10^^  cm"^.  The  profiles  were  convolved  with  a  gaussian  with  fwhm  corre- 
sponding to  an  instrumental  resolution  E//S.E  =  500. 
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kT=1000eV;  N,=5e24/cc 
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Figure  4.21:  These  figures  compare  the  Ly-p  through  -e  of  1%  Ar+i^  in  D2,  where 
the  effects  due  to  plasma  electrons  are  calculated  in  the  dipole  approximation, 
with  a  first  order,  linear  shift  included,  and  with  a  full-Coulomb  expression  of  the 
radiator-plasma  electron  interaction  potential,  with  a  non-linear  shift.  The  pro- 
files were  convolved  with  a  gaussian  with  fwhm  corresponding  to  an  instrumental 
resolution  E/AE  =  500.  The  lineshapes  were  calculated  for  plasmas  character- 
ized by  a  temperature,  kT^lOOOeV  and  electron  density,  a)  rig  =  1  x  lO^^cm"^; 
b)  ne  =  5  X  lO^^cm"^;  it  is  clear  the  approximations  are  no  longer  valid  at  the 
higher  density.  The  lineshapes  calculated  in  the  dipole  approximation  were  dis- 
placed vertically  by  an  arbitrary  amount  for  display  purposes. 


CHAPTER  5 
CONCLUSION 


In  this  work,  we  presented  a  generalized  full-Coulomb  formalism  of  line  broad- 
ening for  multi-electron  ions,  where  the  line  shifts  arise  consistently  from  the  relax- 
ation theory.  We  compared  our  results  with  other  available  models  and  compared 
model  spectra  calculated  with  the  line  shift  with  experimental  data.  Certain  sim- 
plifying approximations  were  tested  that  could  save  computer  time  and  man-hours 
spent  applying  more  involved  models  when  it  may  not  be  necessary.  Finally,  we 
considered  the  impact  of  the  shift  on  line  merging  and  proposed  that  the  shift 
may  have  a  measurable  effect  on  plasma  phenomenon  dependent  on  energy  level 
structure. 

We  discovered  that  incorporating  the  dense  plasma  line  shift  into  line  spectra 
calculations  can  have  a  significant  eflfect  on  the  position  and  overall  shape  of  line 
profiles  from  highly  ionized,  moderate  Z  radiators. 

From  a  strictly  applied  point  of  view  this  is  immediately  useful.  We  found 
that  including  the  shift  in  line  spectra  from  highly  ionized,  moderate  Z  radiators 
has  improved  plasma  diagnostics.  We  have  shown  the  lineshapes  with  the  shift 
are  more  temperature  sensitive  than  those  calculated  in  the  dipole  approximation, 
which  is  also  an  important  issue  in  diagnostics.  We  then  showed  that  employing 
the  dipole  approximation  with  the  first  order  shift  is  an  acceptable  alternative  to 
the  full-Coulomb  calculation  for  highly  ionized,  moderate  Z  radiators,  in  cases 
where  the  second  order  approximation  is  valid  and  the  orbital  size  is  less  than  or 
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approximately  equal  to  the  inter-electron  spacing,  particularly  for  the  computer 
intensive  multi-electron  calculations.  These  results  are  significant  because  plasma 
diagnostics  may  be  extended  to  greater  densities  as  researchers  attempt  to  achieve 
"break  even"  in  fusion  research,  as  well  as  being  applied  to  astrophysical  studies. 

If  we  consider  a  wider  perspective,  we  can  speculate  on  the  impact  that  the 
shifts  will  have  on  plasma  phenomena  dependent  on  the  energy  level  structure  of 
radiators  as  they  interact  with  hot,  dense  plasmas.  Now  that  we  demonstrated 
that  including  mixing  between  upper  states  of  differing  principal  quantum  number 
introduces  non-linear  behavior  of  the  shift  as  adjacent  members  of  a  Rydberg 
series  merge,  we  must  redefine  the  distinction  between  continuum  and  bound  states 
and  among  bound  states  themselves.  Ionization  balance  is  sensitive  to  changes  in 
energy  level  structure  [33,  34,  35,  36,  37],  and  opacity  is  sensitive  to  the  energy 
location  of  line  emission.  Therefore,  understanding  the  effect  of  the  shift  near  the 
continuum  edge  is  very  important.  The  shift  may  have  a  measurable  effect  on 
continuum  lowering,  ionization  balance,  and  opacity  calculations,  as  well  as  line 
ratio  diagnostics. 

There  are  areas  that  need  further  examination.  First,  as  the  nuclear  charge  Z 
decreases,  the  first  order  shift  becomes  smaller  relative  to  the  width.  Higher  order 
terms  have  greater  significance.  We  did  not  explore  this  in  detail;  on  the  other 
hand,  the  observed  shifts  are  small  and  apparently  have  little  eflFect  on  line  merging. 
Second,  at  this  time,  we  cannot  not  easily  compare  our  line  merging  results  with 
those  of  an  all-order  calculation,  which  would  be  a  more  appropriate  choice  for 
studying  adjacent  lines.  Such  a  comparison  would  be  quite  interesting.  Third, 
other  concerns  arise  as  the  density  increases,  the  binary-collision  assumption  and 
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the  effect  of  plasma  electron  correlations  has  been  over-simplified,  we  may  need  to 
consider  how  these  matters  are  treated  in  condensed  matter  physics.  Fourth,  the 
plasma  electrons  do  not  act  in  isolation,  we  have  commented  little  about  the  effects 
of  the  plasma  ions  in  conjunction  with  the  electrons.  Finally,  at  this  time  we  do  not 
have  experimental  evidence  to  confirm  the  line-merging  behavior  predicted  by  our 
theory.  This  is  due  to  our  inability  to  generate  sufficiently  high  plasma  densities 
and  temperatures  in  the  laboratory.  The  few  experiments  where  such  densities 
were  observed  were  focused  on  the  a-lines  which  would  have  been  unaffected  by 
line  merging  in  that  plasma  regime  [62]. 

Future  experiments  are  planned  in  which  we  expect  to  observe  sufficiently  high 
densities  and  temperatures,  thereby  affording  us  an  opportunity  to  verify  our  pre- 
dictions. 


APPENDIX  A 

3N-J  SYMBOLS  AND  IRREDUCIBLE  TENSOR  OPERATORS 


This  appendix  provides  brief  overview  of  a  formalism  that  is  very  useful  for  the 
treatment  of  angular  momentum  dependent  matrix  elements.  These  definitions, 
identities,  and  sum  rules  will  be  employed  throughout  this  work.  First,  we  will 
introduce  a  Wigner  3N-j  symbols  [54,  63,  64]  in  terms  of  an  integral  over  spherical 
harmonics, 
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There  are  several  useful  relations  and  symmetry  properties  that  will  be  employed 
in  this  work,  such  as 
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The  following  relation  can  be  obtained  using  the  orthogonality  of  spherical  har- 
monic functions, 


j  f  0 
m   —m'  0 
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(A.5) 


Also,  the  equation  mi  +  m2  +  mg  =  0  must  be  satisfied  for  a  3-j  symbol  to  be 
nonzero.  Below  we  state  several  useful  sum  rules: 
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where  5{jij2jz)  is  +1  if  the  triangle  relations  are  satisfied. 

Also,  we  will  be  making  extensive  use  of  the  Wigner-Eckart  Theorem  [55,  56], 
which  is  written  below.  We  will  employ  the  theorem  to  express  the  matrix  elements 
of  our  interaction  potential  in  terms  of  reduced  matrix  elements,  as  follows: 


(aim|T«|aym')  =  (-1)^- 
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where  T^^*^'  is  an  irreducible  tensor  operator,  j  and  /  are  the  angular  momentum 
indices  of  the  coordinates  relevant  to  the  operator,  and  a  and  a'  represent  all  the 
other  quantities  (energy,  spin,  etc.)  that  define  their  respective  states. 

A  particular  example  of  reduced  matrix  elements  that  will  be  employed  is  the 
matrix  elements  of  spherical  harmonics.  We  will  be  writing  the  angular  portion 
of  our  radiator-perturbing  electron  interaction  in  terms  of  normalized  spherical 
harmonics,  defined  as. 

Several  relations  pertaining  to  reduced  matrix  elements  of  spherical  harmonics  are 
written  below: 
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and 


(/'||CW||0  =  (-l)^/||C(^)||O.  (A.13) 


We  will  also  use  6-j  and  9-j  symbols  to  evaluate  the  angular  portions  of  multi- 
electron  matrix  elements.  The  symbols  are  defined  as  follows: 
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We  will  make  extensive  use  of  the  following  relations  when  one  of  the  arguments 
is  zero  as  well  as  orthogonality  relations.  With  regard  to  symmetry  properties,  a 
6-j  symbol  is  unchanged  when  any  two  columns  are  interchanged,  as  well  as  if 
any  two  numbers  in  the  bottom  row  are  interchanged  with  the  corresponding  two 
numbers  in  the  top  row.  The  arguments  in  the  top  row  must  satisfy  a  triangle 
relation. 


1^2  -  h\  <li<l2  +  k,      Ih  -  /al  <l2<li  +  k,      Ih  -  kl  <l3<li+  h,  (A.16) 

and  their  sum  be  a  non-negative  integer.  Also,  any  set  of  arguments  that  can  be 
permuted  into  the  top  row  must  also  satisfy  these  conditions. 
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The  orthogonality  relation  for  6-j  symbols  can  be  written  as  follows: 
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We  will  also  make  use  of  the  following  sum  rule  in  section  E.l  of  Appendix  E, 
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The  following  relation  is  employed  extensively  throughout  Appendices  D  and  E, 
when  one  of  the  arguments  of  a  6-j  symbol  is  zero, 
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Returning  to  9-j  symbols,  one  must  multiply  the  symbol  by  a  factor  of  (-1)^, 
where  K  is  the  sum  of  all  the  arguments  when  exchanging  a  pair  of  columns  or 
rows.  The  arguments  in  any  row  or  column  must  satisfy  a  triangle  relation  and 
their  sum  must  be  a  nonnegative  integer.  The  orthogonality  relations  are  written 
as  follows: 
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The  following  relation  is  used  when  one  of  the  arguments  is  zero, 
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APPENDIX  B 
DETERMINING  THE  FORM  OF  B('\  B^^\  AND  M{uj) 

In  the  following  chapter,  we  will  obtain  forms  of  the  various  terms  of  the  electron 
shift  and  width  operator  in  terms  of  their  matrix  elements.  These  terms  of  the 
electron  shift  and  width  operator  have  not  previously  been  presented  together  in 
a  detailed  manner,  and  we  intend  to  show  that  our  results  have  been  obtained  in 
a  consistent  fashion. 

B.l    Matrix  Elements  of  B^^^ 

In  the  series  of  steps  below,  the  Liouville  operator  representing  B^^^  is  expressed 
as  a  commutator  to  get  the  two  terms  shown.  Complete  sets  of  states  are  inserted 
to  obtain  the  correct  matrix  elements  acting  on  an  arbitrary  operator  dependent 
on  radiator  coordinates,  Cl,  as  follows, 

n^[TriLi(r,  l)/o(l)]^^;^v'fiMV 

nJ2[{^ia\V^{r,l)fo{l)C^\lya)  -  {fia\fo{l)m{r,l)\ua)] 
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n  E  [{^oc\Vi{r,l)f,{l)\fi'a'){ii'a'\Q\va) 
-  {fxa\fo{l)\fi'a'}{nW\QV,{r,  l)\ua)], 
n    E    KHViir,  l)MlWa'){fiW\{l\i^'a"){u'a"\iya) 

fi'  i/'aa'a" 

-  {m\fo{Wa'){fxW\Q\u'a''){\u'a''\Vi{r,  (B.l) 
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The  operator,  Q,  does  not  act  on  perturbing  electron  states,  therefore  its  matrix 
elements  are  diagonal  in  the  perturbing  electron  states.  Also,  the  distribution 
function  for  the  perturbing  electron  states,  /o(l),  will  be  diagonal  in  both  radiator 
and  perturber  electron  states,  so  that  we  may  simplify  this  last  expression, 

S  ^ilU'^^''"'     =    "      S  [{f^(^\Vl{rAWoi)fo{a)Cl^>^,5^^>5a'a"Sa"a 
H'l/'  n'u'aa'a" 

-   k{oi)^ti'u'{y'oi\Vi{r,l)\ua)5f,^,,5aa'6a'a"]-  (B.2) 
Hence,  expression  for  B^^}^,^,  can  be  written  as, 

^'iv  =  ^E/o(c^)  1)1  A)^^..'  -  {^'ci\V,{r,  \)\va)6^^.] .  (B.3) 

a 

In  order  to  calculate  this  quantity,  we  will  express  the  perturbing  electron  states 
in  terms  of  a  wavenumber,  {a  k),  representing  solutions  to  the  eigenvalue 
problem  of  free  electrons  in  the  presence  of  an  ion.  Referring  to  Eq.  3.4,  we  define 
our  basis  set  in  terms  of  i/(l), 

We  then  express  the  sum  over  a  as  an  integral  over  k,  and  arrive  at  the  following 
expression  for  B^^\ 


5(1)  =  n  Jdkfik)  [{i^;k\V,{r,  l)|/i'; k)<5,,,  -  O.T.] ,  (B.6) 
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where  the  other  term  (O.T.)  includes  matrix  elements  between  statesz^  and  u'. 
We  shall  make  the  assumption  that  the  plasma  electrons  are  in  local  thermal 
equilibrium  and  employ  a  Maxwellian  distribution  for  /o(l),  so  that 

/o(l)  =  A^e    2-  ,      where     Ay  =^  f    ^    J   .  (B.7) 
We  substitute  these  expressions  into  eq.  B.6  to  obtain  the  following: 

5(1)  =nA^ I (ike-^  [(//; k|Fi (r,  1)1//'; k)5,,,  -  O.T.],  (B.8) 
so  that  equation  B.8  will  be  a  starting  point  to  evaluate  B^^)  in  section  3.4. 

B.2    Matrix  Elements  of  S^^) 

We  will  follow  a  similar  method  as  was  used  in  section  B.l  to  obtain  B^'^'i  in 
terms  of  matrix  elements  of  Vi(r,  1).  We  therefore  return  to  eq.  3.18, 


B(')  =  -nTn 


^i(^,l)/o(r)/o(l)  /  rfre^[^«+^(i)Vie-^[^M+^(i)l/o(r)- 

(B.9) 

Using  the  same  notation  as  in  section  B.l,  we  re-express  Li(r,  1)  as  a  commutator, 
we  get  the  following  expression: 
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E^U'^M'.'   =  -n5:/^dr(^a|Fi(r,l)/o(r)/o(l)e^W'-H^(i)l 

xK(r,l)e-^['^(^)+'^(i«/o(r)-ifi|i/a) 
+  rrfr(/ia|/o(r)/o(l)e^[^W+^(i)]K(r,l)e-^I^('-)+^(^^] 

^h{r)-'(lV,{r,l)\va).  (B.IO) 

Next,  as  we  did  for  B^^\  we  will  insert  complete  sets  of  radiator  and  perturber 
states  and  express  B^^)  in  terms  of  the  matrix  elements  of  the  interaction  Fi(r,  1), 
so  that  we  can  execute  the  integral  over  r, 

=   -"EE  r^^(/^«l^i(^l)l/^V)/o(//")/o(a')e^K"+-o'I 

/i"  aa'  ■'^ 
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xe--K'+-.']/o(^')-i(^V|Fi(r,  (B.ll) 

Once  again,  express  the  perturber  states  a  and  a'  in  terms  of  the  wave-numbers 
ki  and  k2  respectively,  and  substitute  the  sums  over  a  and  a'  with  integrals, 

=  ^^'i  / ^k2E(/^ki|V^i(r,l)|/i"k2)(/i"k2|v;(r,l)|/z'ki) 

li" 

xfo{fJ'")fo{k2)M^iT'      dre^K"+^-'-^'-^]5^^,  +  q.T.  (B.12) 
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We  now  execute  the  integral  over  r  to  obtain  the  following  expression: 


^(2) 


=   -n/c^ki|dk2  5;(A^ki|yi(r,l)|/x%)(/k2|K(r,l)|/x'ki)/o(/i'') 
x/o(^2)/o(//)-^^  ^,  U.u'+O.T.  (B.13) 

In  order  to  simplify  eq.  B.13,  we  define  the  expression,  ^2),  which 

contains  all  the  information  concerning  the  second  order  radiator-perturbing  elec- 
tron interaction.  We  also  express  the  perturber  states  in  terms  of  a  partial  wave 
expansion  where  /  and  m  describe  the  angular  momentum  state,  and  k  is  the  wave 
number,  so  that  we  arrive  at  the  following  expression: 

00       h       00  I2 

D,y>yih,k2)    =   kjklYl    E    E    E  {f^hJi,m,\V,{r,l)\tx"k2j2,m2) 

li=Omi  —  ~li  I2—O  Tn2=-l2 

I2,  m2\Vs{r,  l)\n'ki,  /i,  mi).  (B.14) 

We  will  now  substitute  Boltzmann  distributions  for  /o  (//'),  /o(/^")  and  /o(A;2)  and 
then  substitute  D^f,»^.{ki,  ki)  into  eq.  B.13, 


=   -S^^mX^  J  dkx  j  dki      D^,^,,^,  {ki ,  ^2) 


+0.r.,  (B.15) 


where  the  other  term  (O.T.)  contains  the  matrix  elements  between  u  and  u'  as  well 
as  //  and      which  is  often  called  an  interference  term. 
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The  expression  D^^>i^>{ki,k2)  is  common  to  both  B^^^  and  M{u)  and  will  prove 
to  be  the  most  challenging  portion  to  simplify  and  evaluate.  We  now  proceed  to 
finding  an  expression  for  M{uj)  in  terms  of  matrix  elements  of  Vi(r,  1). 

B.3    Matrix  Elements  of  M{uj) 

Returning  to  eq.  3.8,  we  see  that  M{u)  is  already  at  least  second  order  in  the 
radiator-perturbing  electron  potential  £(r,  1).  We  should  replace  [uj-L{r)-L{l)]~^ 
for  [L{r)  +  L{1)  +  C{r,  in  order  to  keep  the  interaction  to  second  order. 
Instead,  we  will  replace  £(r,  1)  with  (Li(r,  1)),  or  B'-^^  This  is  necessary  because 
the  first  order  shift  can  be  too  large  to  treat  perturbatively.  We  therefore  obtain 
the  following  expression 

M{u)  =  nTri{£(r,  l)[a;  -  L{r)  -  L(l)  -  B^'^]'' f{rl)/:{r,  l)f{r)-'},  (B.16) 

which  can  be  expressed  as  a  Laplace  transform, 

M{uj)  =  -mTn  |£(r,  1)  e'(-'-i{r)-L(i)-B(i)+u+)t f(^rl)Cir,  l)/(r)-^| . 

(B.17) 

Now  we  will  expand  out  the  Liouville  operators  in  terms  of  commutators  in  the 
next  four  equations:       .  '    ■  .      .  <    T  ^ 

fi^v  =  lim  -in  /  dt  e'^'^^''+^^ 
EiH^r,  l)e-(^W+^(iH^<^^)V(rl)£(r,  l)f{r)-'Q\i.a),  (B.18) 


which  leads  to  the  following  expression  as  one  commutator  is  expanded, 
E  =  -in  [  (ite'('^+^^+)* 

a 

-    {na\C{r,  l)e-(^W+^(i)+«'^')/(rl)/(r)-ij^K(r,  (B.19) 

We  now  expand  the  time  development  operators  as  well  as  the  other  Liouville 
operators  to  arrive  at  the  following: 

5]  M(a;)^^^v'f^MV'  =  lim  -in  f  dte'^'^+''+^^y 
n'v'  J  a 

-    (/ia|e-(^('-)+^(i)+«<'^)7(rl)K(r,  l)/(r)-ii)e^(^W+«(iH«^^^)*V,(r,  l)\ua) 

+   (//a|e-W'-H^(i)+«<^')7(rl)/(r)-ir2K(r,  l)eW)+//(iHB(i))t^^(^^  ^^^^^^ 

(B.20) 

In  the  next  steps,  we  will  introduce  complete  sets  of  states  to  obtain  an  expres- 
sion for  M(a;)  in  terms  of  matrix  elements  of  1): 
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(/.a|K{r,  l)e-^(^('-)+^(i)+^'^')7(rl)K(r,  l)/(r)-i| 

-  (//aie-W'-)+^(i)+^*'')7(rl)K(r,l)/(r)-i|/iV)(A'|f^k'/?') 

(z.'/3'|e'(^W+^(i)+^'")'\/(r,l)|m) 

-  (/.a|V^,(r,l)e-W'-)+^(i)+^''')7(rl)/(r)-VV)(A'|i7|^/'/?') 

(i^'/?'|K(r,l)e'^''^'"^+''^'^+^'''^>a) 
+  (/za|e-W'-)+^(i)+^''')7(rl)/(r)-VV)(//V|Q|^.'/?') 

(i.'/?'|K(r,  l)e'(^W+^(i)-^«'''Hy,(r,  l)|j.a)}  ,  (B.21) 

We  will  now  add  another  set  of  states  to  the  first  and  last  terms  to  obtain 

M(a;)^^^v'  =  ^Hm  -in  j  dt  e'^'^+''+^^  J2 

aa' 

E  {m\Vs{r,  l)|/xV')(/xV'|e-^(^('-)+^(^H^'^^)7(rl)K(r,  l)f{r)-'\f^W) 

fi"a" 

-(//a|e-'W'-H^W+^<^')7(rl)K(r,l)/(r)-VV)(i.V|e'(^^^^^ 
-(^«|F,(ra)e-W)+^(i)+^^'')7(rl)/(r)-i|A0^'a'|K(r,l)e^W^^ 
+    E  (/^«k-'(^(^^+^(^)+^*'')7(rl)/(r)-VV) 

i/"q" 

^V|K(r,  l)e'(^('-H^(i)+^^'')>V')Ka"|K(r,  l)|m).  (B.22) 
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We  will  assume  H{r),  H{1),  and  B^^\  along  with  the  corresponding  distribution 
functions,  are  diagonal  in  the  radiator  and  perturber  states.  With  some  simplifying 
and  renaming  of  indices  in  the  next  few  sets  of  equations,  we  obtain  the  following 
expression: 

MMmV'   =    hm  -my"rfte'('^+"+)'5; 

aa' 

{Y:{HV,{r,  l)|/.V)(//V|K(r,  l)|/i'«)/(//V)/(Mr' 

i{E^<-E^n+E^-E^,+B[\^-B^'])t. 

-  {m\Vs{r.  l)Wa'){u'a'\Vs{r,  l)Wcx) f  {fia) f  {fx')-' 

-  {^ia\V,{r,  l)\^i'a'){u'a'\Vs{r,  l)Wa)f{fx' a') 

i{E,-E.+Ea-E^,+Bl'^-B^]^)t 


+  EW^\Vs{r,Wa'){u"a'\Vs{r,l)\i.a)f{^a)f{fi) 


-1 


AB,,-E^^E^,-E.+Bl]]-Bl}^)t^^^^^^  ^g_23) 


The  first  and  last  terms  represent  the  second  order  shift  and  width  operators  from 
upper  and  lower  states,  respectively.  The  middle  terms  are  interference  terms. 
The  calculational  details  will  be  similar,  so  we  will  focus  on  the  first  term.  We 
will  express  the  perturbing  electron  states  in  terms  of  the  partial  wave  expansion 
and  re-express  the  matrix  elements  of  the  interaction  potential  in  terms  of  the 
expression  for  D^^'(ki,k2)  in  Eq.  B.14,  where  Au^>,  =  oj  -  E^,  +  E„  Finally, 
we  employ  the  Boltzmann  distribution  functions  to  obtain  final  our  expression  for 
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M{lo)  as  follows, 


/f  k^h^ 
dki  I  dk2Df,^>{ki,k2)e''^^^^ 

rffe'(^'^'''-+^-^+^^  (B.24) 


APPENDIX  C 
GENERAL  TRANSITION  ARRAYS 


In  this  appendix,  we  shall  describe  in  more  detail  the  functions  of  operators 
Di  through  D^  shown  in  Equations  3.40  and  3.41.  More  information  can  be  found 
in  Section  14-11  of  Cowan's  book  The  Theory  of  Atomic  Structure  and  Spectra, 
1981.  [54]  First,  we  must  introduce  the  labeling  of  angular  momentum  states  in  the 
geneological  CS  basis  with  the  following  notation: 

{\{lT(^iLxSulTc^2L2S2)C2S2,  lToc^L^S,\C^S^, . .  .}C,S,J,M,,  (C.l) 

where  is  the  orbital  angular  momentum  of  a  single  electron  in  the  subshell  m, 
Wm  is  the  occupation  number  of  that  subshell,  represents  all  other  quantum 
numbers  of  that  subshell,  and  q  is  the  total  number  of  relevant  subshells.  Lm  is 
the  sum  of  orbital  angular  momentum  of  all  the  electrons  within  subshell  m,  with 
the  corresponding  value  for  spins  5^.  The  intermediate  sum  of  the  total  angular 
momentum  up  to  and  including  subshell  mis  Cq,  with  a  corresponding  value  for 
spin.  Finally,  the  quantum  numbers  representing  the  total  angular  momentum, 
which  couples  £,  and  5„  are  J,  and  Mq. 

Most  of  the  factors  described  in  this  appendix  are  derived  when  coupled  angular 
momentum  states  are  uncoupled  in  order  to  evaluate  the  matrix  element  of  a  single 
angular  momentum  coordinate.  We  often  will  be  applying  the  identities  found  in 
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eqs.  11.38  and  11.39  in  Cowan,  and  they  are  written  below, 


(«Ui«2i2j||TW(l)||aij>^i^/) 
{^ihc^2mW^'\2)\\a\j[a',j'^3') 


3'    k  j[ 


{a^n\\T('\l)\\a[j[)  (C.2) 


Ji  32  3 
k    3'  32 


{c^232\\W^'\2)\\a',j',),  (C.3) 


where  T^''\l)  and  W^(*^)(2)  are  irreducible  tensor  operators  acting  on  differing  co- 
ordinates of  a  coupled  angular  momentum  state.  The  indices  ji  and  j[  represent 
the  angular  momentum  states  of  the  first  coordinate,  with  similar  notation  for  the 
second  set.  The  total  angular  momentum  is  represented  by  j  and  /.  All  other 
quantum  numbers  of  the  coordinates  are  represented  by  ai  and  a2.  These  expres- 
sions will  be  employed  throughout  Appendices  D  and  E.  A  brief  explanation  of  6-j 
symbols  can  be  found  in  Appendix  A. 

In  calculating  a  general  transition  array,  we  assume  that  it  is  a  one-electron 
transition,  in  other  words,  an  electron  moves  from  one  occupied  subshell  in  the 
first  state  to  another  subshell  in  the  final  state,  while  all  other  electrons  remain 
in  the  same  subshells  during  the  transition.  We  employ  the  convention  that  an 
electron  moves  from  a  "lower"  subshell  to  a  "higher"  subshell.  We  can  find  the 
reverse  transition  by  evaluating  the  complex  conjugate  of  the  matrix  element. 


('^,|V|^,>  =  ((^,|V|^,))*. 


(C.4) 
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We  must  introduce  some  notation  that  will  be  used  to  describe  a  "general 
transition."  A  one  electron  transition  can  be  described  in  its  most  general  sense 
as  follows: 

/-I . .  ./^  . .  ./fc-i . . ^      .  ..in-i  ...Ik.  ..iw,^  ^^^^ 

where  represents  the  orbital  angular  momentum  of  an  electron  occupying  the 
mth  subshell,  Wm  is  the  occupation  number  of  the  m  subshell,  q  is  the  total  number 
of  relevant  subshells.  If  one  were  considering  a  hydrogenic  atom,  n  -  k  =  1,  and  all 
other  occupation  numbers  would  equal  zero.  This  notation  will  be  used  throughout 
Appendices  D  and  E. 

The  factor  Di  arises  from  coordinate  permutations  as  well  as  addressing  cases 
where  the  occupancy  of  a  relevant  subshell  is  greater  than  one.  This  term  is  analo- 
gous to  ELi  t^s---  for  matrix  elements  between  states  with  the  same  configuration, 

(C.6) 

where 

Ap  =  k-1+   Yl  ^m,  (C.7) 

m=i+l 

The  factor  D2  arises  from  coefficient  of  fractional  parentage  expansions  when 
the  transition  involves  subshell  occupied  by  more  than  two  electrons.  Coefficients  of 
fractional  parentage  arise  when  the  angular  momentum  of  a  transitioning  electron 
is  uncoupled  from  the  angular  momenta  of  the  other  electrons  occupying  the  same 
subshell  as  the  transitioning  electron.  For  the  case  when  only  one  or  zero  electrons 
shares  a  subshell  with  the  transitioning  electron,  the  sum  of  the  remaining  angular 
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momenta  is  trivial,  so  such  a  factor  reduces  to  one. 


D2  =  (/raii:iS,{|/ria:L:s:)(/J-ia,L,S,{|/Ja;.L;s;.), 


(C.8) 


where  the  notation  above  is  defined  and  described  in  more  detail  in  section  9.5 
in  Cowan.  The  particular  examples  we  will  explore  will  have  no  more  than  two 
electrons  in  the  relevant  subshells  i  and  j. 

The  factor  D3  arises  from  the  uncoupling  of  the  total  spin  Sq  and  the  total 
orbital  angular  momentum     from  the  total  angular  momentum  J'g,  using  eq.  C.2, 


D3   =  55„5^(-l)^'+'^'+-^''+'[JV^,]^/M 


kJ q 


>  . 


(C.9) 


This  expression  will  be  applied  to  cases  where  the  matrix  elements  act  between 
states  with  the  same  configuration  and  those  that  do  not. 

The  factor  is  a  product  of  terms  arising  from  the  successive  uncoupling  of 
orbital  angular  momenta  for  subshells  included  following  the  subshells  relevant  to 
the  transition. 


A     =       n  '^a.i.S.,<L'„S'^55.,5;„(-l)^'"-+^-+^-+' 


a     t  c 


(C.IO) 
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Like  D3,  this  is  applied  to  cases  where  the  matrix  elements  act  between  states 
with  the  same  configuration  and  also  between  those  that  do  not  have  the  same 
configuration. 

The  factor  D5  arises  from  recoupling  the  angular  momenta  and  spin  of  the 
lower  subshell  relevant  to  the  transition  L,  with  the  intermediate  total  orbital 
angular  momentum,  including  the  next  lower  subshell,  £i_i  after  the  last  electron 
coordinate  has  been  uncoupled, 


1     r  T 

''I 


>  X  spms. 


(C.ll) 


The  factor  arises  if  the  relevant  subshells  of  the  transition  are  not  adjacent 
in  the  series  as  the  angular  momenta  are  coupled,  and  it  accounts  for  "jumping" 
over  the  intervening  subshells.  The  re-ordering  the  the  subshells  involves  using 
a  product  of  Wigner  6-j  symbols  for  each  intervening  subshell  for  both  orbital 
angular  momentum  and  spin,  as  shown  below: 


De   =     n  6^ 


I     c  r 


*  X  spins. 


(C.12) 


The  factor  D7  includes  the  reduced  matrix  element  between  the  relevant  sub- 
shells  as  well  as  terms  to  recouple  the  orbital  angular  momenta  of  those  subshells 
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to  the  total  orbital  angular  momenta  and  to  recouple  the  spins,  so  that 


Dy  = 


Sj     Sj-i  s 


(CIS) 


> .  (C.14) 


APPENDIX  D 
EVALUATION  OF  ANGULAR  PORTION  OF  B^^l 


In  this  appendix  we  will  simplify  the  angular  portion  of  B^^\  using  symmetry 
properties  to  reduce  the  factors  outlined  in  Appendix  C.  The  example  of  the  first 
order  shift,  B^^^  is  an  excellent  introduction  to  these  methods.  Please  refer  to 
Appendix  C  for  a  description  of  all  notation  as  well  as  eqs.  C.l,  C.2,  and  C.2. 
From  Section  3.4,  we  learn  that  the  trace  over  perturber  electron  states  allows  a 
non-zero  result  only  for  the  monopole  term,  ^  =  0,  of  the  multipole  expansion. 
This  simplifies  the  calculation  considerably.  We  return  to  the  examination  of  5'^) 
from  eq.  3.49  for  the  case  when  states  a  and  a'  are  in  the  same  configuration. 
Later  in  this  appendix  we  will  refer  to  eq.  3.52  for  cases  where  they  are  not  in  the 
same  configuration.  Let  us  examine  eq.  3.49, 


I    s  =  l 


'[J]  1/2 


'^Uj)D{s){nMwsA'\s)C^'\s)\M^^^  (D.l) 

In  the  case  where  the  two  states  of  the  matrix  elements  are  in  the  same  config- 
uration, we  must  first  uncouple  the  total  spins,  S  and  S',  and  total  orbital  angular 
momenta,  C  and  from  the  total  angular  momentum,  J.  This  can  be  accom- 
plished using  eq.  C.2,  which  is  just  D^.  We  can  employ  eq.  A.19  when  an  argument 
of  a  6-j  symbol  is  zero.  We  must  remember  that  a  6-j  symbol  is  unchanged  when 
any  two  columns  are  interchanged,  as  well  as  if  any  two  numbers  in  the  bottom 
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row  are  interchanged  with  the  corresponding  two  numbers  in  the  top  row,  so  that 


D.  = 


3a 


=  ^5„5i(-l)^'+^'+*^''+*[JVJ,]^/M 

J',  C[  0 


(D.2) 


We  know  J,  as  simply  J,  where  q  represents  the  total  number  of  subshells 
summed  over.  Also,  the  three  arguments  in  the  top  row  must  satisfy  a  triangle 
relation  and  their  sum  must  equal  an  integer.  Therefore  the  factor,  (-i)2'C,+25,+2j, 
drops  out.  We  are  left  with  the  following: 


q^^q-'-q'-q-JqJq  jl/2' 


^S„  S'  ^CC  ^77' 
'^qi'-'a    *'«'-o  'JQ'Jc 


(D.3) 


We  continue  to  uncouple  the  orbital  angular  momenta  of  subshells  until  arriving 
at  subshell  s.  These  steps  are  contained  in  D4.  We  will  factor  out  the  first  term  in 
the  product  and  simplify  using  the  same  methods  as  in  Dj.  Then,  we  will  continue 
for  the  entire  product  series,  as  shown  below: 
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C'q        0  Cq_, 


9-1 

X   n  ^a.L.s.,a;„L'^s;„<J5.,5;.(-l)^"-^^'"+^-+' 


^m— 1     -^m  •^m 


9-1 


x^r,-..;_.(-i)"-^^''+"iA-i,A]-^/^  n  (  ) 

m=s+l 
[£  9-1 

1^9-1J  m=s+l 


(D.4) 


Thus  far,  the  factors  resulting  from  uncoupling  the  angular  momenta  have 
reduced  to  a  simple  expression.  We  have  uncoupled  the  reduced  matrix  element 
up  to  subshell  s,  which  means  we  have  the  following  expression: 


m=5+l 
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which  can  be  uncoupled  further  by  successively  applying  eq.  C.3, 


(a._i£,_ia,L,£,p'°(.)C7W(5)||«l_,4_,a;L;£;) 


>  X 


t    L'^  Is 


>K/,p'°(s)C(°)(s)||n,g.  (D.6) 


Setting  f  =  0,  we  find  that  these  expressions  reduce  very  much  like  D4.  The 
6-j  symbols  and  their  accompanying  factors  drop  out  leaving  Ls  =  L'^  and  a  factor 

of  [CsV^VU/', 


^S, ,S'         ^s.  s'  5 .  r' 


1/2 


s.,s'/L.,Li^pj7^K/,P'°(s)C(°)(s)||n,/,).  (D.7) 
By  combining  this  factor  with  factors  from  D3  and  D4,  we  get  the  following  result: 

[7-11/2  , 

-  n  ^a.L.S„.a'„.;„S'.  •  -  •  (D.8) 

m=l 

[711/2 

D{s)   =  -  •     ■  (D.9) 


Besides  reducing  to  a  simple  factor,  we  have  shown  a  and  a'  must  be  the 
same  state  as  well  as  the  same  configuration.  One  can  now  return  this  value  to 
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Section  3.4.  We  will  now  consider  the  off-diagonal  elements  of  5^^).  Di  cannot  be 
simplified  and  will  remain  as  it  is  in  eq  C.6.  If  we  make  the  assumption  that  there 
are  no  more  than  two  electrons  within  a  subshell,  D2  reduces  to  one.  As  with  the 
diagonal  cases  the  product  of  D3  •  D4  reduces  to  the  values  obtained  in  eqs.  D.3 
and  D.4.  However,  D5  and  Dg  do  not  offer  any  simplifications.  We  must  consider 
some  special  cases.  We  will  consider  three  cases  with  the  following  notation 


■  ■ 

l°  -- 

1  'm'j 

where  n  is  no  more  than  two. 


D.l    Case  1: 

Examining  D5  for  the  first  case,  we  become  aware  that  L'^,  the  total  angular 
momentum  of  the  ith  subshell  for  state  a'  may  only  be  li  or  zero,  depending  whether 
the  occupancy  is  one  or  zero.  It  must  be  remembered  that  this  is  also  true  for  the 
corresponding  spin  values.  As  a  result,  D5  can  have  one  of  the  following  values: 

Si  =  l;  =  0, 


■  i-l 


\r'  T  11/2  J  ^'-^  0 
i^i)        \  }  X  spms. 


(D.IO) 
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which  we  may  simplify  using  eq.  A.  19, 


■  i-l 


li      0  £' 


►  X  spms 


i-i 


(-l)^^-'^'^^'-[>C;,/.]V25^,_,^;(-l)^.--^^-+'-[£,_i,/,]-V2    JJ  ...    X  spins 


(n--)^r._,z:;(-l) 

\m=l  / 


X  spms 


\m=l 


or  L;  =  ii,s;  =  s, 


■  i-l 


>  X  spms 


(D.ll) 


(D.12) 


Therefore,  D5  simplifies  for  the  case  resembling  a  hydrogenic  model,  but  not  for 
more  complex  cases.  When  considering  Dq,  we  notice  1^  =  0  for  any  intervening 
subshells.  As  a  result,  Dq  reduces  to  one.  We  now  turn  to  Dj,  we  remember  that 
we  have  already  discovered  from  Section  3.4  that  k  =  Ij.  we  may  use  eq.  A.22  to 
simplify  D7, 
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X  f - 1 ") +2'i  +l;       _ j 


(D.13) 


We  know  Cj  must  equal  C'j  so  that  the  triangle  relations  of  that  row  are  satisfied. 
Although  this  is  an  improvement,  we  must  examine  the  other  terms  to  determine  if 
we  can  simplify  the  expression  further.  Dj  may  be  simplified  because  Lj  =  Sj  =  0, 
and  L'j  =  =  /j,  with  corresponding  values  for  spin.  We  also  remember  that 
i  -  1  =  z 


0     «Sj_i  s 

=  <^5,,5;<^/.,./£,£;(-i)^-^^^^[5,_i,s;]^/25s;,<j5,5,_,(-i)^-'^^^^[s;,^^ 
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which  can  be  reduced  further. 


=    <^5,,5j5s^.,<55,5^_//,,/,^£,.£;.^£,£,.-p^.  (D.15) 

We  may  now  multiply  the  factors  together  to  get  the  following  expression, 

s;  =  l;  =  0, 

m—j+l 

(D.16) 


[/,]V2 


L;  =  ii,s;  =  s, 


I       li        C,i  Li 

xspins  X  (^5.5.  •  •  • .  (D.17) 


D.2    Case  2: 


With  respect  to  case  2,  we  may  use  our  results  for  D5  when  L-  =  0  in  eq.  D.ll, 
which  is  very  simple,  and  the  same  results  for  Dg,  that  the  factor  equals  one.  We 
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now  remember  Lj  =  Ij  -  U  and  Sj  =  s  and  consider  D7, 


D7  =  6s„s'^{-^f',-^^^.^^'\S,_„^^^'^ . 


(D.18) 


X 


L'j      Lj  Cj 


S'j_i    s  Sj^i 


>  . 


(D.19) 


(D.20) 


We  have  arranged  the  argument  order  so  that  the  two  6-j  symbols  resemble 
each  other  and  the  result  for  D5  from  Case  1.  In  fact,  result  for  the  entire  product 
will  be  very  similar, 


•  X  spins  X  5c- c'.  •  ■ 


(D.21) 


D.3    Case  3: 


In  the  third  and  final  case,  we  find  D5  and  D^  have  their  simplest  solutions 
from  eqs.  D.ll  and  D.15,  although  we  cannot  say  j  -  I  =  i,  rather  j  -  1  =  m. 
However,  we  must  consider  Dg,  using  the  fact  that       ^  Im,  Cm-x  =  A,  and 
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^m-i  =     ^-nd  the  corresponding  values  for  spin, 

^6    =    (^a   L   S    a'  V  S'  (-l)''+''"+^'+^'- 


'■m     >~'m  '~'m 


>  X  spins. 


Combining  our  results,  we  obtain  the  following  expression: 


(D.22) 


("•A:)(-l)^''t^(-l)^-+'-^''"+^''"[£„C]V2j  ^' 


X  spins  X 


(D.23) 


As  one  can  see,  the  transitions  with  doubly  occupied  manifolds  are  more  com- 
plicated that  those  with  no  more  than  singly  occupied  subshells,  but  these  three 
cases  can  be  reduced  to  forms  much  simpler  than  the  original.  Note  that  all  the 
subshells  other  than  j  and  i  have  the  same  quantum  numbers  (a^,  L^,  and  S^) 
for  states  a  and  a'.  This  confirms  our  conjecture  that  since  k  =  Ij  the  off-diagon 
terms  must  be  non-zero  only  between  states  within  a  Rydberg  series.  Many  of  the 
solutions  we  have  obtained  for  D5  and  Dq  will  be  applied  to  the  second  order  terms 
as  well. 


APPENDIX  E 

EVALUATION  OF  ANGULAR  PORTION  OF  £>/,/",/' (i^i, /iTs) 

In  this  appendix,  we  will  evaluate  the  angular  matrix  elements  of  the  operator 
second  order  in  the  radiator-perturbing  electron  potential,  A,i",i'(^i>  ^^2)-  Please 
refer  to  Appendix  C  for  a  description  of  all  notation.  In  Figure  E,  a  schematic 
has  been  drawn  to  show  the  cases  one  must  sum  over  to  calculate  the  second 
order  term  for  the  case  where  no  more  than  two  electrons  are  allowed  to  interact. 
There  may  be  other  electrons  in  other  manifolds,  but  we  will  only  consider  those 
cases  where  the  contributions  from  other  manifolds  are  minimal.  The  a"  states 
contained  in  the  rectangular  boxes  are  of  the  same  configuration  as  the  states  a 
and  a'  while  not  necessarily  identical.  The  terms  contained  in  rounded  boxes  are 
"near-hydrogenic,"  in  the  sense  that  an  electron  is  undergoing  a  transition  from 
a  singly  occupied  subshell  to  an  empty  subshell  without  crossing  any  occupied 
subshells. 
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In  the  case  where  the  two  states  of  the  matrix  elements  are  in  the  same  con- 
figuration, we  may  refer  to  Appendix  D  for  guidance.  We  have  the  same  factors; 
although  we  cannot  set  the  multipole  index  t  to  zero.  However,  we  may  take 
advantage  of  sum  rules  to  simplify  our  expression. 

E.l    Single  Configuration,  Non- Transitioning  Matrix  Elements 

In  this  section,  we  will  simplify  the  expression  for  A,i",i'(A;i,  ita)  in  eq.  3.69  first 
by  reexpressing  the  sum  over  individual  electrons  to  a  sum  over  subshells  using 
eq.  3.39.  We  define  the  quantity,  A'(*)  •  A^^J,  as  follows, 

a"  s=l r=l 

{'ifa',j\\wrA'^'\r)C^'Hrm,„,j„y.  (E.l) 

This  expression  can  be  reduced  to  a  collection  of  sums  dependent  on  whether 
s  is  higher,  lower,  or  equal  to  r,  as  shown  the  the  expression  below. 

r=l s<r 

+ 1  Dir,  r)wnA'^l,^^,r  +  E  E  D{r,  s)w.w^Af^,^^,A'^l^,^  ,(E.2) 

'■=1  r=l  s>r 

where 

A'nM   =   M\A'\s)C^'\s)\\nJs).  (E.3) 
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In  order  to  obtain  the  factors  D{r,  r)  and  D{s,  r),  we  must  uncouple  the  angular 
momenta  as  we  did  in  Appendix  D.  We  will  continue  to  make  our  simplifying 
assumption  that  there  are  no  more  than  two  electrons  within  a  subshell.  We  now 
uncouple  the  total  spins  and  total  orbital  angular  momenta  from  the  total  angular 
momentum,  J  by  reducing  the  product  D3  ■  Ds,  where  D3  is  defined  in  eq.  C.9, 


^q  \J q 

Jq  t 


x*5i',5'(-i)-^'-^^'^''-lj-;,j;r/M 


^'q     ^'q  Jq 

Jq   t  q 


>  . 


(E.4) 


We  know  J,  as  simply  J,  where  q  represents  the  total  number  of  subshells  summed 
over.  We  will  employ  the  completeness  relation  written  in  eq.  A.  17  that  takes  ad- 
vantage of  the  fact  that  5,  =  S'^  =  S'^  and  J  =  J'.  We  also  use  the  symmetry 
relations  that  a  6-j  symbol  is  unchanged  if  any  two  columns  are  interchanged,  as 
well  as  if  any  two  numbers  in  the  bottom  row  are  interchanged  with  the  corre- 
sponding two  numbers  in  the  top  row. 


Ks:,'5s>>,s'(-lf^-^'^[Jq,J''] 


JH 

^'q 

7" 

Oq 

>  < 

Jq  t 

Jq  t 

^q  J 

'-'qt'-'q      '-'q  t'^q  '^q,'~q 


[Cql 


(E.5) 
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This  simplifies  well,  and  we  will  find  that  D4  ■  D4  simplifies  well  also,  until  we 
must  treat  the  relevant  subshell.  We  write  D4  ■  D4  as  follows: 

9  9 

m=s+l n—r+\ 


£11        f  nil 


>  < 


ri  J I  ri 

nil        f  rii 


(E.6) 


Taking  advantage  of  the  Kronecker  deltas  from  ■  D^,  we  factor  out  the  terms 
related  to  the  gth  subshell  and  take  the  sum  over 

D,D,  =  5.,,,s„ai'L^'s^'<55„5i'(-l)^'-^^'+^'^'[£„£;f 


Lq 

^9 

II 

C" 

C" 

t 

Vi  ^ 

t 

c 

9-1 

/  9-1 

n 

I  m=«+l 


n 

I  m=r+l 


(E.7) 


We  can  take  advantage  of  the  relations  L,  =  L'^  =  L'^  and  C'^  =  to  help  us 
simplify,  using  eq.  A. 17, 


CJi 


6r    ,  r> 


9-1 

'  [^9-1]  Vm=s+1 


9-1 

n 

m=r+l 


(E.8) 
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We  can  successively  reduce  each  set  of  6-j  symbols  until  we  reach  either  subshell 
s  or  r.  For  the  moment  we  will  assume  r  is  the  higher  subshell.  First,  let  us  write 
the  reduced  D4  ■  D4  up  to  subshell  r, 


C" 


C" 


\.^r\  m=s+l 


nil  ±  nil 

n  <^«nL„s„,<L'^s'^<5Q;;L'^s'^,a'„L'„s;,(^/;„,/:;.<55„,5;'*s;',5;-  (E.9) 


9 

E 

n—r+l 


If  r  =  s,  the  product  series  drops  out  and  we  may  examine  the  remaining  matrix 
elements.  If  r  7^  s,  we  must  retain  this  product  series  as  a  factor  of  D{s,r).  We 
can  find  some  symmetries,  although  the  expression  seems  difficult  to  simplify.  For 
example,  we  will  show  D{s,r)  =  D{r,s).  Before  we  can  show  this,  let  us  examine 
the  factor  that  arises  from  uncoupling  the  angular  momenta  of  the  last  electrons 
in  subshells  r  and  s  respectively, 

D{s,  r)  =  D3-  D3D4  ■  A^s.,s^'55;',5^(^<.,_i£,_,,a;'_,£;'_,^a;'_,£';_j,a;_,£;_^ 
<^s„si'%',s;  E 

L  V 

(_l)^,-fL+L.+t(_  j)-/.-L'-L;-tj^^^  C'X'-'IL,,  L':]"^[L';,  L',]"^ 


L  h 

Ls 

L' 

Ir      L[  1 

>  < 

>  < 

t 

J 

t 

J 

Is 

4 

t 

L'l      Ir  ] 

ar 


where  the  sums  over  L  and  L'  represent  the  the  intermediate  sums  over  angul 
momenta  before  that  of  the  last  electrons  within  subshells  s  and  r  respectively  are 
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included.  In  this  work,  we  have  chosen  to  consider  cases  where  there  are  no  more 
than  two  electrons  in  a  given  subshell,  therefore,  L  =  L'  =  0  if  the  subshell  is 
singly  occupied  or  L  =  and  L'  =  1^  if  doubly  occupied,  so  the  sum  reduces  to  a 
single  term. 

We  can  consider  two  cases,  t  —  s  and  r  ^  s.  If  r  =  s,  we  may  sum  over 
L'l  =  L';  and  =  C;  Also,  from  eq.  E.9,  we  know  =  £;  and  =  C'^.^  since 
all  subshells  below  r  are  in  the  same  state.  We  express  D{r,  r)  as  follows: 


D{r,r)   =  I>3-/)3/)4-i^4EE(-l)^^"'''l>Cr,4'][^r,Z';T/'[L';,L;]i/2^ 


t    c;  L'; 


>  < 


'  T' 

r-l 

L  If 

> 

L 

lr 

L'r 

► 

>  < 

►  < 

t  c; 

K 

J 

t  v; 

lr 

* 

t 

K 

[Cr]  [lr. 


D{r,  r) 


[lr] 


5nn'- 


L  lr 

K 

L 

> 

< 

>  < 

t  Lf 

t 

(E.ll) 
(E.12) 


We  see  that  D{r,  r)  reduces  to  simple  terms,  we  now  examine  D{s,  r)  The  next 
step  is  to  show  that  D{s,r)  =  D{r,s).  We  cannot  simplify  much  but  we  may 
reduce  two  6-j  symbols  to  one  symbol  using  eq.  A.  18.  We  first  write  down  the  6-j 
symbol  dependent  on  in  eq.  E.IO  and  its  accompanying  factors  and  factor 
from  ■  Di  We  can  conclude  from  eq.  E.IO  that  a"  is  in  the  same  state  as  a' 
for  all  subshells  below  r  ,  which  we  can  use  to  simplify  the  expression,  so  that 
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-  £r_i  where  m  =  r.  Also,  from  eq.  E.9,      =  L". 


t     C  L'l 


}  < 


t  Lj-  jCy— 1 


>  < 


4-1  t  c; 


(E.13) 


As  we  substitute  this  expression  into  eq.  E.IO  we  remember  L'^  =  L'^  andC'L 
£'  when  m  is  less  than  r. 


D{s,r)  = 


[J] 

5s„si'5s'/,s;EE(-l)''^-^"'''(-l)^^^^^^'^'"'^^-'+' 

L  L' 


£s-l     -^-s  £s 

t      C  L' 


>  < 


Cr       L'j.  C'j._i 


t    L[  I, 


>  < 


L    Ij.  L'j. 


t  Lj.  Ij. 


r-l 

X     n  '^a.L.S.,<L'^S'^(55.,5;4(-l)^'"-+^-+^^+' 
m=s+l 


C-m—l  L. 


m  '"m 


t  C 


(E.14) 
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If  we  exchange  ->  L',,  Lr  -)■  L;,  and  so  forth,  we  find  D{r,s)  =  D{s,r)  as 
long  as  any  subshell  between  r  and  s  contains  no  more  than  one  electron.  This  is 
possible  because  Lm  =  L'^^  Im  or  the  subshell  orbital  angular  momentum  equals 
zero. 

Let  us  consider  two  practical  examples,  a  Li-like  satellite  of  a  He-like  transition 
nln'Vls  n'l'W  and  a  Be-like  transition  nl2l'\s'^  2l'2l"\s'^.  If  subshell  s  is  the 
Is  orbital  and  subshell  r  is  the  nl  orbital  the  result  simplifies  considerably  because 
Ls  =  L'g  =  Cg-i  =  Cg  =  0.  We  will  now  reduce  the  four  6-j  symbols  and  their 
accompanying  factors. 


< 

>  < 

) 

t  C^r—\ 

L 

> 

< 

II 

5to(-l)-^^- 

0      Cr—\  Lj- 


>  < 


L'       lr  Lt 


0    Lf  If 


[Cr] 


—     5f()Sc     ,  r'     6r   T'  r 


(E.15) 


If  we  evaluate  the  factor  arising  from  D^-  D^,  we  get  the  following  solution: 
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D{ls,nl) 


4-1  0 


''''''''' [Cr-ulr]y^lC:p~ 


Let  us  define  the  quantity  D'{s,r), 


so  that 


D'(5,r)  =  (-iy^+'^££j^D(.,r), 


£>'(ls,nO- Waa'(-1)'^+'' 


(E.16) 


(E.17) 


(E.18) 


We  introduce  D'{s,  r)  to  eliminate  the  factors  JT,]  and  (-1)''-+'^  from  the  final 
expression  for  ^2-  For  the  case  when  subshell  s  is  the  n'l'  orbital  and 

subshell  r  is  the  nl  orbital  we  may  take  advantage  of  the  relation  Cs-i  -  0  and  the 
knowledge  that  there  are  no  intervening  subshells.  The  product  series  from  ■  Dt 
reduces  to  one.  Since  the  subshells  are  singly  occupied  L  =  L'  =  0,  we  arrive  at 
the  following  expression: 
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D{nl,  n'l') 


Cr     L'  C' 


0 
< 

Ls 

II 

t 

J 

0 

< 

Is 

> 

>  < 

t 

h 

J 

t  jCfg  Ly 


r    V  I 


5aa'(-l)-^^-^'^+'l>C.] 


D'{s,r)   =  <5„„,(-l)-^'[/„/,]i/2 


(E.19) 


(E.20) 


When  t  =  0,  D'{s,r)  reduces  to  the  solution  in  eq.  E.18.  Although  the  terms 
between  the  same  configuration  can  seem  complex,  the  matrix  elements  between 
differing  configurations  reduces  to  a  very  simple  expressions  for  the  configurations 
under  consideration.  We  must  address  up  to  six  cases  to  show  this  is  true. 


E.2    Transitioning  Matrix  Elements 

For  non-diagonal  matrix  elements  of  multi-electron  radiators,  the  general  tran- 
sition arrays  are  set  up  so  that  one  may  calculate  a  matrix  element  where  the 
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moving  electron  travels  from  a  'lower'  subshell  to  a  'higher  subshell'.  The  simplest 
cases  to  calculate  are  those  where  the  relevant  electron  is  in  lower  subshells  in 
states  a  and  a',  while  it  is  in  a  higher  subshell  in  state  a",  or  vice  versa.  One 
cannot  immediately  assume  the  relevant  lower  subshells  are  the  same.  We  will  also 
address  the  reverse  case.  The  most  complex  case  is  when  either  a  or  a'  is  in  a 
lower  subshell  than  a"  while  the  other  is  in  a  higher  subshell.  This  is  relevant,  I 
believe,  when  looking  at  interactions  between  manifolds. 

We  will  examine  each  factor  in  Di---  Dj  to  see  what  simplifications  can  be 
made.  Some  can  be  simplified  considerably  without  making  all  the  assumptions 
considered.  We  will  show  six  most  general  cases,  where  four  of  the  cases  shown  in 
the  schematic  can  be  combined  into  two.  Not  including  the  case  where  a  —  a'  =  a", 
we  have  the  following: 


Case  1,  ■  •  •  lllp'.-' . . .  ^  . . .  /o/o/j  . . .  ^  . . .  /O/i /J-i 
Case  2,  •  •  •    /J  )■  •  •  •  l^'H]  >■■■  l^^l]  ■  ■  ■ 


(E.21) 


(E.22) 


Cases,  •  -.llH^tJ^  • . .  ^  . .  ./O/O/^/j  •  •  •  ^  •  •  ./O/^/^/O  . . . 


(E.23) 


and 


Case 4,  •  ■  ■  •  ^  •  •  •  •  ^  •  ■■ir'l'jl]'  ■  ■  ■ 

Cases,  •  --l^  ■■■h^-  --III'-'  ■■■In^-  --l^l^  ■■■In 
Case  6,  •  •  •  •..->...  /^/^/J/o  •  •  •  ^  .  •  •  /o^/O/j, . . . 


(E.26) 


(E.25) 


(E.24) 


where  n  and  k  are  no  greater  than  two.  There  are  two  cases  we  are  not  considering 
at  this  time, 

•  ••/r/?---^/---/r/]---^---/],/r^/;---,  (E.27) 
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and 

•  •  •  l^l'j . . .  ^  . . .       . . .    ^  . . .  /o/^fc-i/j, ....  (E.28) 

We  shall  show  that  in  general  k  =  li>  and  Ij  -  1^.  Therefore,  if  i  and  i'  are 
not  the  same  subshell,  it  must  be  a  cross  term  between  manifolds.  We  will  not 
be  considering  such  terms  in  this  work.  The  first  three  cases  above  involve  an 
electron  undergoing  a  transition  from  a  lower  subshell  to  a  higher  one,  and  back 
down,  while  the  second  three  cases  involve  an  electron  undergoing  a  transition  from 
higher  to  lower  and  back.  In  order  to  calculate  these  terms  we  must  remember  to 
take  the  complex  conjugate  of  one  of  the  matrix  elements,  depending  on  which  one 
is  reversed  relative  to  Cowan's  notation. 

E.3    Transition  Cases  1-3 

We  will  use  Di  as  it  is  written  in  eq.  C.6  and  consider  it  on  a  case  by  case  basis. 
We  will  make  the  same  assumption  that  there  are  no  more  than  two  electrons  in 
a  given  subshell  so  that  D2  reduces  to  one.  The  factor  D3  ■  D3  is  treated  in  an 
identical  manner  as  in  section  E.l  and  we  will  use  the  result  in  eq.  E.5. 

We  move  on  to  D4  and  once  again  take  advantage  of  the  sum  rules.  Since  the 
electron  undergoing  the  transition  is  in  the  higher  subshell  for  state  a",  we  can  take 
the  products  in  D4  over  the  same  subshells  above  j,  since  the  upper  subshell  j  will 
be  the  same  for  both  matrix  elements.  However,  we  must  consider  the  possibility 
that  subshells  j  and  j'  are  not  the  same  in  the  alternate  case.  We  now  examine 
the  first  case.  This  process  can  be  applied  very  much  like  the  case  in  the  previous 
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section,  so  that  we  arrive  at  the  quantity, 

In  order  to  reduce  D^-  D^,  we  must  remember  that  we  cannot  assume  that  the 
relevant  lower  subshells  U  and  /(  are  the  same  for  the  first  and  third  cases.  Let 
us  assume  that  i'  >  i,  i.e.  the  electron  undergoing  the  transition  is  in  a  higher 
subshell  in  state  a'  than  it  is  in  state  a.  We  will  factor  out  the  terms  related  to 
subshells  i  to  i', 


l^jJ  m=i+l 


^Crn,C'Ja,nL„,Srn,a'^U^S'^5a'^L'^S'^,a'^L'^5''  ■  (E.29) 


(i'-i 


ri         T  It  rn 


X  spms 


(E.30) 


At  this  time,  there  are  not  enough  common  elements  to  simplify  this  expression 
in  a  way  that  would  be  useful.  Like  D5  •  D5  he  product  ■  also  does  not  seem 
to  simplify  at  this  point.  We  will  factor  out  the  terms  in  the  product  series  that 
relate  to  subshells  between  i  and  i'. 
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^m=i+l 


r  r  rti 

*-'m  *-m 


>  X  spins 


(E.31) 


m=i'  +  l 


/     r"      r  7     f"  fi 


*-'m  •'-rn 


r  I         ri  nil 
'^m 


>  X  spins. 


(E.32) 


Although  there  may  seem  to  be  some  possible  points  of  approach  to  simplify 
Z>5  •  and  ■  D^,  we  must  examine  the  individual  cases.  We  finally  write  the 
formidable  term  D7, 
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There  are  a  few  simplifications  that  can  be  made  taking  advantage  of  the  fact 
that  I'j  =  Ij,  Sj  =  S'j  =  S'J  and  Sj  =  S'j. 
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Before  we  go  any  further,  we  will  remember  that  the  jth  subshell  contains  no 
more  than  one  electron  in  states  a  and  a';  this  is  the  equivalent  of  A;  <  2.  As  a 
result,  we  can  use  the  relations  Lj  =  L'^  and  =  S^.  We  can  sum  over  S'j  using 
eq.  A.  17  and  C'j  and  L'j  using  the  orthogonality  relation  for  9-symbols  given  in 
eq.  A.20.  However,  we  must  first  arrange  the  arguments  to  match  eq.  A.20.  The 
Wigner  9-j  symbols  are  multiplied  by  a  factor  of  (-1)^,  where  K  is  the  sum  of 
all  the  arguments,  when  exchanging  a  pair  of  columns  or  rows.  We  reduce  the 
expression  as  follows: 
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on   r"  T f 

^1  S'^j 


(E.35) 


This  is  an  extraordinary  result,  which  simplifies  our  calculations  considerably. 
We  also  have  shown  that  /j  equals  Z^/,  which  indicates  that  o;  may  be  in  the  same 
configuration  as  a'  or  be  in  configurations  within  a  Rydberg  series.  Let  us  now 
address  the  individual  cases  to  obtain  further  simplifications. 

E.3.1    Case  1 

In  the  first  case,  we  determine  Di  ■  Di,  knowing  n  =  1  and  there  are  no 
intervening  subshells  between  i  and  j. 


The  product  D2  ■  D2  reduces  to  one,  D3  ■  D3  reduces  to  the  expression  in  eq.  E.5 
and  D4  •  D4  reduces  to  eq.  E.29.  We  must  examine  D5  ■  D5  to  see  if  we  have 
any  simplifications.  We  may  simplify  by  knowing  that  neither  subshells  i  or  i' 
are  occupied  in  a"  so  we  use  the  following  relations,  L'/  =  L"  =  0  and  k  =  k. 
We  know  there  are  no  intervening  subshells  between  i  andz'  so  we  write  down  the 
following  expression: 

D-,-D^  =  ^a;L;s;,a;'L;'s7(^£;s;,£;'5|'(-l)^""'^^'(-l)~^'''-i~'^'>^ 


A-A-(-l) 


jt-i-jfc+i 


{l-k)  =  k. 


(E.36) 


[£:',L.]V2[£:'„4]^/M  ^ 


0  C'l, 


> 


(E.37) 
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which  simplifies  to  obtain  the  following, 

C^£._.£'/<5,L.(-l)^-^+^'+^"[>C,-l,L,]-l/25£,,  r>,Si^,,>, 

t'  — 1      l'        »  l' 

\m=l  / 

=  <^a;L;s;,a;a;'s;'<^£;5^£;'5;'(-l)''^'-^"'''^"+''^(^£._i£;'^/,LA 
xspins  j  JJ  •  •  •  j 

\m=l  / 

=  <^q;  l;  s;  ,07  l;'s;'  <^£;5;  xJ' s;'    _  1  r;'  (^ii  <^£;,  _  ^  r;^  (^i^,  l;, 

xspins  (n---V  (E.38) 


We  remember  that  d-i,  d,  and  =  satisfy  a  triangle  relation  so  the  quantity 
-1^^  where  K  equals  the  sum  of  these  arguments,  drops  out.  Although  all  these 
quantities  are  integers,  the  corresponding  spin  values  may  not  be. 

When  evaluating  ■  D^,  we  remember  that  there  are  no  intervening  subshells 
between  i'  and  j  so  that  the  second  product  series  reduces  to  one,  also,  ^'  is  unoc- 
cupied in  state  a,  so  Li>  =  0.  The  product  De  ■      reduces  as  follows: 


^a.,L.,s.„a;;L;;s;;(-l)''+^"-+^:'[£,-i,£::]i/2 


►  X  spms 


0     Ci>  C'l, 
'^a.,L,s,..;;L:;s:;(-l)''+^'^^;'[£.,£:;]^/2 

<5£;'/:;;<JA£.,(-l)''+^"+^'[£r,£r^/'  x  spins, 


(E.39) 
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which  simplifies  to  obtain  the  following, 

De-Dg   =  5a.,L.,Si,,a;',L;',s'.'/£'/£'.'/£i£;,  X  spins  (E.40) 
Now  we  may  combine  Di  ■  Di  through  Dj  ■  D7, 


7"  f"        r"     Q"  r'l  111 

'-g  '-j+i   Oj    i-j  .i'j 

_  1     (J/i  Li       _  ^  £J'// .,  l;,  ^a;,  L,,  Si,  ,q'.',  L'.',  S'.',  <^£'.' X  SpinS 

9 

m=j+l 

^£^,£'^^a^L^S^,a;;,L'^S'4^a'^L'„S'„,a'^L'4S'^      11  '  '  ' 


-  k  l^s 


(E.41) 


This  reduces  to  a  very  simple  result.  In  fact,  we  will  show  that  for  "near- 
hydrogenic"  cases,  where  =  0  and  n  =  A;  =  1,  it  reduces  to  the  hydrogenic 
result  of  Woltz  [13,  14]. 

E.3.2    Case  2 

As  in  the  first  case,  we  determine  Dx  •  Dx,  knowing  A;  =  1  and  there  are  no 
intervening  subshells  between  i  and  j, 


A  •  A  =  (n  •  1)  =  n. 


(E.42) 
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The  product  D2  ■  D2  reduces  to  one,  ■  D3  reduces  to  the  expression  in  eq.  E.5 
and  DrD^  still  reduces  to  eq.  E.29.  As  we  examine  D5  D5  in  eq.  E.30,  we  remember 
i  =  i'  so  the  first  product  series  in  drops  out.  Within  the  second  product  series  we 
employ  the  relation  d-i  =  C'l_^  =  C\_^  to  simplify  the  following  expression: 


li        C^i  L{ 


A-1   L'l  C'l 


/i-i 
X  spins  JJ 


\m— 1 


(E.43) 


We  seem  at  an  impasse,  however,  from  eq.  E.34,  we  have  the  relations  £_,_i  = 
and  j  —  1  = 


■  i-i 


^i.L;<Js.s'.  n 


(E.44) 


\m=l 


As  in  the  first  case  De-De  reduces  to  one,  and  D7  D7  is  obtained  from  eq.  E.34. 
Combining  our  results,  we  get  a  result  very  similar  to  the  first  case, 


EE  -  EE  E  1:d  d 

J"  C"       f"     9"  r"  I"  r'l 


(E.45) 


153 

E.3.3    Case  3 

This  case  is  very  similar  to  the  first  case,  however  there  is  an  intervening  oc- 
cupied subshell  m.  It  will  make  itself  apparent  in  the  quantity,  ■  D^.  We  first 
obtain  an  expression  for  Di  •  D\, 


Di-Dr  =  {-iy-\n-k)  =  1. 


(E.46) 


The  product  D2  •  D2  reduces  to  one,  D3  ■  D3  reduces  to  the  expression  in  eq.  E.5 
and  D4  ■  D4  still  reduces  to  eq.  E.29.  We  can  use  our  result  for  the  first  case  for 
D5  ■  D5  in  eq.  E.38.  We  now  examine  Dq-  De.  Since  subshell  m  is  singly  occupied 
we  may  use  the  relation  ^  L'^  =  L';^  =  Im-  Also,  as  in  the  first  case,  subshell 
i'  is  unoccupied  in  state  a  so  we  use  the  relation,  Lj.  =  0,  to  obtain  the  following 
expression: 


5a„L^s^,a'AL'^s^<Ja'„L'„s;„,a;;.L'^s;;(-l)^''~^-' 


>  X  spms 


X  spins. 


I     C  C" 


>  < 


li      JC-^i  JC^i 

I    a  c" 


(E.47) 


The  first  6-j  symbol  with  drop  out  as  it  did  in  Case  1,  but  we  need  to  use  the 
relation  that  £j_i  =        where  i  -  1  =  m. 
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^ai,Li,S,,,a",L",S",Sc''C",SCi,Ci  X  SpinS 


OfV„f,_';      n"  T"        On'    T>   <»'    n"  1  >>  <?" 

(-i)^.'-^'.'[£,,c]^/'[>c:mC]^/' 

1%    C^i    £,{1      I  li    d'^i  Cfj^i 
ICC"      \  1      C  C" 
^Q;,Li,Si,,Q'.',L'.',s;'/£'/£'.',<^£i,£i  X  spins 

«^a^L^s^,a'4.z,'^s;;<^a'„L'„s;„,a;;.L'^s;;.<j£,,£'  X  spins. 


►  X  spins 


(E.48) 


Combining  our  results,  using  eq.  E.34,  we  once  again  get  a  result  very  similar  to 
the  first  case, 

EE---EE   E   E^-^=W<^aa'.  (E.49) 

At  this  point,  the  author  finds  it  quite  amazing  that  result  reduces  to  something 
so  simple.  This  allows  one  to  write  a  code  for  simple  multi-electron  cases  without 
keeping  accounts  of  all  the  intermediate  angular  sums. 


E.4    Transition  Cases  4-6 

In  these  cases,  it  is  the  electron  in  the  jth  subshell  that  is  undergoing  a  transi- 
tion. We  find  the  result  for  D/^  ■  D4  somewhat  more  complex  because  we  may  not 
always  assume  ;  equals  /.  Let  us  rewrite  D4  •      with  that  understanding, 
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m=j+l  n=j'+l 


rii  T II  nil 

'-'m  *-'m 


r 


t  C 


m— 1 


nil         J II  nil 

CJ      t  c 


>  . 


(E.50) 


As  one  can  see,  the  order  is  different  as  in  the  first  three  cases,  because  we  must 
keep  track  of  the  matrix  elements  in  which  the  transitioning  electron  moves  from 
a  lower  subshell  to  a  higher  subshell,  or  vice  versa.  Once  again,  we  can  use  the 
relationship  from  eq  E.5,  =  to  reduce  the  previous  expression  in  the  same 
manner  as  in  section  E.3  up  to  subshell  /, 

[A] 


nil      Til  nil 

^3  S' 


n   <^«-i".s„.<i'^s;;<^5„,5;;.<^a'^L;;s'^,a'„L'„s'„%,5;„  • 


(E.51) 


We  cannot  simplify  this  further  until  we  examine  specific  cases. 

As  we  examine  D^  D^,  we  may  use  the  fact  that  subshell  i  and  i'  are  the  same, 
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■  t-i 


\tn=l 


r"      T  r 

h  L'l 


>  < 


h     C'l  L'l 


>  X  spins. 


(E.52) 


Here  we  can  make  a  useful  assumption,  that  the  occupancy  of  the  zth  subshell  is 
no  more  than  one  in  states  a  and  a',  which  is  equivalent  to  saying  n  =  n'  <  2.  We 
may  use  the  relations,  Li  =  =  0  if  n  =  0  and  =  =  Zi  if  n  =  1,  to  sum  over 
L"  and  S".  We  obtain  the  following  expression: 


•  i-l 


n<2. 


(E.53) 


In  order  to  simplify      •      and  Dt  Dj,  we  must  examine  individual 


cases. 


E.4.1    Case  4 


Once  again  we  calculate  Di  ■  Di, 


A  =  (-l)°(n)  =  n.  (E.54) 

he  product  D2  ■  D2  reduces  to  one,  D3  •  D3  reduces  to  the  expression  in  eq.  E.5 
id  £>5  •  D5  reduces  to  eq.  E.53.  However,  we  must  make  some  simplifications  to 
•4  •  A-  We  know  that  subshell  /  is  unoccupied  in  state  a",  so  we  may  use  the 
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relation,  £",  =  0  in  eq.  E.51, 


rn  nil 
^9"   j'  +  l 


0 

Cji   t  Cj 


n 

m=j'  +  l 


=  s. 


a,,  L,, S,,  ,a'^,  v.,  S';,  Ss-,  ,5;,  Sc,  ,C'.,  Scj  ,Cy  Sc'J ,C'^, 


X  n  (•••) 

m=j'+l 


(E.55) 


We  now  examine      ■  De,  since  L>5  •  •  •      is  completely  simplified.  We  note  that 
there  are  no  occupied  intervening  subshells  between  subshells  i  and  j  or  between 
and  /.  We  remember  that  subshell  j  is  unoccupied  in  state  a"  so  L'j  =  0, 


li+L'J+C'/+C'. 


/    r'  f" 


£J 


'  X  spms 


—   ^a'.L'.s'o!''L'fS'!Sc'.c'.Sc"c"  X  spins. 

J      J     J       J       J      J  J  ^ 


(E.56) 


Returning  to  D7  •  £)7,  this  product  does  not  simplify  as  easily  as  in  the  first 
three  cases.  We  begin  by  writing  out  the  full  expression: 
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y 

1  s 

1  h 

f 

C" 

h 

>  < 

T" 

t 

4 

C" 

t 

For  this  case  we  will  remember  that  subshells  j  and  /  are  unoccupied  in  state  a", 
so  we  may  use  the  relations,  S"  =  S",  =  L'-  =  L"-,  =  0, 


Si 

>  < 

's',_ 

1  ^3' 

Cll 
^3' 

0 

s 

0 

an 

1  s 

q 

f 

'-3 

c 

li 

0 

h 

>  < 

0 

Iji 

(E.58) 

C" 

3 

t 

nil 

^'3' 

t 
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Since  arguments  in  a  column  or  row  of  a  9-j  symbol  must  obey  a  triangle  relation, 
L'j,  =  Iji,  Lj  =  Ij,  and  =  C"  =  C",.  By  employing  eq.  A.  19  and  A.22,  we  get 
the  following  expression: 


D-j  ■  Dj    —  5s j  ,5'/ ^S'  ,5".  ^StS^S'.'S"      ,  s^S", s" 


( 

t 

3 

h 

C" 

>  < 

0 

J 

DtDj  = 


r'      I  r" 

^SjS^S''S''h'.,sh",  S'' h'.,  L,  hj U  5c'!c" ^C"L", 

J     '       ]'         j'    J        j'  J        J  J       3     1        ]  j' 


t 

q 

>  * 

C'j, 

t 

•> 

[  h 

h 

k 

We  will  employ  some  of  the  relations  discovered  previously,  Cj  =  Cf 
£=/!'  —  C 


(E.59) 


q,  and 


Dr-Dr  = 


^S3,s^^s^„s^,[q,c^^'[q,c,]'^' 


it      £i  lj 


q  t  q 


C'J 


(E.60) 
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Once  again,  the  product  Dy  ■  Dj  reduces  down  to  a  simple  solution;  however, 
there  is  a  factor  of  (2/^  +  1)  in  the  denominator,  rather  than  (2/^  +  1).  This  factor 
is  dependent  on  the  orbital  angular  momentum  of  the  electron  undergoing  the 
transition  when  it  is  in  state  a  and  a',  which  is  consistent  with  the  first  three 
cases.  We  now  obtain  an  expression  for  the  entire  product, 

EE---E^-^  =  n!^W.  (E.61) 
J"  C'      c",  I'jJ 

E.4.2    Case  5 

First,  we  calculate  Dx  ■  Di, 

D,  ■  D,  =  {~l)'-'-'+\n  ■  k)  ^  k  (E.62) 

The  product  D2  ■  D2  reduces  to  one,  D3  ■  D3  reduces  to  the  expression  in  eq.  E.5 
and  D5  •  L>5  reduces  to  eq.  E.53.  Since  j  =  f,  we  may  use  the  solution  for  the  first 
three  cases  for  A  •  D4  in  eq.  E.29.  The  product  ■  De  reduces  to  one  because 
there  are  no  occupied  subshells  between  i  and  j.  This  leaves  us  with  DtDj.  We 
can  use  the  relation  =  C'j  from  D4  •  D4,  and  d  =  £^  and  <S.  =  S',  from  D5  ■  D5 
to  obtain  the  following  expression: 
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5" 


Si   Sj  Sj 

1 

Si 

< 

S'J 

S'j  SI'  s 

Cll 

[ 

SI'  s 

> 

CI 

q  I,  ^ 

q  Cj  t  ^ 

q 

q  t 

7  — 


\qx\.q\{q.q\'i'h,^'^ 

q  t  q\{ q  t  q 

u  q  \\  Ci  li  q 

Lj  ij  q      q  /  q 


(E.63) 


We  obtain  the  same  solution  as  Case  4,  shown  in  eq.  E.61,  barring  D^  -  Di. 

EE-ED-D   =  (E.64) 


J"  c>^  c'p 


E.4.3    Case  6 

For  this  case,  the  product  D2  ■  D2  reduces  to  one,  D3  ■  D3  reduces  to  the 
expression  in  eq.  E.5  and      ■      reduces  to  eq.  E.53.  The  product      ■  reduces 
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to  the  result  for  case  4,  in  eq.  E.55  The  only  difference  between  this  and  case  4  is 
that  there  is  an  intervening  occupied  subshell.  We  first  obtain  and  expression  for 

D,  ■  Di  =  {n-k)  =  l.  (E.65) 

We  then  examine      ■  Dg  by  writing  down  the  full  expression, 


n  '^a'^L;;S'^,a.L„S„(-l)-'-~''--''--'-'^" 
m—i+l 


I     r  r" 

h      ■'-m-l  ■'-m-l 

r 


J II  nil 


►  X  spins 


/-I 


X  n  '^<lj(s;:,<l;.s;(-1)''+^"^^"-+^" 

n=t+l 


r  II        nil  ni 


»  X  spins. 


(E.66) 


We  rearrange  the  terms,  factoring  out  the  term  relevant  to  subshell  j,  and  remem- 
ber that  subshell  is  unoccupied  for  state  a",  so  we  may  use  the  relation  L"j=Q, 


/    C  C" 


0  q 


>  X  spins 


"m        "^m  '""i  i/m 

S     On"  T" 

-I  —  T  " 

-£ 

// 

m-l 

1'  ^m 

1/2 

f  /  £ 

1      ''I  '^m—l 

nil 

k 

■'-'m-l 

1 

>  X 

r 

C" 

£' 

m 
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^q"L"S'.'  a'  L'S'        C  ^C"  C"  ^  SpinS  X  5n"  T"  V  n    T  S 


X  < 


T "    r"  r 


>  < 


'^m 


►  X  spins. 


We  obtain  the  following  expression  by  summing  over  C", 


^Dq-  De   —   Sa" L"s'! ,a'.L'.s'.Sc>c'.Sc''c'!  X  spins 


C'J 


J   3  }  '  J  :  ] 


(E.67) 


<^a'^L;;s;;„a„L„S„.5a'^L'^S'^,a;„L'„S'„(^£„£'^(^5„5;„  •  (E.68) 


Since  subshells  j  and  /  contain  no  more  than  one  electron,  the  results  for  DtDj 
are  the  same  as  for  Case  4,  except  that  subshell  j  -  1  is  subshellm  as  opposed  to 
subshell  i.  However,  Cm  =  C'^  and  Sm  =  S'^,  so  we  obtain  the  same  solution  as 
eq.  E.60,  so  that  we  arrive  at  the  following  solution  for  D  ■  D, 


J"  C"  C", 


Saa'. 


(E.69) 


For  a  multi-electron  atom  with  no  more  than  two  electrons  having  the 
principal  quantum  number,  we  arrive  at  a  single  solution, 


same 


Q"  con  fig"  l'qJ 


(E.70) 


where  1^  is  the  orbital  angular  momentum  of  the  electron  undergoing  the  transition 
while  it  is  in  state  alpha.  The  sum  over  states  a"  includes  the  sum  over  all 
total  and  intermediate  angular  momentum  states,  but  in  the  process  of  reducing 
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our  expressions,  we  either  executed  those  sums  or  they  were  reduced  to  delta 
functions.  As  a  result,  we  need  only  to  sum  over  all  possible  the  configurations 
of  a",  remembering  that  the  case  when  alpha"  is  in  the  same  configuration  was 
addressed  in  section  E.l.  • 
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